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Introduction 



This workshop is part of the activities in the project Study of Strongly Interacting Matter (acronym 
HadronPhysics2), which is an integrating activity of the Seventh Framework Program of EU. This 
HP2 project contains several activities, one of them being the network PrimeNet having the 
focus on Meson Physics in Low-Energy QCD. This network is created to exchange information on 
experimental and theoretical ongoing activities on mainly rj and rf physics at different European 
accelerator facilities and institutes. 

The present workshop included the three general topics: 

1. 77, rj' and 4> decays from experimental and theoretical perspectives. 

2. Meson production in photon-induced reactions and from NN collisions as well as e + e~ pro- 
duction from pp collisions. 

3. Interaction of 77 and rf with nucleons and nuclei. 

The talks covered the very recent achievements in the respective fields from the experimental 
facilities KLOE at DAPHNE, Crystal Ball at MAMI, Crystal Ball and TAPS at ELSA, and 
WASA-at-COSY, as well as from different theory institutes. Electronic versions of the talks can 
be found on the PrimeNet homepage http://www.fz-juelich.de/ikp/primenet 

The detailed program was arranged by a program committee led by Dieter Grzonka and Susan 
Schadmand. 

The workshop was held Sept 26-28, 2011, at the Forschungszentrum Jiilich, enjoying kind hospi- 
tality and support from the Forschungszentrum. 

Financial support is gratefully acknowledged from the European Commission under the 7th Frame- 
work Programme through the 'Research Infrastructures' action of the 'Capacities' Programme; 
Call: FP7-INFRASTRUCTURES-2008-1, Grant Agreement N. 227431. 
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INTRODUCTION: RARE rj DECAYS AND 
TESTS OF FUNDAMENTAL SYMMETRIES 

ANDREAS WIRZBA* 



1.1 Introduction 

The organizers asked me to give an introductory talk for this meeting. Because of time limitations 
and also my ignorance about other subjects, I interpreted this task in a somewhat restricted sense 
by only concentrating on the physics of rare 77 decays. This is of course only a small subset of 
the various topics covered in this meeting. Many other aspects would have deserved dedicated 
introductions as well. Let me also remark that I heavily relied on Ref. [1] in preparing this talk, 
a joint paper with Andrzej Kupsc about tests of fundamental symmetries in rare 77 decays. 

After these disclaimers, I would like to start with two questions: Why is the r\ meson interesting 
in general and why are the rare r\ decays especially interesting? Well, the r\ is special, because 
of essentially three reasons: it is a Pseudo-Goldstone boson of the strong interaction part of the 
Standard Model (SM), i.e. its mass and all its interactions vanish in the chiral limit, the limit of 
vanishing quark masses and external momenta. It is a neutral meson with respect to all its internal 
quantum numbers. Finally, it is the heaviest of the octet of pseudoscalar mesons and therefore 
many decay channels are energetically accessible. The last two arguments do hold of course for its 
heavier "nonet-brother" , the rj meson as well which, however, is not a Pseudo-Goldstone boson, 
because the V a(X) symmetry is explicitly broken by the non-perturbative instanton physics. 

Although both mesons possess many open decay channels, all the decays — with the exception 
of the allowed decay rj — > i]irir — are suppressed in one way or the other: by the conservation of 
Nocther currents (in the chiral limit) which in turn may be broken by (quantum) anomalies, by 
the conservation of parity, C-parity, CP, G-parity etc. Since the 77 and 77' are eigenstates of G, 
P, C and CP, namely I G (J PC ) — + (0 h ), both mesons can serve as "test laboratories" for the 
conservation or rather the breaking of these discrete symmetries. Moreover, since there is mixing, 
namely the singlet-octet mixing of the r\ and rj mesons because of the U a(\) breaking, and - 
to a lesser degree - the mixing with the 7r° meson because of isospin breaking, there is also some 
access to these mixing and breaking phenomena as well. 

In responding to the second question one should keep in mind that the pertinent SM-mcchanisms 
of rare decays have to be small to begin with. Moreover, symmetry-violating decays are rare by 
nature, as otherwise we would not speak about a symmetry. Thus, by studying rare 77 decays, 
we have the chance to encounter a window of opportunity for experimental searches for physics 
beyond the SM. Since many different final states are energetically accessible, a variety of different 
symmetry tests is possible. These tests include searches for rare of forbidden decays, for asymme- 
tries among the decay products in not-so-rare decays, and for light dark-matter particles if there 
is a discrepancy between SM expectations and experimental results. Since the r\ and 77' mesons 
belong to the class of flavor-neutral (light) mesons, especially flavor-conserving symmetry violating 
mechanisms are tested - complementary to the flavor-changing CP tests with kaons or jB-mesons. 



1.2 CP tests 

Most of the CP tests in 77 decays are modeled in analogy to GP-violating decays for the Kl 
meson. Whereas the latter are flavor-changing, the former decays test - as mentioned above - a 
flavor-conserving CP violating mechanism. 



*Institut fur Kernphysik, Institute for Advanced Simulation and Julich Center for Hadron Physics, 
Forschungszentrum Julich, D-52425 Julich, Germany; a. wirzba(8fz-juelich.de 
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1.2.1 T] — >■ 7T7T 

As reviewed in Ref. [2], theory predicts very tiny branching ratios for the rj — > tttt decays. Since 
the flavor is conserved, these CP-violating decays have to be second order weak decays in the 
SM (governed by the CKM matrix), i.e. they are proportional to G^sm 2 ^- In addition, 
there are dynamical suppressions such that the upper bound on the branching ratio is given 
by B(CKM S m) < 2 x 1CT 27 . The contribution of the the CP-violating QCD 9 term to this decay 
is indirectly constrained by the empirical upper bound on the electric dipole moment d n of the 
neutron which in turn limits the QCD vacuum angle, such that B(9qcd) < 3 x 10~ 17 , whereas the 
'extra Higgs' models allow a slightly larger upper bound, £? (extra Higgs) < 1.2 x 10~ 15 . According 
to Gorchtcin [3] , the empirical bound on d n can be used directly to determine an upper bound for 
the branching ratio of rj — > tttt, no matter what the underlying CP- violating mechanism will turn 
out to be: for a larger value of the branching ratio than 3.5 x 10~ 14 the CP-violating part of the 
rj- loop contributions alone would already violate the bound on d n . 

The empirical numbers for the upper bounds are much larger: B(r] — > 7r+7r~) < 1.3 x 1CP 5 
from the KLOE collaboration [4] and B(rj -> 7r°7r°) < 3.5 x 1CT 4 from GAMS-4tt [5]. 

1.2.2 rj decay asymmetries 

The decay rj — > 7r + 7r~7 is governed by the CP-allowed M\ amplitude which can be traced back 
to the chiral triangle and box anomalies. Still, this decay can be used to test the CP violation 
in rj decays, if an interference with any CP-violating (direct) E\ term can be made visible. In 
principle, there might exist the interference with the Bremsradiation- modified decay rj — > ir + ir~ 
as well, but as shown in Sect. 1.2.1 this is very unlikely because of the tiny branching ratio of the 
rj — > tttt decay. Geng, Ng and Wu [6] suggested to measure the (linear) polarization state of the 
outgoing photon, since differential cross section measurements are insensitive to the interference 
patterns. In the same year, Gao [7] proposed to concentrate - in analogy to the CP experiment in 
the Kl system - on the rj — > 7r+7r~e + e~ decay and to measure rather the asymmetry between the 
pion and lepton decay planes which is sensitive to the inference term between the M\ and E\ decay 
amplitudes. Gao estimated the theoretical bound on the asymmetry as A t hoo < 0.02 by invoking an 
unconventional flavor-conserving CP-violating electric dipole mechanism (originally suggested in 
Rcf. [6]), which is especially sensitive to the ss quark content of the rj and therefore not constrained 
but other CP tests. The KLOE collaboration [8] measured A cmp = (-0.6 ± 2.5 ± 1.8) x 10~ 2 as 
empirical value which is both compatible with the theoretical bound and with zero of course. 

1.2.3 New types of CP tests without kaon-analog case 

Whereas the upper cases were modeled according to Kl decays, Nefkens [9] suggested a CP test 
involving the rj meson which has no analog in the kaon or P-meson sectors, namely the decay 
rj — > 4ir° which is characterized by a very low excess energy of 7.9 MeV. It is forbidden if all pion 
pairs are in relative S or P-wave states. In the SM, it can only proceed by P-wave or higher 
interactions which are - because of phase space limitations - extremely suppressed. For more 
information on this and other four-pion decays of the rj and rj' , please consult the talk of Bastian 
Kubis [11] in this conference. Since the decay is characterized by a unique 8 photon signal and as 
there is a very low background from direct pion production, this decay in fact is the most sensitive 
of any of the measured rj decays: B(rj — > At: ) < 6.9 • 10~ 7 [10]. The experimental limit on the 
branching ratio of the rj' — > 47r° decay, however, is larger: 5 • 10~ 4 [12]. 



1.3 Tests of the C symmetry in rare 77 decays 

Any decay of the rj (or rj') meson into neutrals including an odd number of photons in the final state 
would be a signal for C violation. The simplest example, rj — > 37 is heavily suppressed as each 
2-photon-pair has to be at least in a relative J = 2 state. The J = case would correspond to a 
radiative — > transition which is forbidden because of angular momentum conservation [13], and 
the J = 1 case would violate the Bose symmetry. The experimental upper limit for the branching 
ratio is 1.6 • 10~ 5 [14]. The next simplest case, rj — > n°j is again an example for the forbidden 
radiative — > transition. The empirical upper limit on its branching ratio is 9 • 10~ 5 from the 
Crystal Ball experiment [15]. The Crystal Ball collaboration [16] gives also the best limits on the 
branching ratios B(rj -> 7r°7r°7) < 5 • 10~ 4 and B(rj -> 3tt°7) < 6 • 10~ 5 , respectively [16]. 
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1.3.1 The decay rj — > 7r°e+e as C symmetry test 

The branching ratio for the C-allowed process r\ — > ir°-f*-f* — > 7r°e + e~, which involves two virtual 
photons, was estimated to be of the order 10~ 8 [17]. Therefore any measurement with a branching 
ratio B(r\ — > 7r°e + e~) larger than this value would be a clear signal for C violation. Note, however, 
that the decay mechanism with one off-shell photon is not only suppressed by C, but also by a 
vanishing transition form factor in the on-shell limit, since the "parent" decay r\ — » 7r°7 is - as 
mentioned above - a radiative — > transition. The current upper limit of the branching ratio 
B(i] — > 7r°e + e~) is equal to 4.5 • 10~ 5 90% CL, from an experiment performed in the 70s of the 
last century [18]. For further details, especially about the measurement of the WASA-at-COSY 
collaboration, please consult the presentation of Florian Bergmann [19]. 

The C invariance can also be tested in, e.g., r\ — > n + iT~y and in 77 — > 7r + 7r~7r° decays by ^ 
asymmetry measurements in the r\ rest frame. Most of the measurements date again back to the 
70s of the last century. For the asymmetries in the r\ — > 7r + 7r~7r° decay there exist recent upper 
bounds from the KLOE collaboration with 10 -3 sensitivity [20] . 



1.4 Dark particle searches 

Note that rare dilepton decays of r\ mesons and pseudoscalar mesons P — 7r°, r/, ?y' in general can 
be used as candidates for light dark-matter particle searches. Since the tree-level decay via a 
virtual Z boson is minuscule, the leading order decay process is P — > 7*7* —> . The latter 
is suppressed not only by two powers of the fine-structure constant a, but also by the helicity 
mismatch of the outgoing leptons, i.e. B(P — )■ l + l~)/B{P — > 77) ~ (f where mp and mi 

are the masses of the decaying meson and the leptons — e ± (or also /x* in the 77,77' cases), 
respectively. For the 7r° — > e + e~ decay, the branching ratio is of order 0(1O~ 7 ). In fact, the best 
experimental [21] and theoretical [22] values of B(ir° — > e+e~), namely (7.48 ± 0.29 sta t ± 0.25 syst ) • 
10~ 8 and (6.23 ± 0.09) • 10~ 8 , respectively, differ by three standard deviations. 

Kahn and collaborators [23] suggested as a possible explanation of this excess the tree-level 
exchange of an off-shell neutral vector boson U (pioneered by Fayet[24]) of mass mjj ~ (10 - 
100) McV. Acccording to the light dark matter models of Refs. [25, 26], this U boson mediates 
the annihilation reaction \X ~^ e + e~ of a neutral scalar dark matter particle x of (1 - 10)MeV 
mass, such that the produced positrons could account for the bright 511 keV line emanating from 
the center of the Galaxy. Kahn et al. could resolve the mismatch between experiment and theory 
of the 7T° — > e+e~ decay by assuming that the neutral vector meson U would couple to the u 
and d quark fields of the tt° and the e + e~ dilepton pair with a common axial- vector coupling 
g A = g\ - g d A = g e A of strength g A = (2.0 ± 0.5) • 10" 4 • mt//(10MeV). If the octet axial- 
vector quark-coupling is assumed to be of the same order as above, the U boson contribution to 
the branching ratio BR(?7 -4- e+e~) is estimated as <~ 10~ 9 , which is of the same order as the 
predictions of Dorokhov and collaborators [22] and much smaller than the experimental bound 
2.7 • 10~ 5 of the CELSIUS/WASA collaboration [27]. However, the branching ratio B(r/ -> n+ n~) 
is then predicted to be nearly an order of magnitude larger than the measured value B(r\ — > 
/-i + /j _ ) ~ (5.7 ± 0.9) • 10~ 6 [28], unless the axial- vector coupling of the U meson to the muon is 
smaller than g e A or the octet axial- vector quark coupling is smaller than g u A — g A or both. This 
limits of course the predictive power of the above presented U boson exchange mechanism. 

Another variant of dark matter U gauge bosons, which acquires a mass scale of 1 MeV to 
a few GeV via a 'dark Higgs mechanism' and which couples vectorially to all SM charged fields 
with a strength e <~ 10~ 3 or less, was suggested by Reece and Wang [29] (see also [30]). This type 
of [/-boson would also couple to SM weak neutral currents; however, the corresponding coupling 
constant would be suppressed by a factor of order m^/m|. Therefore, this variant cannot explain 
the missing excess of the 7r° — > e + e~ decay. First searches for U bosons of the Reece and Wang 
type are on the way, carried out by the KLOE collaboration. The relevant channels are, e.g., the 
decays (f> — > rjU, U — > e + e~ and 77 — > -yU, U -4- e + e~ . A signature for a U boson would be a peak 
in the pertinent Dalitz decay <j> — > r]e + e~ or 77 — > -ye + e~, while the background for such searches 
is the standard Dalitz decay. 

Finally note that there may be attractive alternatives for U boson searches which are not 
affected by resolving a narrow peak structure on a large conversion background: namely the C- 
forbidden 77 decays, e.g. the above mentioned decay 77 — > 7r°e + e~ would be an ideal candidate 
from the background point of view. 
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SEARCH FOR EXOTIC EFFECTS 
IN THE 77 -> e+e- DECAY 

MARCIN BERL WSKI* 
for the WASA-at-COSY Collaboration 



1.1 Introduction 

The goal the preliminary analysis is presented here the search for the very rare r\ meson decay 
r] — > e+e~ . From the experimental point of view this decay is very hard to observe due to the highly 
suppressed branching ratio (BR). In the Standard Model (SM) this fourth-order electromagnetic 
transition, predicted (see Ref. [1]) at the level of ~ 10~ 9 , is described by the leading order Feynman 
amplitude in Fig. 1.1. The branching ratio is thus suppressed with respect to the r\ — > 77 decay 
by: 

BRtheoiv -> e+e"] - BR[q — > 77] • a 2 • (m e /m n ) 2 , 

where a is the fine structure constant and the (m e /m^) 2 dependance stems from helicity con- 
servation [2]. Currently the best limit BR exp < 2.7 • 10~ 5 at CL = 90% comes from the CEL- 
SIUS /WAS A experiment [3] and it is five orders of magnitude bigger then the value expected from 
SM calculations. 




Figure 1.1: The conventional mechanism for neutral meson decay into an lepton-antilepton pair. 



1.2 The 77 — y e + e decay as a probe of new physics 

The 77 meson decay into an e+e~ pair became one of the most interesting topics in low-energy 
hadron physics, since there is a possible admixture of non-SM processes that can enhance the 
BR of this decay. The most recent and very precise experimental value of BR(tt° — > e+e~) = 
(7.49 ± 0.29 ± 0.25) • 10 -8 determined by the KTeV collaboration exceeds the latest theoretical 
calculations from Dorokhov and collaborators [4] by 3 standard deviations. The suggestion from 
Kahn et al. [5] is that the possible explanation can be the exchange of an off-shell neutral boson 
U of mass ran ~ (10 — 100)MeV. The U boson was previously proposed by Boehm et al. [6] in 
the light dark matter models to mediate the annihilation of a neutral scalar dark matter particle 
X of (1-10) MeV mass in the reaction XX ~ * e + e~. This models were used as an explanation 
for a 511 keV line emanating from the center of the Galaxy due to excess positrons produced 
in this annihilation process [5]. An observation of a signal above the level predicted in the SM 
calculations could be evidence for an unconventional process which enhances this decay rate. 

1.3 Data analysis 

The experiment was performed using the WASA detector [7] at the COoler SYnchrotron COSY, 
located in the Forschungszentrum Jiilich, Germany. The 77 mesons were produced in the reaction 
pp — > pprj at 1.4 GeV proton beam energy. We used a special trigger in which a high energy deposit 
in both sides of the electromagnetic calorimeter was demanded. This type of trigger should equally 
favor electromagnetic decays of 77 mesons like r] — > 77, 77 -4- e+e~7, 77 — > e+e - . Using the 77 — > 77 

*National Centre for Nuclear Research, Marcin.Berlowskiafuw.edu.pl 
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decay channel we estimate the number of 77 mesons collected in the 2 week data sample to be 
~ 4.4 • 10 7 . This number is approximately forty times bigger than the data sample used in the 
previous analysis [3]. It was found that the main sources of background are pp — > pp7r + ir~ due to 
a 100 time larger cross section, pp —> p(A(1232) —> p[j* — > e + e~]) having the same particles in 
the final state as 77 — > e + e~ and pp —> pp(r) —> e + e~-f) due to the same final state particles if the 
mass of the virtual photon is large and the photon is not observed. 

One of the analysis highlights, the charged particle identification by energy deposited in the 
electromagnetic calorimeter (SEC) is shown on the Fig. 1.2. 




Energy deposit in the left side of SEC Energy deposit in the left side of SEC 



Figure 1.2: Energy deposit in each side of the electromagnetic calorimeter for: MC simulated 
77 —> e + e~ signal (left panel), MC direct production of two charged pions (right panel). 



1.4 Results and conclusions 

After the full data analysis we identified 182 event candidates in the entire data sample in the 
plot of missing mass deducted from two protons before subtracting non-eta background. At this 
level the Monte Carlo signal acceptance is estimated to be 4.8%. The background coming from 
pp — » ppir + ir~ and pp — » p(A(1232) — > p[y* — > e + e~]) is subtracted by a polynomial fit on the 
missing mass plot and the expected background is coming from the 77 — > e + e~7 decay channel 
with 5.6 ± 1.1 events. The remaining number of events after background subtraction is consistent 
within errors with the expected number of events coming from 77 — > e + e~7. No events of r\ — > e + e~ 
corresponds to the BR upper limit better than the previously published one [3] by an order of 
magnitude. 

Based on 77 — > e + e~7 we have proven the ability to detect and reconstruct e + e~ pairs with 
different masses. At this stage of the analysis we are able to go an order of magnitude below the 
present limit [3] for the BR of 77 — > e + e~ decay. Four times bigger statistics of the same reaction 
with the same trigger is still available for analysis. We plan an additional eight weeks of data 
taking using the same incident proton energy and improved experimental conditions. 

This publication has been supported by the European Commission under the 7th Framework Pro- 
gramme through the 'Research Infrastructures ' action of the 'Capacities ' Programme. Call: FP7- 
INFRASTRUCTURES-2008-1, Grant Agreement N. 227431. 
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ANALYSIS OF THE C- VIOLATING DECAY 

77 — ► 7T° + 7* — > 7T° + e+ + e~ 
WITH WASA-AT-COSY 

FLORIAN BERGMANN* 
for the WASA-at-COSY Collaboration 

1.1 Introduction 

The decay 77 — > 7r° + 7* — » 7r° + c + + c~ is forbidden since it violates the C-parity conservation: 

C = C %0 ■ C 7 . = (+1) • (-1) =-1^+1 = C„. (1.1) 

So far this decay has not yet been observed, but an upper limit for the branching ratio has been 
determined in previous experiments under the assumption of a phase space distribution for the 
emission of the ejectiles: 

BR(r7 ^ vr +c+ +c~) <4-l(T 5 [l] (1.2) 

Recent high statistics measurements with WASA-at-COSY [2] aim at lowering the existing upper 
limit to test C-parity conservation with increased sensitivity. For the C-violating decay rj — > 
7T° + 7* — » 7T° + c + + c~ different decay model assumptions exist which have to be taken into 
account for the analysis, e.g.: decay via a virtual photon, vector meson dominance (see Fig. 1.1). 






+ I 




+ 



Figure 1.1: Feynman diagrams of the decay 77 — > 7r° + 7* — > 7r° + e + + e assuming different 
models. 



1.2 Analysis 

For the analysis two data sets are used. One with 10 7 registered p + d — > 3 He + 77 events from a 
beam time in 2008 and the other one with 2 • 10 7 such 77-events from a beam time in 2009. 
The first step in the analysis is a preselection on p + d — > 3 He + X events via 3 He identification 
using the AE-E method in the forward detector system which consists of a tracking device and 
an array of thin plastic scintillators. The next step is a preselection on the signature of the decay 
r) — > tt° + e + + e~ — > 7 + 7 + e + + e _ , requiring two neutral particles in the central detector 
(27), one positively and one negatively charged particle in the central detector (e + , e _ ) and one 
charged particle in the forward detector ( 3 He). In order to determine further selection conditions 
all possible background reactions like (multi-)pion production and other 77-decay channels are 
simulated. For choosing further possible selection conditions the underlying mechanisms for the 
decay r\ — > ir° + c + + e~ have to be taken into account. This is illustrated in Fig. 1.2 for the 
invariant mass distribution of the lepton pair. 

Before searching for the optimal set of cuts one has to determine the number of events from each 
background channel which are still left after the preselection. For this purpose the measured data 
(2008) has been fitted by a Monte Carlo cocktail using the same fit parameters for all differential 
spectra. This results in an accurate description of the data (See Fig. 1.3). 
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Figure 1.2: Simulated invariant mass of the lepton pair for different decay mechanisms. Left: r\ 
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simulations of decay mechanisms which, e.g., consider vector meson dominance are planned. 
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Figure 1.3: Missing mass of the measured data fitted with a Monte Carlo cocktail. 



For the determination of the optimal selection conditions the evaluation function G — Sr ■ -g^ is 
used for all cuts, where Sr is the relative signal and Br the relative background after applying 
the cuts. Only Monte Carlo data are used to find the optimal set of cuts via varying the cuts and 
maximizing the evaluation function G. 

1.3 Preliminary Results and Outlook 

The 2008 data set has been analyzed by A. Winnemoller [3]: After applying the optimal set of 
cuts 9920 out of 10 6 simulated r\ — > 7r° + 7* —> ir° + c + + e~ events remain and according to the 
simulations the background is reduced by a factor 1.7 • 10 s . In agreement with the simulations one 
event candidate remains in the measured data after applying all cuts. 

Currently further 2 • 10 7 p + d — » 3 Hc + r\ events from 2009 are analyzed considering different 
model assumptions. Additionally, significantly more p + p— >p + p + rj events are available for the 
analysis and further p + p— >p + p + X bcamtimes are planed. 

This publication has been supported by the European Commission under the 7th Framework Pro- 
gramme through the 'Research Infrastructures ' action of the 'Capacities ' Programme. Call: FP7- 
INFRASTRUCTURES-2008-1, Grant Agreement N. 227431. 
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THEORETICAL PERSPECTIVES ON 
77, 77' 4tt DECAYS 

FENG-KUN GUO* , BASTIAN KUBIS*^ AND ANDREAS WIRZBA* 



1.1 Introduction: anomalous processes 

Processes in low-energy QCD that involve an odd number of (pseudo-)Goldstone bosons (and 
possibly photons), which are, therefore, of odd intrinsic parity, are thought to be governed by the 
Wess-Zumino-Witten term [1] via chiral anomalies. While the so-called triangle- and box-anomaly 
terms are well-tested in processes such as e.g. ir°, i] — > 77 (triangle) or ~/ir — > irir, 77 — > irirj (box), 
the pentagon-anomaly remains more elusive; the simplest possible process that is usually cited 
is K + K~ — > 7T + 7r~7r°, which however has not been experimentally tested yet. A different set of 
processes involving five light pseudoscalars are the four-pion decays of r\ and rj '. Experimental 
information about these is also scarce, only upper limits on branching ratios exist [2]; however, 
this may change in the near future for at least some of the possible final states. 

In principle, the decays rj' — > Att seem not terribly forbidden by approximate symmetries: they 
arc neither isospin-forbidden, nor electromagnetic. 77 — > Air, in contrast, is essentially suppressed 
by tiny phase space: only the decay into An is kinematically allowed (M v — AM^o = 7.9 MeV, 
M v — 2(M W ± +M 7T a) = — 1.2 MeV). The fact that anomalous amplitudes always involve the totally 
antisymmetric tensor e^ap furthermore shows that no two pseudoscalars are allowed to be in a 
relative S-wave; the decays rj — > 2(7r + 7r~) and rf — > 7r + 7r~27r° are therefore P- wave-dominated. 
As furthermore Bose symmetry forbids 27r° to be in an odd partial wave, rf — > Att and r\ — > Att 
even require all ir° to be at least in relative D- waves [3]. This, combined with the tiny phase space 
available, leads to the notion of rj — > 47r° being CP-forbidden [2], although strictly-speaking it is 
only S-wave CP-forbidden. 



1.2 Charged final states: 77' — > 2(7r + 7r~), r}' — > 7r + 7r _ 27r° 

In the following, we assume the standard r/r]' mixing in terms of singlet (770) and octet (j]s) 
states, i.e. \rj) = l/3\rj ) + 2V2/3\rj s ), \rf) = 2V2/3|r7 ) - 1/3|?7 8 ). The flavor structure of the 
chiral anomaly does not allow for direct contributions of the Wess-Zumino-Witten term to the 
decays 77, 77' — > An. The leading amplitude to the two channels involving charged pions in the 
final state occurs at one loop or 0(p 6 ) in terms of chiral power counting, and is given as the 
sum of kaon loops and local counterterms, see Fig. 1.1. For instance, the decay amplitude for 
77s -> 7T+ {pi)ir~ (p2)tt + (P3)tt~ (p 4 ) is given by [4] 

A( m -+ tt+tt-tt+tt-) = ^0-^ 2P - p l[F % {s 12 ) + F s (s 34 ) - F 8 (s 14 ) - F 8 ( S23 )} , 



1 I AM 2 

Jkk(s) = ^ {l - <Jk arccot a K } , a K = y - 1 , (1.1) 

where Sjj = (pi + pj) 2 . The corresponding 770 amplitude consists of local contributions only. The 
scale dependence of the counterterm C^ 1 " is known [5] and cancels the one of the loop contribution. 
?7 8 — > 7r + (pi)n° {p2)^~ {P3)n° (pa) is given by the same amplitude Eq. (1.1) up to a sign. 

In order to estimate the finite part of the coupling constant C-J - ^ 1 * , we resort to resonance 
saturation. While this has been studied in a more general setting recently [6], we employ the 
hidden- local-symmetry (HLS) formalism to calculate vector- meson contributions, see Fig. 1.1. We 
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only need two (combinations of) coupling constants occurring in the anomalous sector [7] , c\ — c 2 
and C3. The two chiral constants contributing to five-pseudoscalar amplitudes at 0(p e ) (only one 
of them contributes to rj —> 4tt) are given in terms of these by 

Cr(M P ) = -2CE r (M p ) = i[C [ 2% Xpf ■ 

The HLS couplings are often assumed to be given roughly by c\ — C2 « C3 « 1. A more extensive 
phcnomcnological study finds c\ — c 2 ~ 1.21, C3 « 0.93 [8]. Comparing the relative importance 
of counterterm (=vector-meson) contributions to kaon loops in Eq. (1.1), we find that the decay 
amplitude is entirely dominated by the vector-meson terms. Furthermore, as the maximum allowed 
kinematical range in rj — > 4ir is yfsjj < M v > — 2M X w 680 McV, it is necessary to retain the full 
p propagators for a phenomenologically reliable description. This, at the same time, takes care of 
the essential P-wave 7T7t final-state interactions. 

We calculate the branching ratios as functions of the HLS couplings, before inserting the 
two sets of numerical values quoted above. In the isospin limit, 2 x BR(?y' — > 2(7r + 7r~)) = 
BR(r/ — > 7T + 7r~27r°); we adjust for the correct phase space by using the pion masses M x ± and 
(M„± + M 7V o)/2, respectively. We find 

BR (7/ -> 2(7r+70) = [0.15(ci - c 2 f + 0.46(ci - c 2 )c 3 + 0.37^] x 10~ 4 = {0.97, 1.04} x 10~ 4 , 

BR(t/ -> vr + 7r-27r°) = [0.34(ci - c 2 f + 1.07(c x - c 2 )c 3 + 0.85cf] x 10" 4 = {2.26, 2.44} x 10~ 4 . 

(1.3) 

These numbers are to be compared to the available experimental upper limits, BR(?7' — > 2(ir + ir~)) 
< 2.4 x 10" 4 and BR(rf -> tt+tt^tt ) < 2.6 x 10" 3 [2]. 



1.3 Neutral final states: 77' — >• 47r°, 77 — > 47r° 

As we have mentioned in the introduction, the P-wave mechanism described in the previous section, 
proceeding essentially via two p intermediate resonances, cannot contribute to the 47r° final states. 
In fact, we can show that the D-wave characteristic of r], rj — > 47r° suppresses these decays to 
0(p 10 ) in chiral power counting, that is to the level of three loops in the anomalous sector. This is 
in particular due to the flavor and isospin structure of the anomaly, which does not contain five- 
meson vertices including 2ir° at leading order (0(p 4 )), and to the chiral structure of meson-meson 
scattering amplitudes, which only allows for S- and P-waves at tree level (0(p 2 )). As a complete 
three-loop calculation would be a formidable task and is beyond the scope of our exploratory 
study, we instead construct a decay mechanism that, we believe, ought to capture at least the 
correct order of magnitude. This is built on the observation that we can calculate the leading 
imaginary part of the decay amplitude at 0(p 10 ), which is given by a 7r + 7r~ intermediate state, 
see Fig. 1.2. As indicated, we thereby insert the full vector-meson amplitude for r], rj — > 7r + Tr _ 27r° 
described in the previous section. We then reconstruct the "full" D-wave irir final-state amplitude 
by requiring that it is proportional to the Omnes function f^C 5 ) (assuming dominance of isospin 
1 = 0), which near threshold behaves like 



Imfi!j( a )«yi-— ^xf°(s) , (1.4) 

where t\ is the appropriate hit partial wave of isospin 1 = and angular momentum 1 = 2. 
Well above threshold, however, the Omnes function can well be approximated by the /2(1270) 
resonance, fl^s) rj M? /(M? — s). The resulting decay mechanism is indicated in Fig. 1.2. Note 
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Figure 1.2: D-wave mechanism for 77, 7/ — > 47r°. 



that it is far from a real / 2 dominance picture, as the p is not "short-ranged" on the scale of the 
/ 2 mass. For the branching ratios, we find 

BR(r/ -> 4tt°) = [0.4(c a - c 2 ) 2 + 1.6(ci - c 2 )c 3 + 1.6c§] X 10~ 8 = {3.6, 3.8} x 10~ 8 , 
BR(r) -> 4tt ) = [0.3(ci - c 2 ) 2 + 1.0(ci - c 2 )c 3 + 1.0c 2 .] x 10~ 30 = {2.4, 2.5} x 1CT 30 . (1.5) 

While the D-wave mechanism therefore suppresses the decay rj — > 47r° by 3-4 orders of magnitude 
compared to the charged final states, the suppression of r\ — > Air due to the combination of D-wave 
and tiny phase space is enormous, and any experimental signal for this channel could indeed be 
interpreted as a sign of CP violation. 

1.4 CP violation through the QCD 0-term 

One possible source of strong CP violation is the so-called QCD #-term, linked to the U(1)a 
anomaly, that can also be treated on an effective-Lagrangian level [9]. It induces effects such as 
an electric dipole moment of the neutron, and also allows for the decay rj — > 2n. In addition, it 
induces a CP- violating S-wave mechanism for the decay 77 — > Air [4] , 

A(v ( ^^)=- ] fl§ l x9 . (1.6) 
The resulting branching ratio is given by 

BR(t7 ^ 4tt ) = 5.2 x 1(T 5 x 0% , (1.7) 

which, were 9o a quantity of natural size, would demonstrate the enhancement of the CP-violating 
S-wave mechanism compared to the CP-conserving D-wave one, see Eq. (1.5). However, current 
limits from neutron electric dipole moment measurements suggest 0q ^ 10 -11 [10], which puts also 
this result beyond any experimental reach. 
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PLUTO: A MULTI-PURPOSE 
EVENT GENERATOR FRAMEWORK 
FOR RARE 77 DECAYS 

IN GO FROHLICH* 



1.1 Introduction 

In this contribution, the new design and recent developments of the simulation package Pluto [1] 
have been presented. Pluto's main mission is to offer a modular framework with an object-oriented 
structure, thereby making additions such as new particles, decays of resonances, new models up 
to modules for entire changes easily applicable [2]. Overall consistency is ensured by a plugin- 
and distribution manager. Particular features are the support of a modular structure for physics 
process descriptions, and the possibility to access the particle stream for in-line modifications. Ad- 
ditional configurations can be attached by the user by a dedicated script language ( "PlutoScript" ) 
without re-compiling the package, which makes Pluto extremely configurable. 

1.2 Rare 77 decays 

The rare r\ decays have been implemented into Pluto via means of a dedicated "plugin". This 
ensures that changes introduced by the new implementation do not affect simulations done by 
other users. 

The models follow basically calculations done by Wirzba and Petri [3], but leaving out tiny 
effects like magnetic form factors. 

1.2.1 77 — > e+e~e+e~ 

This decay implementation is done via two virtual photons: 

PReaction my_reaction( "_T1=2 . 2" , "p" , "d" , 

"He3 eta[dilepton [e+ e-] dilepton [e+ e-]]"); 

In the first step of the event generation Pluto samples the mass of the 2 dileptons. This is 
very similar to the implementation of the 77 Dalitz decay of the base package. The shape of the 
dilepton spectrum is based on an implementation provided by A. Kupsc [4]. 

In a second step, the angular distributions are sampled via the rejection method. This can be 
done as the correction is only about 20%, which means that only a small fraction of events have 
to be resampled. Moreover, the angular distributions of the "simple" generator are isotropic and 
therefore equally populated. 

1.2.2 77 — > ir + ir~e + e~ 

This decay implementation uses the a "genbod" technique: Here, the first step of the 4-body decay 
is done via a modified genbod algorithm, where the mass distribution M(q 2 ) = F(q 2 )/(8 ■ q 2 ) is 
folded into the genbod model (with q as the four-momentum of the virtual photon) . This is done, 
because a generation of events in a 4-body space is numerically difficult and would result in large 
computing times. On the other hand, a pure rejection method (using an unmodified genbod), as 
in the previous case, is not possible, because the pole of M(q 2 ) is located in a region where no 
events are generated by genbod. 

For the form factor by default a "simple VMD model" is used: 

F(q 2 ) = m 2 J(m 2 v q 2 ) (1.1) 
*Goethe-Univcrsitat Frankfurt, Froehlich9physik.uni-f rankf urt . de 
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with m v = 0.77GeV. 

In total, the following matrix element is used: 

A 2 = M(q 2 ) * * ft sin 2 (^) * X(m 2 rnS ^,q 2 )(l + /3 2 * sm 2 (0 e ) * sin 2 (0)) (1.2) 

with as the invariant mass of the pion pair, and the factor B x = \ 1 — ^ aL . In the case of the 
dilepton pair, s ee is the invariant mass of the virtual photon. 

1.3 Scripting with Pluto 

The idea behind scripting with Pluto is to combine objects like histograms and ntuples via a 
flexible batch language with the main event loop, model objects, and the internal data base of 
Pluto. The commands of the batch script are compiled by Pluto in advance and can be executed 
as often as required without additional delay. This execution can be done inside the event loop 
which allows for modifications of particles, e.g. to write acceptance filters or include smearing 
effects. 

Furthermore, it is also possible to call many public available Pluto-methods of the utility 
class (for random sampling) and all methods of the Pluto particle class, which is inherited from 
TLorentz Vector. As the script can combine ROOT-objects (histograms and ntuples) with the event 
loop, particle properties can be converted very quickly to histograms and ntuples. And finally, 
the script can be packed together with acceptance and resolution matrices to form a universal, 
common format for a detector description. 

1.4 Examples 

Finally, the application of the script method is discused. Exchanging, e.g., the widely discussed ui 
form factor is a quite simple action 1 : 

PSimpleVMDFF *ff = new PSimpleVMDFF ("vmd_ff_dd@w_to_dilepton_piO/formf actor" , 

"VMD form f actor", -1); 
ff->AddEquation("_ff2 = 0.17918225/( (0.4225- _q2) * (0 . 4225- _q2) + 0.000676)"); 
makeDistributionManager () ->Add(f f ) ; 

The basic syntax is based on the ROOT TFormula class, but with a lot of extensions. For 
details, reference is made to the full manual which is available at the Pluto web page [1]. 

Considering, that the form factor model of rate eta decays is a separated model, the form 
factor can be read via: 

PSimpleVMDFF *ff = (PSimpleVMDFF*)makeDistributionManager()-> 
GetDistribution ( " eta_ee_pipi_vmd_f f " ) ; 

and therefore exchanged via similar means 
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DALITZ PLOT ANALYSIS 
FOR rj tt+tt tt AT KLOE 

LI CALDEIRA BALKESTAHL* 
on behalf of the KLOE-2 Collaboration 



1.1 Theoretical Motivation 

The decay r\ — > ir + ir~ir a is an isospin- violating process sensitive to the light quark mass double 
ratio: 



Q 2 = 



-, where m = —(ma + m u ) 



(1.1) 



The rj — > 7r + 7r~7r° decay rate is proportional to Q~ 4 and it has been calculated at leading 
order (LO) and next-to-leading order (NLO) in xPT. The values obtained, Flo ~ 70 eV and 
Fnlo — 160 ± 50 eV , are significantly lower than the experimental value T exp — 296 ± 16 eV , 
obtained from a fit to all the available data [1]. The difference between LO and NLO calculations 
and the discrepancy with data at NLO are evidence of the slow convergence of the xPT series, 
due to strong pion rescattering effects in the final state. These effects can be treated by means 
of dispersion relations as recently proposed in Ref. [2]. The relations contain four parameters to 
be fixed either matching the LO xPT results or by a fit to experimental data. Using the first 
approach G. Colangelo and coauthors [2] obtained Q = 22.3 ± 0.4, in good agreement with most 
of the results presented in Fig. 1.1, where the gray elliptic band corresponds to the 2-er range of 
the Q-value above. Other plotted constraints come from lattice calculations. 
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Figure 1.1: Ratio of the light quark masses. The gray band corresponds to Q = 22.3 ± 0.8 [3]. 
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1.2 Data analysis 

The KLOE detector, at the DA$NE collider, consists of a large volume drift chamber and 
a lead-scintillating fibers calorimeter. The detector usually takes data at the <^-meson peak 
(s/s = 1.019 GeV). The experiment took data in years 2001-2002 and 2004-2006, integrating 
a total luminosity 2.5 fb _1 at the (/>-peak and 240 pb _1 at the center-of-mass energy of 1 GeV. 
KLOE-2 has been approved to collect ~ 20 fb _1 in 3-4 years of data taking at the upgraded 
DA$NE facility, which is expected to operate with a factor ~ 3 improvement in luminosity. 

The KLOE collaboration published the Dalitz plot analysis of the 77 — > 7r + 7r~7r° decay in year 
2008 [4]. In the paper, based on ~ 500 pb _1 of integrated luminosity, the Dalitz plot density 
shown in Fig. 1.2 was fitted by: 

\A(X,Y)\ 2 ~ 1 + aY + bY 2 + cX + dX 2 + eXY + fY 3 + gX 3 + hX 2 Y + IXY 2 (1.2) 

obtaining: 

a = -1.090 ± Q.005(stat)toi1l(sys) b = 0.124 ± 0.006(siai) ± 0.010(sys) 

d = 0.057 ± 0.006(stat)t^° %(sys) / = °- 14 ± 0.01 (stat) ± Q.02(sys) 

and the other parameters consitent with zero. The c and e parameters are expected to be zero 
from C-invariance. 




Figure 1.2: Dalitz plot density measured by KLOE[4] on the basis of ~ 500 pb 1 integrated 
luminosity. 

The results have been recently used by S. Lanz and coworkers to fix the parameters of the 
dispersion relations in the framework of the analysis presented in Ref. [2]. The resulting Q-value 
is Q = 21.31 ± 0.60 [5] and the Dalitz plot density obtained is in better agreement with the 
measurement performed on the rj — > 7r°7r 7r° sample [6] than previous calculations based on the 
matching of the dispersion-relation parameters with one-loop x?T results. The KLOE results 
have also been used in analytic dispersive analysis [7] for the value of R (R = )> both in a 

direct fit to the results, and in a correction to xPT NNLO low energy constants. 

More data on rj — > 7r + 7r~7r are needed to understand the origin of the residual tension between 
data and theoretical calculations, concerning both, the comparison between neutral and charged 
channel, and the Q- value. 

On the experimental side, a new analysis of the r\ — > 7r + 7r~7r° sample (Fig. 1.3) is in progress 
to improve on the statistical sample and to overcome some limitations of the previous analysis. 
For the latter, the selection efficiency is measured directly from minimum bias events, while in 
the old analysis it was evaluated from Monte Carlo simulation, giving a large contribution to the 
systematics. Further improvement on the precision of the Dalitz plot distribution is expected: 
i) replacing the kinematical fit, requiring 4-momentum conservation and the speed of light for 
photons, with constraints from the two-body <j) decay kinematics; ii) eliminating the request of 
vertex reconstruction and the cut on the total number of tracks in drift chamber; Hi) introducing 
the particle identification via time-of-flight measurement with the calorimeter; iv) performing the 
analysis on an independent, larger data sample. 
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Figure 1.3: Energy distribution of the photons from <fi 777 — » ir + ir~jjj. 
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1.1 Introduction 

The decay 77 — > 7r + 7r~7r is important since it allows access to light quark mass ratios. Since 
electromagnetic corrections are small [1, 2] this decay is a probe for the strong isospin violation. 
The calculated tree level decay width, r iree ss 70 eV [3] , deviates more than a factor of four from 
the PDG value Y = 296 ±16 [4]. Higher order diagrams are needed which take into account tttt 
final state interactions [5, 6, 7]. Alternatively, tttt final state re-scattering can be implemented in 
dispersive approaches [8, 9, 10, 11, 12] as well as other approaches [13]. The decay width scales as 

r=f%Vf, (i.i) 



Q 

where Q 2 = (ml — rh 2 )/(m 2 d — m 2 ), m = \{m u + mj), and the decay width f and Qd = 24.2 are 
calculated in the Dashen limit [14]. Q gives access to light quark mass ratios. The ratios serve also 
as an important input for lattice QCD [15]. To derive Q, Y has to be known reliably from theory 
which can be tested by comparing with the experimental Dalitz plot distributions. As Dalitz plot 
variables 

x = ytf±^,y = ^-l. (1.2) 

are used. Here T + , T_ and T denote the kinetic energies of tt + , n~ and n° in the 77 rest frame, 
and Q v is the sum of the three kinetic energies. The standard way to parametrize the Dalitz plot 
density is a polynomial expansion around x = y = 

\A{x, y)\ 2 ex 1 + ay + by 2 + dx 2 + fy 3 + gx 2 y + ... (1.3) 

where a,b,...,g are called the Dalitz Plot parameters. The most precise experimental result is 
based on a Dalitz plot containing 1.34-10 6 events obtained by KLOE [16]. Parameters b and 
/ in [16] deviate significantly from ChPT predictions. Independent measurements are therefore 
important and WASA-at-COSY [17] aims at providing two independent data sets with 77 produced 
in pp and pd reactions. 



1.2 Experiment 

In 2008 and 2009 WASA-at-COSY collected 1 • 10 7 and 2- 10 7 77 mesons respectively, in the reaction 
pd — > 3 HeA at a kinetic beam energy of 1 GeV. The missing mass with respect to 3 Hc is used to 
tag the 77 in the reaction pd — > 3 HeA. To obtain 77 -4- 7r + 7r~7r° candidates, at least two particles 
of opposite charge and at least two photons with an invariant mass close to tt° are required. To 
further reduce time coincidental events, cuts on vertices and timing requirements are included. 
To reduce background channels, mainly coming from pd — > 3 Hc7T7r reactions, conditions on the 
missing mass calculated for 3 He7r + 7r~, 3 He7r° and ir + ir~ are applied. The preliminary analysis 
gives 260 000 7/ — > n + ir~ir candidates from the 2008 data sample. The majority of the events 
populating the Dalitz plot is 3 He7r + 7r _ 7r° and 77 — > tt + tt~tt . The r\ content in each Dalitz plot 
bin is obtained by performing a four-degree polynomial fit over the direct 37r background and the 
fitted polynomial is subtracted from the signal region. The number of 77 events in each bin is 
corrected for acceptance. 

"Institute for Physics and Astronomy, Uppsala University, Uppsala, 751 20, Sweden, 
patrik . adlarsonOf ysast .uu.se 
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1.3 Result 

The preliminary acceptance corrected experimental results with statistical errors for the x, y 
projections of the Dalitz plot are compared in Fig. 1.1 to the distributions of the r\ — > 7r + 7r~7r° 
with Monte Carlo weighted with the set of Dalitz plot parameters calculated for the leading order 
(LO) in ChPT (solid blue line), next-to- next-to-leading order (NNLO) (dashed-dotted green line) 
and the central values obtained by KLOE (dashed red line). The Monte Carlo results have been 
normalized to the WASA experimental data. From the preliminary Dalitz plot projections the 
WAS A data points are in reasonable agreement with the KLOE experimental result. The work 
on the pd data will be continued in order to obtain the Dalitz plot parameters for rj — > Tr + ir~ir°. 
An effort will be made to understand the systematical effects. 



3.5 



2.5 



K 1-5 



u 1 



0.5 



x10 



.x10 



WASA 
- LO parameters 
-NNLO parameters 
-KLOE parameters 



-0.5 



0.5 




Figure 1.1: Acceptance corrected data points with statistical errors, projected on x (left) and y 
(right), compared to Monte Carlo weighted with the Dalitz plot parameters for the leading order 
calculation (LO) in ChPT (solid blue line), next-to- next-to-leading order (NNLO) (dashed-dotted 
green line) and the central values obtained by KLOE (dashed red line). The Monte Carlo results 
are normalized to the experimental data. 



Acknowledgements 

This publication has been supported by the European Commission under the 7th Framework 
Programme through the 'Research Infrastructures' action of the 'Capacities' Programme. Call: 
FP7-INFRASTRUCTURES-2008-1, Grant Agreement N. 227431. 



Bibliography 



[2] 
[3] 
[4] 
[5] 
[6] 
[7] 



D.G. Sutherland, Phys Lett 23, 384, (1966). 

C. Ditsche, B. Kubis, U-G. Meifiner, arXiv:0910.0210vl [hep-ph], (2009) 

W.A. Bardeen, et at, Phys Rev Lett 18, 1170, (1967). 

K.Nakamura et at, J Phys G 37, 075021, (2010). 

C. Roiesnel, T. Truong, Nucl Phys B 187, 293, (1981). 

J. Gasser, H. Lcutwylcr, Nucl Phys B 250, 539, (1985) . 

J.Bijnens, K.Ghorbani, JEEP 11, 030, (2007). 

A.Anisovich, H.Leutwyler, Phys Lett B 375, 335, (1996). 



22 



P. Adlarson 



PrimcNet Workshop 2011, Jiilich 



[9] J.Kambor, C.Wiesendanger, D.Wyler, Nucl Phys B 465, 215, (1996). 

[10] J.Bijnens, J.Gasser, Phys Scnpta T99, 34, (2002). 

[11] G.Colangelo, S.Lanz, E.Passemar, PoSCD09 047,(2009), arXiv:0910.0765vl [hep-ph]. 

[12] K.Kampf, M.Knecht, J.Novotny, M.Zdrahal, arXiv:1103.0982vl [hcp-ph] (2011). 

[13] S.P. Schneider, B.Kubis, C.Ditsche, JEEP 1102, 028, 201 (2011), arXiv:1010.3946v2. 

[14] R. Dashen, Phys Rev 183, 1245, (1969). 

[15] G.Colangelo, et al, arXiv:1011.4408v2 [hep-lat], (2010). 

[16] F. Ambrosini et al (KLOE Collab.), JEEP 05, 006, (2008). 

[17] H. H. Adam et al, arXiv:0411038 [nucl-ex] (2004). 



23 



BERNE-LUND- VALENCIA 
DISPERSIVE TREATMENT OF 77 

STEFAN LANtf 



3tt 



The decay rj — > 3-7T is of particular theoretical interest because it can only proceed through isospin 
breaking. If the strongly suppressed electromagnetic contributions are neglected, the decay am- 
plitude is proportional to (m u — md) or, alternatively, to the quark mass double ratio 



Q 



with 



mt 



1/ 

2^ 



m d ) 



(1.1) 



Furthermore, this process involves two theoretical puzzles. The first one is the fact that the 
predictions for the decay width from current algebra, but also from one-loop chiral perturbation 
theory (xPT), fail to reproduce the experimental value. It is understood that this is due to large 
final state rescattering effects that can be ideally treated using dispersion relations. Gasser and 
Lcutwyler, having included one final state rescattering process in the decay amplitude [1], found 
for the decay width of the charged channel T c = 160 ± 50 eV. This was in perfect agreement with 
the then-current PDG average T c — 197 ± 29 eV and the case seemed to be settled. However, 
the authors already pointed out a new experiment hinting at a larger decay width, which would 
"resurrect the problem" they "claim[ed] to solve". Figure 1.1 clearly shows that over the following 
years, the PDG average grew out of the error bars of the one-loop prediction. The second puzzle 
is related to the slope parameter a in the neutral channel, where xPT, but also an older dispersive 
analysis [2], fail to reproduce the experimental finding. 
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Figure 1.1: PDG average for the decay width of the two rj — > 3w channels since 1968. Gasser and 
Leutwyler used the value for the charged channel from 1982. 



In 1996 two dispersive analyses of rj — >• 3ir have been published [2, 3]. Our analysis follows 
closely the method proposed in the second reference, where a detailed description and further 
references can be found. Many developments in recent years, as well as current experimental and 
theoretical activity by other groups, have made a new analysis worthwhile. The inputs that the 
dispersive analysis relies on have been improved considerably: several groups have published accu- 
rate descriptions of the irn scattering phase shifts [4, 5, 6], and a modern precision measurement 
of rj — > 37r is available [7], which can be used in order to fix the subtraction constants. Even more 
experimental data is to be expected in the near future [8]. In addition, several theoretical works 
on this decay have appeared in the last few years (e.g., [9, 10, 11, 12]). In particular, there is also 
an analytical dispersive treatment [13] that was presented at this workshop. 

*in collaboration with G. Colangelo, H. Leutwyler and E. Passemar 
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The dispersive analysis relies on a decomposition of the amplitude into isospin amplitudes 
Mi(s) that are functions of one variable only, 

M[s,t,u) = M (s) + {s- u)Mi(t) + (s - t)Mi(«) + M 2 (i) + M 2 (u) - |m 2 (s) . (1.2) 
From unitarity and analyticity follow dispersion relations for the functions M/(s), 

oo 

i , , s 2 f ds' s'mS (s')Mo(s') , , < 

= <»» + 7 y ^ i» o(s0 | (s /-;-; £ ) ^ ■ 

and similarly for Mi(s) and M 2 (s). The functions M/(s) are angular averages over all of the Mj 
such that the dispersion relations are coupled. The Omnes function is given by 

"-w-^ii^"!' (L4) 

The dispersion relations contain a total number of four subtraction constants: ao, Aj and 7o 
in the equation for Mq and one more, /3i, in the equation for M%. These are left free by the 
dispersion relations and have to be determined otherwise. The dispersion relations are solved 
numerically by an iterative procedure, starting at some initial configuration for the Mi(s). If 
the subtraction constants are determined by a matching to the one-loop result from x?T (as in 
Ref. [3]), it turns out that there is a clear deviation from the experimental result in Ref. [7] for 
large s (see Fig. 1.2). Alternatively, we can determine the subtraction constants by a combined 
fit to the measured Dalitz plot and to one-loop xPT around the Adler zero of the amplitude, 
where the series is expected to converge best. Since the amplitude is proportional to Q~ 2 , its 
normalisation cannot be fixed from the experimental data. For this, Chiral Perturbation Theory 
remains the only source of information. The results from a preliminary fit are shown in Fig. 1.2. 
As expected, these curves agree well with the result from KLOE. 

Comparing the decay width that is calculated from the dispersive amplitude with the present 
PDG value of r c = 296 ± 16 eV, we can extract a value for Q. Using the subtraction constants 
from the pure matching to one-loop we get Q — 22.7 ± 0.7, thus updating the analysis from 
Ref. [3]. The preliminary result using the subtraction constants from the fit to the KLOE data is 
Q = 21.3 ± 0.6. In Fig. 1.3, these results are compared to a number of other values for Q from the 
literature. Note, in particular, the rather strong disagreement with the result from the analytic 
dispersive calculation by Kampf et al. The origin of the discrepancy is not yet fully understood, 
but it lies in the different procedure that is used in order to fix the normalisation, since Q is very 
sensitive to the latter. While we fit to the one-loop amplitude along the line s = u in the vicinity of 
the Adler zero, the authors of Ref. [13] fit along the line t — u. Doing the same with our dispersive 
representation results in a very strong violation of the Adler zero and is thus not a viable course. 




Figure 1.2: The squared amplitude for the decay 77 — > 7r°7r + 7r~ along the lines Y — (left) and 
X = 0.125 (right). X and Y are the usual Dalitz plot variables, the dashed lines represent the 
limits of the physical region. The uncertainty band is marked by the dotted lines. One can 
clearly see the disagreement between the dispersive result that relies only on one-loop xPT and 
experiment for Y < (which corresponds to large s). 
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Figure 1.3: Comparison of several results for the quark mass ratio Q. The two last dots represent 
our preliminary results. 

For the slope parameter in the neutral channel, we get a — —0.045 ± 0.01, in only mild 
disagreement with the PDG average. This value is in excellent agreement with the finding of Kampf 
et al., again hinting that the discrepancy in Q is not due to the shape but to the normalisation of 
the amplitude, since a is independent thereof. 

Note that all results presented here are preliminary. Last refinements on the numerical analysis 
are in progress and the final results will be made available soon. 
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1.1 Introduction 

As is obvious from the number of contributions dealing with 77 — > 37r decays presented in this 
workshop, these processes in recent years attracted and still attract both theoreticians and ex- 
perimentalists. The main reasons are the following. These processes proceed via isospin breaking 
effects. Moreover, since the electromagnetic contribution to them is strongly suppressed, to a 
good approximation they are proportional to (m„ — md) mass difference and thus by comparing 
the result of measurement of their decay rate with the theoretically predicted ones we have a 
direct method for determination of this quantity, in our work appearing as the quark mass ratio 
R = ™"Z™ 1 where m = \(m u + md). Secondly, recent measurements (cf. e.g. [1]) of the neutral 
decay -q — > 3tt° start to be precise enough to observe the effect of cusp occurring in this decay. 
Finally, there is still a theoretical challenge to explain satisfactorily the observed discrepancy be- 
tween the experimentally measured Dalitz plot parameters describing the energy dependence of 
the amplitudes and their values predicted by the two-loop computation in the framework of chiral 
perturbation theory (ChPT) [2]. 

This has inspired various studies exploring possible explanations of the discrepancy [3] , namely 
higher order final state rescatterings, influence of slow convergence of tttt scattering or of 77 — >• 3tt 
amplitude itself, unexpectedly large electromagnetic contributions, effects of resonances or with 
that connected possibility of incorrectly determined next-to-next-to-leading-order (NNLO) low- 
energy constants of ChPT CjS. On the experimental side there is substantial activity as well (for 
instance KLOE and WASA-at-COSY at this workshop [4] promised new data on charged decay, 
for which there exists currently only one precise analysis on at least 4 Dalitz plot parameters, [5]). 

The complication of those analyses which are trying to avoid the above mentioned possible 
problems of ChPT by using alternative approaches is that currently there exists no such alternative 
containing explicitly quark masses. This means that even if one found using them a correct 
description of the energy dependence of the r\ decays, in order to extract any information on the 
quark masses (i.e. to fix the normalization), one unavoidably needs to match it back to ChPT. In 
the case one used in his approach different assumptions than those of ChPT, this matching brings 
another sources of errors into the game. 



1.2 Analytic dispersive parametrization 

For the task of finding the correct value of R taking all this into account, we have employed in [6] 
our analytic dispersive parametrization valid to two-loop order, which is constructed using just 
basic assumptions of quantum field theories together with some observed hierarchy of various con- 
tributions to the amplitude and therefore can reproduce very well also the NNLO chiral amplitude. 
It is possible to include full isospin breaking effects connected with different masses of charged 
and neutral pions but at the current level of precision we work in the limit m 7r ± — m^o. In that 
limit the neutral decay is connected with the charged one and the parametrization of the charged 
decay is in this case in the form 

M x (s, t, u) = A X M 2 + B x (s-s) + C x {s ~ s) 2 + D x [(t - s) 2 + (« - s) 2 ] + E x (s - s) 3 

+ F x [(t - s) 3 + (u - s) 3 } + U(s, t,u), ' ' 

where the unitarity part U(s,t,u) contains the same parameters A X} B X , C Xl D x together with a 
few additional parameters describing tttt scattering (which are held fixed) and is given analytically. 

*IPNP, Faculty of Mathematics and Physics, Charles University in Prague, zdrahalSipnp.mff.cuni.cz 
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Using this parametrization and the data for charged decay from KLOE [5], we perform the 
following two distinct analyses. 



1.3 Analysis I: Correcting ChPT 77 — >■ 37r result 

By comparing the values of Dalitz plot parameters from the experimental determination with 
their predictions from NNLO ChPT, we have found that despite the discrepancy for the individual 
parameters, they fulfill two relations which are valid in the case the imaginary part of the chiral 
amplitude is determined correctly (the third relation which is simply the value of the neutral 
parameter f3 has not been experimentally measured yet). This means that the data indicate 
the possibility that all the difference between the physically measured amplitude and the one 
stemming from the ChPT computation can be described by a small real polynomial, i.e. it can 
also be included into the values of the parameters d. This inspired our first analysis, where we 
have added such polynomial to our parametrization of the ChPT amplitude and fitted it from the 
KLOE data. From that one could extract constraints on such phenomenological values of (7, and 
more importantly for our task we have obtained corrected value of R = 37.7 ± 2.9, where the error 
is estimated from the (slow) convergence of chiral expansion in the first leading orders. 



1.4 Analysis II: Direct fit to experimental 77 — > 37r data 

Our second analysis resigned on explaining the origin of the Dalitz plot discrepancy and instead as- 
sumed that independently on the way one obtained the correct amplitude from QFT computations, 
it would fulfill the general properties we have employed in the construction of our parametrization, 
i.e. it will be reproduced well by our parametrization. We have therefore fitted the experimental 
data on the parametrization. 

As was pointed out in the introduction, in order to obtain any information on the quark masses, 
we need to fix the normalization of the amplitude by matching to ChPT. Thanks to the fact we 
have analytic expressions, it is sufficient to find just one point where we perform the matching, 
i.e. the point where the chiral expansion of the amplitude converges fast. However, since there 
is no theoretical argument that would choose one such point for the matching, we have used the 
assumption that if such point exists, it will not be an isolated point but a region where the behavior 
of the NNLO chiral amplitude and the behavior of the amplitude fitted from data would be similar 
and also various chiral orders behave in that region well. We have found a prescription for such 
matching fulfilling these requirements and obtained as a result of this analysis R = 37.8 ± 3.3, 
where the error is fully dominated by the error of the fit to KLOE data. 



1.5 Conclusions 

We have used two analyses based on different assumptions, which can be however fulfilled simul- 
taneously by the genuine physical amplitude. Since both of the analyses lead to the compatible 
values for R, we expect that they provide us with some information about the amplitudes and 
also that the determined value of R is reasonable. 

In [7] we have shown how one can combine this determination with the isospin symmetric 
studies on lattice and those using sum-rules techniques in order to obtain individual quark masses. 
For example using FLAG [8] values, we obtain 

m u = 2.23(13) MeV, m d = 4.63(16) MeV, (1.2) 

where these values are taken in the MS renormalization scheme at /j, = 2 GeV. 

Let us conclude with the comment that it will be very important in the forthcoming experi- 
mental measurements of the charged r\ — > 37r decay to verify the validity of our assumptions since 
the current analysis is based on only one experimental information on 5 Dalitz plot parameters 
which is at our disposal [5]. 
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r] ->3tt IN RESUMMED %PT: 
DALITZ PLOT PARAMETERS 

MARIAN KOLESAR*\ JIRI NOVOTNY* 
AND SEE A S TIEN DESC TES- GENON^ 

We analyze the Dalitz plot parameters of r\ — > 37r decay in the framework of resummed chiral 
perturbation theory. This approach allows us to keep the uncertainties in the NNLO and higher 
orders under better control and estimate their influence. The compatibility of assumption of 
reasonably small higher order remainders with experimental data and NNLO standard xPT results 
is investigated. We calculate the effect of lowest order resonance exchange and irir rescattering 
bubble corrections up to three loops. 

The Dalitz plot parameters are usually defined as: 

rj -5- 7r°7r + 7r- : | A\ 2 = A%(1 + ay ± by 2 ± dx 2 ± . . . ) (1.1) 
r/^3^°: \A\ 2 = ^(1 + a* + . . . ) (1.2) 

where x ~ u — t, y <~ s — s, z ~ x 2 +y 2 and s is the Dalitz plot center s = 1/3(M 2 +3M 2 ). 
There is recent experimental information available [1, 2, 3, 4, 5, 6, 7, 8, 9], results from collab- 
orations producing values for both channels can be found in table 1.1, along with NNLO xPT 
results [10]. 

As can be seen, a discrepancy can be suspected for the charged decay parameter b and the 
neutral decay parameter a. Our aim is to try to understand whether the theory, by which we 
mean xPT as a low energy representation of QCD, really does have difficulties explaining the data 
and if so, try to identify the source of the problem. 

There is a long standing suspicion that chiral perturbation theory might posses a slow or 
irregular convergence in the case of the three quark flavor series [11]. An alternative approach, 
dubbed resummed xPT [12, 13], was developed in order to express standard assumptions in terms 
of parameters and uncertainty bands. The procedure can be very shortly summarized as: 

- standard %PT Lagrangian and power counting 

- possible irregularities in the expansion assumed 

- only expansions derived directly from the generating functional trusted 

- explicitly to NLO, higher orders collected in remainders 

- remainders not neglected, treated as sources of error 

- manipulations done in non-perturbative algebraic way 

Our calculation closely follows the procedure outlined in [15], more details can be found in [16]. 
A comprehensive work is in preparation [17]. 

In accord with the method, leading order low energy constants (LECs) are expressed in terms 
of convenient free parameters 

Z=%, X=^|§? , r = ^, R=p^±, (1.3) 
Fl F*M* m {m d -m u ) 
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a 


b 


d 


ct 


Crystal Barrel [1, 2] 
KLOE [3, 4] 


-1.22 ±0.07 
-1.090 ±0.020 


0.22 ±0.11 
0.124 ± 0.012 


0.06 ± 0.04 
0.057 ±0.017 


— 0.052 ± 0.020 
-0.0301 ± 0.0050 


NNLO xPT [10] 


-1.271 ±0.075 


0.394 ±0.102 


0.055 ± 0.057 


±0.013 ± 0.032 



Table 1.1: Some of the available experimental and xPT results. 
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where rh=(m u + m t ;)/2. To first order in isospin breaking the Dalitz plot parameters do not 
depend on R. At next-to-lcading order, the LECs L4-L8 are algebraically reparametrized using 
chiral expansions of two point Green functions. The 0(p 6 ) and higher order LECs, notorious for 
their abundance, are collected in a relatively smaller number of higher order remainders. 

We use several approaches to deal with the remainders. The first one is based on general 
arguments about the convergence of the chiral series [12], which leads to 

G = G^ + G^+A^, ~±0.1G, (1.4) 

where G stands for any of our 2- or 4-point Green functions, which generate the remainders. This 
is in principle an assumption. Hence we test the compatibility of this assumption of a reasonably 
good chiral convergence of trusted quantities with experimental data in a statistical sense. 

The framework of resummed %PT is well suited to include additional information about higher 
orders from various sources [15]. An independent estimate of the remainders can on one side be 
an important check of the validity of the statistical remainder estimate or could try to explain 
any deviations from this assumption. Here we employ a specific higher order calculation, namely 
3-loop bubble contributions to final state irir rescattering and an estimate of the influence of the 
lowest lying resonances using resonance chiral theory [14]. 

Using the approach outlined above, we have found the experimental values consistent with the 
assumption of small higher order remainders for the charged decay parameters. We cannot confirm 
the suspected discrepancy in the case of parameter b, as small uncertainties in higher orders could 
accommodate the difference. 

In the case of the neutral decay parameter a, we confirm the positive-sign prediction of xPT 
and find the experimental value incompatible with the assumption of good convergence properties 
in the center of the Dalitz plot for any scenario of spontaneous chiral symmetry breaking controlled 
by the free parameters X, Z and r. 

Contributions from resonance exchange have turned out to be small. On the other hand, 7T7t 
rescattering might explain the discrepancy in a, especially for low values of the pseudoscalar decay 
constant in the chiral limit. That could signal a failure of convergence of chiral series already at 
around 500MeV. 
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Figure 1.1: Dalitz plot parameters b and a at r=25: top Z=l, bottom Z=0.5 . Horizontal lines: 
dashed - KLOE measurement [3]; dotted - NNLO %PT [10]. Light band - statistical remainder 
estimate. Dark solid line - result including lowest lying resonance multiplets. 
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Figure 1.2: Dalitz plot parameters b and a at r=25: top Z=l, bottom Z=0.5 . Horizontal lines: 
dashed - KLOE measurement [3]; dotted - NNLO %PT [10]. Light band - statistical remainder 
estimate. Dark band - result including irir rescattering, depending on ren. scale /z=0.5-^lGeV 
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MODEL-INDEPENDENT APPROACH TO 

77 — >- 7r + 7r~7 AND 77' - 

FELIX STOLLENWERK^ 



7T + 7T 7 



1.1 Introduction 

We present a new, model-independent method [1] to analyze radiative decays of mesons to a 
vector, isovector pair of pions of invariant mass square below the first significant irir threshold 
in the vector channel. It is based on a combination of chiral perturbation theory and dispersion 
theory. This allows for a controlled inclusion of resonance physics without the necessity to involve 
models, like vector meson dominance, explicitly. 

1.2 The approach 

For the radiative decays at hand, there is a significant deviation between the predictions [2] of 
chiral perturbation theory (ChPT) and data, the source of which is mostly the non-perturbative 
7T7T final state interaction. Those are universally collected in the pion vector form factor Fy(s„), 
whose analytic properties can be exploited in the form of a dispersion relation to give 



Fv(s W7r ) = exp 



ds- 



Su(s) 



Ami 



2 (s 



ie) 



(1.1) 



where (r 2 ) and Sn(s) denote the mean square charge radius of the pion and the elastic 7T7t phase 
shift in the vector-isovector channel, respectively, and is the invariant pion mass. 



N — 10 - 




-0,2 0,2 0,4 0,6 0,8 



s [GeV ] 



Figure 1.1: The (red) solid band shows the form factor derived from Eq. (1.1), the (blue) dashed line 
the result from one-loop ChPT — both with identical values for the pion radius. The time-like data 
are from Refs. [3, 4], whereas the space- like data are from Ref. [5]. The short (long) thick, horizontal 
bar denotes the kinematic range covered in the decay of the 77 (rf). In the 77 range, the form factor 
can be approximated by \F v (s^)\ = 1 + (2.12 ± O.Ol^™ + (2.13 ± 0.01)4^ + (13.80 ± 0.14)4^. 

Eq. (1.1) fixes the form factor only up to a multiplicative function without right-hand cuts. As 
the transition amplitude features the same analytic structure as Fv{s 7r7! ) apart from a kinemati- 
cally and chirally suppressed left-hand cut, we may factorize the amplitude according to 

\AP(s^)F v {s^)\ 2 , (1.2) 
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and expand P(s- K17 ) in a Taylor series around s Tn = 0: 

P( 8int ) = l + a^a mt + 0(4J) . (1.3) 

Assuming that the universal part Fy(s niT ) comprises all non-perturbative effects, the reaction 
specific part j4P(s„) should be accessible to a perturbative treatment. Indeed, we find that a 
linear polynomial (1.3) with 

a = (1.96 ±0.27 ±0.02) GeV~ 2 ; a' = (1.80 ± 0.49 ± 0.04) GeV~ 2 (1.4) 

is already sufficient to describe spectral decay data, see Fig. 1.2. In principle, p — to mixing distorts 
the rf spectrum. If one assumes the mixing strength to be identical to that in the vector form 
factor (based on the quark model the effect is much smaller), one finds the wiggly (green) line. 
Clearly, the data is not yet sufficient to extract information on the magnitude of the effect. 




Figure 1.2: Experimental data and error weighted fits for 77 (left, data are from Ref. [6] (filled 
squares) and Ref. [7] (open circles)) and rf (right, data are from Ref. [8]) to 7r + 7r~7 according to 
Eqs. (1.2) and (1.3) with s vv — m„(/)(m„w — 2E 1 ). The wiggly (green) line in the right panel 
denotes the possible impact of p — u> mixing under the assumption that it appears here with the 
same strength as in Fy. 

The postulated decay amplitude (1.2) can be matched to one-loop 17(3) extended ChPT. 
Concerning the s^-dependent part, this connects the slope parameters a^'' to the low energy 
constant 02 giving 

a\= 9.70 ±0.7; a( = 9.23 ±1.4. 

Thus, the data sets are consistent with the assumption that the only non-perturbative part of the 
amplitude originates from the pion vector form factor, and Eq. (1.2) provides a controlled way to 
disentangle perturbative and non-perturbative effects in order to methodically extend perturbative 
calculations to the resonance region. By invoking chiral Ward identity and large N c arguments, 
essentially half of the parameter cr'' can be physically interpreted. Further possible ways of 
interpretation are subject of current research. 
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STUDY OF rj — ► tt+tt^ DECAY 



CAMILLA DI DONATO* 
for the KLOE and KLOE-2 Collaborations 

1.1 Introduction 

We measured the ratio R v = T(r] — > ir + ir~j)/T(ri — > 7r + 7r~7r ), analyzing an integrated luminosity 
of 558 pb^ 1 , collected by KLOE at the DA$NE e + e~ collider. The 77 — > tt + it~j process is sup- 
posed to proceed both via a resonant contribution, mediated by the p meson, and a non resonant 
direct term, connected to the box anomaly. The presence of the direct term affects the partial 
width value. 

The Chiral Perturbation Theory (ChPT) provides accurate description of interactions and de- 
cays of light mesons [1]. The decays rj — > ir + ir~j and 77' — > 7r + 7r~7 are expected to get con- 
tribution from the anomaly accounted for by the Wess Zumino Witten (WZW) term into the 
ChPT Lagrangian [2]. According to effective theory [2] the contribution of the direct term 
should be present together with VMD. In case of 77 — > 7r + 7r~7 the p contribution is not dom- 
inant, this makes the partial width sensitive to the presence of the direct term, while in case 
of 7/ — > 7r + 7r~7 the partial width is dominated by the resonance and the direct term effect 
should be visible in the dipion invariant mass distribution. The present world average of the 
77 — > 7r + 7r~7 partial width, T(rj — > Tr + TT~-f) = (60 ± 4) eV [3], provides strong evidence in favour 
of the box anomaly, compared with value expected with and without the direct term, respec- 
tively (56.3 ± 1.7) eV and (100.9 ± 2.8) eV [2]. Recently CLEO [4] has measured the ratio 
R v = 0.175 ± 0.007 S i at ± 0.006 sys t, which differs by more than 3cr from the average result of 
previous measurement [5, 6], R v = 0.207 ± 0.004 [7]. Here we present a new measurement with 
the highest statistics and the smallest systematic error ever achieved. 

1.2 Event selection 

The analysis has been performed using 55 8 pb" 1 of the 2004-2005 data set acquired by KLOE at 
the DA$NE e+e~ collider (y/s ~ 1.02 GeV). The final state under study is 7T + 7r 77, since the 
77 mesons are produced together with a monochromatic recoil photon {E-, — 363 MeV) through 
the radiative decay <j> — > 777. In the considered data sample about ~ 25 x 10 6 77's are produced. 
The main background comes from <j> — >■ n + n~ it , tt° — > 77 decaying to the same final state. Other 
backgrounds are 4> — > 777 — > 7r + 7r~7r° — > ir + iT~i^ with one photon lost, and <f> — > 777, r\ — > e+e _ 7 
when both electrons are mis-identified as pions. 

As first step of the analysis, a preselection is performed, requiring at least two tracks with oppo- 
site charge pointing to the interaction point (IP) and at least two neutral clusters in time (not 
associated to any track), having energy E c i > 10 MeV and polar angle in the range (23° — 157°). 
Tracks are sorted according to the distance of the point of closest approach from the IP. The first 
two tracks are selected. We require the most energetic cluster (7^) to have E c i > 250 MeV and we 
identify it as the photon recoiling against the 77 in the (j> — > 777 decay. Moreover we ask for 7^, inside 
the calorimeter barrel (55° — 125°), to avoid effects of cluster merging between barrel and end-caps 
of the calorimeter. Other cuts are imposed to clean up the sample; cut on cluster-track collinearity 
and identification by time of flight (TOF) are used to reject electrons. The cut effectively rejects 
Bhabha background and other processes with electrons in the final state. To select 77 decays we 
exploit the <f> VI two body decay kinematic computing the 74, energy, using only the 70 polar 
angle: 

- - - J. - TP m l~ m v 

Pt-Pn+lh *-r*- 2 (E4-fo\co8<p) 

where ip is the angle between the average </>-meson momentum measured run by run with high 
accuracy and 7^. This allows us to improve the energy measurement of the recoil photon to 0.1%. 

*INFN Naples, camilla.didonatoana.infn.it 
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We can determine the direction of the photon from rj decay using <j) and 7r-mesons information: 

P~f v =P4>~ P*+ ~ Pit- ~ P14, 

The photon direction is compared with the direction of each neutral cluster, Aip = ip ciu — . 
If no cluster within Aip < 8.5° is found the event is rejected. The cluster with the minimum 
value of Aip is selected for further analysis. In order to reject the cf> — > ir + ir~ir° background, the 
angle between the two photons in the ir° reference frame, evaluated using the 4> and the 7r-mesons 
momenta, is calculated and rejected with an ang ular cut ip* y * 7 < 165°; in order to reduce the 
systematics the angle is evaluated in the transverse plane. Finally we select events requiring 539.5 
MeV < M„+„-^ < 554.5 MeV (Fig. 1.1). 



30000 



Experimental data 
Sum of all background 
Other r\ decays MC 
Other background 
Signal: T|->7r + Jt7 




600 650 

M(7t + 7T y) 



700 

[MeV] 



Figure 1.1: The n + ir 7,, invariant mass distribution: Data-MC comparison. Dots are data, 
Magenta is MC signal 7/ — > 7r + 7r _ 7, Red is all MC background contribution 



1.2.1 Normalization Sample: rj — > 77+77 77 

The process 4> —> 777 with 77 — > 7r + 7r~7r° represents a good control sample, due to the similar 
topology. Moreover the ratio T(rj — > ir + ir~j)/r(rj — » 7r + 7r~7r°) is not affected by the uncertainties 
on the luminosity, the <j) — > rj"f partial width and the <p production cross section cancel in the 
ratio. We use the same preselection as for the 77 — > 7r + 7r~7 signal and calculate the missing 
four-momentum: 

where the variables in the formula represent the four-momenta of the <f> meson and the products 
of the decays. For the 77 — > 7r + 7r~7r° signal, the missing mass peaks at the ir° mass value and we 
select events with \M m i SS — m v a \ < 15 MeV. The remaining background is rejected very efficiently 
by using an angular cut applied to the two photons from the tt° decay, cp^l > 165°; the angle is 
evaluated in the transverse plane. We select ^(77 — > 7r + 7r~7r°) = 1115805 ± 1056, with a selection 
efficiency of e = 0.2276 ± 0.0002 and a background contamination of 0.65%. 



1.3 Results 

The total selection efficiency of the 77 — > 7r + 7r~7 signal is e = 0.2131 ± 0.0004. Background 
contribution and the signal amount in the final sample are evaluated with a fit to the E m i SS — P m iss 
distribution of the 7r + 7r~70 system with the shapes from remaining background and signal MC in 
the range \E miss - P miss \ < 10 MeV. We find ^(77 tt+tt^) = 204950 ±450 with a background 
contamination of 10%. Combining our results we obtain the ratio: 

FT?? — y 7T^7T~ 

Rv = ^T 1 1 ttt = 0.1838 ± 0.0005 stQt ± 0.0030 sast (1.1) 

I (77 — >• 7T+7T 7T U ) 
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Source of uncertainty Relative error 



V^ 71 1 < 165° ± 2° ±0.6% 

A V > 8.5° ± 2° ±0.4% 

\M 7r + r -_ / - M„| < 7.5 MeV ± 2 McV ±0.6% 

\M m i BB - M^ol < 15 MeV ± 4 MeV ±0.4% 

^>165°±2° ±0.1% 

B 1 <n > 10 MeV ± 2 MeV ±0.1% 

> 250 McV ± 4 McV ±0.1% 

" Fit E miss - P mz33 ±0.6% 

Total 1.6% 



Table 1.1: Summary table of systematic uncertainties. 



to be compared with world average value T(r/ — > tt + tt^j)/T(i] — > 7r + 7r~7r ) = 0.202 ± 0.007 [3]. 
Our result is in agreement with a recent CLEO [4] measurement, which differs by more 3 a from 
the average of previous results [5, 6, 7]. 

The systematic uncertainties due to analysis cuts have been evaluated by varying the cuts on 
all variables and re-evaluating the branching ratios. The relative variation for each source of 
systematic is in Tab. 1.1. The total error is taken as the quadratic sum of all contributions. 

The M^+ v - dependence of decay width has been parameterized in different approaches, in 
which VMD has been implemented in effective Lagrangians [2, 8]. Preliminary comparison between 
dipion invariant mass, with the most simple approach as from [8], shows a good agreement with 
data. Fit with model independent parameterizations, as from [9], is in progress. 

1.4 Conclusions 

Using a data sample corresponding to an integrated luminosity of 558 pb _1 , we select 204950 
rj — > 7r + 7r _ 7 events and 1115805 rj — > 7r + 7r~7r°, from the <f> — > tjj decays. The corresponding width 
ratio is: R n — 0.1838 ± 0.0005 s t a i ± 0.0030 S j, S i Our measurement is in agreement with the most 
recent result from CLEO [4], which is R v = 0.175 ± 0.007 stat ± 0.006 syst . 

Combining our measurement with the world average value T(r] — > n+ir^n ) = (295 ± 16) eV [3], 
we get r(?7 — > 7r + 7r~7) = (54.2 ± 0.3) eV, which is in agreement with the value expected taking 
into account the direct term [2], providing a strong evidence in favour of the box anomaly. 
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1.1 Introduction &; motivation 

The decay channel i] — > it + tt~j provides the opportunity to study QCD anomalies at the chiral 
limit. The decay width and the shape of the Ely-distribution of this channel are sensitive to the 




Figure 1.1: Left: Wess-Zumino-Witten-Lagrangian (Lwzw) including the triangle anomaly (first) 
and box anomaly (second) term [1]. Right: Measured E 1 distribution for the reaction: pd — » 
3 He [77 — > 7r + 7r _ 7] [4]. The measurement has been performed with high statistics and was fully 
efficiency corrected. The measured E 7 distribution is described by multiplying the simplest gauge 
invariant matrix element (blue curve) with a form factor (dashed and red curve) which depends 
on a single parameter a [5]. 

box anomly term which is part of the Wess-Zumino-Witten-Lagrangian (see Fig. 1.1). However, 
the theoretically predicted decay width and _E 7 -distribution do not agree with the experimental 
results, if hnal state interactions are not included by unitarized extensions of the Lwzw- The 
experimental observables for testing these extensions are (i) the branching ratio [2, 3] or (ii) the 
distribution of the single photon energy [3]. A recent measurement of the photon energy (see 
Fig. 1.1) with the WAS A detector can be found in [4]. 

The aim of this work is to measure the branching ratio and the single photon energy distribution 
in one experiment using the reaction: pp — > pp[rj — > 7r + 7r~7]. The data have been acquired during 
an 8 week experiment in spring 2010. At least 10 6 rj —> 7r + 7r~7 events are expected. 

The WASA detector (see poster by F. Goldenbaum) is devided into two parts: the forward 
detector and the central detector. 

Scattered projectiles or charged recoil particles arc identified in the forward detector, a tracking 
detector with a set of range hodoscopes, which uses the ^-method for particle identification and 
reconstruction. 

Decay products are measured with the central detector. Charged particles are identified by 
inspecting the deposited energy as function of the momentum as measured with a tracking device 
in a magnetic field. Neutral particles are detected in the calorimeter. 

*IKP1 Forschungszentrum Juelich, Germany, d. lerschafz-juelich.de 
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1.2 Data analysis 

As a first step, the reaction r\ — > Tr + ir~ir° has been analysed. On the one hand this reaction is im- 
portant for the determination of ' 011 ^ e °t ner hand it contributes to the background. 
The 7T° is also used for monitoring the calorimeter calibration. 5% of the available preselected 
data have been analysed with the recent analysis chain. 

The missing mass spectrum deduced from two protons is shown in Fig. 1.2. The 77-peak is 
visible, along with background, from direct pion production. The following cuts on the invariant 
and missing mass spectra of the decay products were employed to extract the ^-signal from the 
direct pion background: (i) The invariant mass of two photons had to be within a 0.11 GeV/c 2 
window around the 7r°-mass (ii) The missing mass deduced from two photons had to be larger 
than the rest mass of the two chraged pions (iii) The missing mass squared deduced from the two 
charged pions had to be positive. 

The energy and momentum balance (Fig. 1.2 right) is used to check the effect of each cut and 
investigate the background contributions. 




Figure 1.2: Left: Missing mass deduced from the initial state and the two final state protons. The 
red line indicates the 77-mass at 0.548 GeV/c 2 . This spectrum is obtained after proton reconstruc- 
tion and applying the cuts mentioned above. Right: Plot of the missing momentum as function 
of the missing energy (derived by the momentum/energy of all initial and final particles in the 
reaction). 



1.3 Summary &; outlook 

The rj-peak is clearly visible after applying cuts on invariant and missing mass related to the decay 
products. Background channels shall be identified in a second step using Monte Carlo simulations 
and analysing more data runs. 

The goal is to investigate the reaction i] —} 7r + 7r~7 and to determine the branching ratio. 
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RADIATIVE DECAYS 
OF PSEUDOSCALARS AND VECTORS 

CARL A TERSCHL USEN*^ and STEFAN LEUPOLD* 



1.1 Introduction 

Due to the running coupling constant in QCD, one cannot use perturbation theory for the low- 
energy regime. Instead one can use effective theories with hadrons instead of quarks as relevant 
degrees of freedom. We use a new counting scheme [1] which treats both vector mesons (V") and 
pseudoscalar mesons (P) as soft, i.e. of the order of a typical momentum q: 

m v , m P , ~ q. (1.1) 

In addition, we use large- N c arguments (with N c denoting the number of colors) to suppress ad- 
ditional flavor traces. We will use this counting scheme to describe radiative decays into dileptons 
A more detailed introduction is given in these proceedings in [2]. 



1.2 Decays V -> Pl + l~ and P -> Vl + l~ 

Using the counting scheme (1.1), the leading-order Lagrangian for the decays V — > Pl + l~ and 
P — > Vl + l~ can be determined as 

£vcc. = -^jh A e^ a Hv{[V^,d T V Ta } + - ^Lb A e^ a ?tr{[V^,V a0 } + [$, Xo ] + } 

- ^ n A e^ 3 tr {V^V a[j } m - tr {V^Q} d»A v . (1.2) 

Thereby xo = diag {m 2 ., m 2 , m 2 K } , Q = diag {2/3, — 1/3, — 1/3} and the matrices V^„ and $ 
describe the vector mesons (represented by antisymmetric tensor fields) and the pseudoscalar 
mesons, respectively, i.e. 

p% + uj^ V2p+ V s/2K+ v 
V := ( V2p- V -p%+co^ V2K° U ), (1.3) 
V2K~„ V2K°„ V2^ 

r° + ^r/ 8 V2tt+ V2K+ 
<I> := I V27T- ~^ + ^s V2K° (1-4) 



V2K- V2K° 



75^ 



All three terms proportional to h A , b A and n A are of order q 2 : The h A term because of the two 
derivatives, the b A term because of the quark-mass insertion proportional to xo an d the n A term 
because of the separate traces (note that 771 is proportional to tr{$}). Obviously, this leading- 
order Lagrangian (1.2) allows only for an indirect decay via a virtual vector meson. Furthermore, 
for the description of r\ and 77' mesons one has to take 77-77' mixing into account, 

77 = cos 9 ?78 — sin 9 r\\ , (1-5) 
77' = sin 9 ?7 8 + cos 9 771 with 9 w -20°. (1.6) 

The decay of the virtual photon into a dilepton is described by usual QED. 
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1.2.1 Parameter determination 

To determine the parameters Jia, bA and tia, the partial-decay widths for the two-body decays 

uj — > 7r°7, uj — > 777, <fi — > 777, — > 7/7, 77' — > CJ7 (1-7) 

are compared to experimental data. First, the value Ha is mainly determined via the decay 
uj — > 7r°7 and, therefore, is fixed as 

ft A = 2.32 (1.8) 

as in previous calculations [3]. The parameters bA and ua are then determined such that all 
considered two-body widths are approximately equally well described. It turns out that one does 
not find fixed values for this parameters but intervals 

b A € [0,0.19], n A S [0,0.05]. (1.9) 

For decays into dileptons no additional parameters are needed and we have predictive power. 

1.2.2 Results 

In Fig. 1.1, the form factors for the transition uj — > 7r° and 77' — > cj are plotted (solid lines). Both 
are compared to the standard vector meson dominance (VMD) form factors (dashed lines). For 
both transitions, there is a deviation from VMD visible. Whereas the uj — > n° form factor is mainly 
determined by the term proportional to Ha, there is a sizable influence of bA and tia on the rf — > uj 
form factor yielding a larger uncertainty compared to uj — > tt° . Obviously, data taken by the NA60 
collaboration for the uj — > it transition form factor [4] is much better described with our approach 
than with standard VMD. Additionally, the calculated results for the integrated rates are within 
acceptable intervals and agree well with the available experimental data [5] if existent: 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.04 0.08 0.12 0.16 

ny r [GeV] m n - [GeV] 

Figure 1.1: Left-hand side: oj — > 7r° form factor compared to dimuon data taken by the NA60 
collaboration [4]. Right-hand side: 77' — > uj transition form factor. 



1.3 Decay P 

For the transition of a pseudoscalar meson into a (real) photon, we use a more phcnomcnological 
approach by including the Wess-Zumino-Witten action [6] . On the one hand, a pseudoscalar meson 
can decay into two (real or virtual) photons via two virtual vector mesons. In the counting scheme 



41 



PrimcNct Workshop 2011, Jiilich 



C. Terschliisen 



(1.1), this is in leading order described by the Lagrangian C vcc . (1.2). On the other hand, the 
leading-order ChPT Lagrangian, the effective Wess-Zumino-Witten Lagrangian [6] 

Avzw = Jjy tr {QH} d,A a d„A fj + 0(4> 2 ), (1.11) 

allows for direct decays. Formally, this term is of order q 4 in our counting scheme (1.1) but nev- 
ertheless important as the low-energy leading-order term. Furthermore, the relative sign between 
both Lagrangians was fixed as negative by comparing the calculated rj — > 7 transition form factor 
to experimental data. 

The rj — > 7 transition form factor is plotted on the left-hand side in Fig. 1.2 (solid line) in 
comparison to the standard VMD calculation (dashed line) and data taken by the NA60 collabo- 
ration for the decay into a dimuon [4] . There is nearly no deviation between the two form factors 
visible, both describe the data very well. On the right-hand side in Fig. 1.2, the rj — > 7 transition 
form factor is plotted (solid line). Again, it is compared to the standard VMD calculation (dashed 
line). Here, a small deviation between our calculation and standard VMD is visible. Note that 
both p° and u are visible in this form factor. 
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Figure 1.2: Left-hand side: rj — > 7 form factor in comparison to dimuon data taken by the NA60 
collaboration [4]. Right-hand side: rj' — > 7 form factor. 



1.4 Summary and outlook 

The counting scheme (1.1) was used to calculate transition form factors in leading order. Thereby, 
the oj — » 7T° form- factor data was much better described than with standard VMD, the 77 — » 7 
form- factor data as good as with standard VMD. The presented partial decay widths are within 
acceptable intervals and agree well with the available experimental data. As a next step, next-to- 
leading order calculations have to be performed. 
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WASA-at-COSY 

CARL-OSCAR GULLSTROM* 
for the WASA-at-COSY Collaboration 

1.1 Introduction 

7T° is the lightest known hadron and can only decay via the electromagnetic interaction. The most 
common decay is into 27 (BR 98.82 %) but the interesting physics occurs when decays include 
leptons. The aim for WASA-at-COSY is to measure the decays 7r° — > e + e~7, tt° —> e + e~e + e~ and 
7T° -> e + e~. In tt° -> e+e~7 (BR 1.2%) and vr° ->■ e+e-e+e" (BR 3 • 1CT 5 ) the transition form 
factor reveals the electromagnetic structure of the decaying meson. The e + e~ mass distribution 
probe the existence of a new hypothetical boson that couples to the virtual photon [1]. The 
current upper limit for ir° — > e + e~j is set by the SINDRUM collaboration [2] based on 100.000 
events above 25 MeV e + e~ invariant mass. For ir° — > e + e~ the KTeV Collaboration found 795 
events in Kl — > 37r° [3]. They measured a BR of (7.49 ± 0.29) x 10~ 8 which exceeds the standard 
model prediction (6.23 ± 0.09 x 10~ 8 ) [4]. This leads to some speculation about new physics. One 
promising teory is a new boson U that couples both to quarks and leptons [5]. Excess positrons 
from the U boson could also explain the intensity and shape of the 511 keV line from the galactic 
center [6]. 

1.2 Analysis 

During spring 2010 a one week test measurement has been performed with WASA-at-COSY in 
Julich. In the experiment proton-proton collisions were used as ir° sources. The kinetic beam 
energy was set to 550 MeV in order to maximize the ir° cross section below threshhold for two 
pion production. The advantage of excluding tt~ production is a clean e~ sample for negatively 
charged tracks. 




'0 0.05 0.1 0.15 0.2 0.25 0.3 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.1 8 0.2 



/Mass e Vy (GeV/c*) Mass e V (GeV/c 1 ) 

Figure 1.1: Left: Invariant mass distribution of e + e~j. Right: Invariant mass of e + e~ pair from 
the e + e~7 peak. Black: Data; Green: MC ir° — > e + e~7; Blue: Background from pair production 
in the detector 7r° —> (7 —> e + e~) 7; Red: Total MC sample ir° —> e + e~j. 

1.2.1 7T° e+e"7 

In 7T° — > e + e~7 1.5 million reconstructed e + e~7 events are available. Out of those 50% is back- 
ground from external conversion of one photon in the ir° — > 27 decay. The main cut for reducing 
such events is the vertex cuts since the main part of the photons convert at the beamtube 30 mm 
from the interaction point. After cuts 40.000 ir° —> e + e~7 events are left with low background and 
high invariant mass resolution ( ~ 4%). The Dalitz distribution in Fig. 1.1 shows no U boson peak 
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and follows the existing 7r° — > e + e~7 transition formfactor. Tighter vertex cuts at the expense of 
efficiency are possible to cxlude the pair production peak at 20 MeV, so that an upper limit of 
U-boson exclusion could be set in this region. 

1.2.2 7T° e+e" 

In order to identify the ir° — > e + e~ candidates in the high-energy tail of the 7r° — > e + e~7 events 
one needs good momentum resolution. The vertex points serve as a check how well momentum was 
reconstructed by the central tracker. If the curvature is bad also the point of closest approach the 
the center is bad. To reach the desired resolution even tighter cuts are set than in the ir° — > e + e~7 
case. Fig. 1.2 shows a peak of 20 tt° — > e + e~ event candidates. MC simulations show that there 
should be about 25 % background from ir° — > e + e~7 with a lost photon. The external conversion 
is completely removed in this region by the vertex cut. 
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Figure 1.2: Left: Missing mass of 2 protons in the Forward Detector vs. invariant mass of e + e 
pair. Right: invariant mass of e + e~ pair in the central detector. Blue filled: Data; Green: 



MC e + e 7 background; Red: MC simulation of 7r° 
background. 



e + e plus e + e 7 and pair production 



1.3 Conclusion and Outlook 

For the one week test run we reconstructed 40.000 it —> e + e~-f. This data could be used to 
extend the upper limit on U Boson searches to lower masses. No deviation from the ir° transition 
form factor has been found. The 15 event candidates of ir° — > e + e~ after background subtraction 
motivate a longer production run in the future. The test run has not been optimized with respect 
to trigger and calibration and hence it is possible to reach 50 7r° — > e + e~ events per week. 

Bibliography 

[1] M. Reece and L.-T. Wang JHEP 07.051 (2009). 

[2] R. Meijer Drees et al. (SINDRUM 1 Collaboration), Phys.Rev.Lett. 68(1992) 3845. 

[3] E. Abouzaid et. al, Phys. Rev. D 75 (2007) 012004. 

[4] A. Dorokhov, M. A. Ivanov, Phys. Rev. D 75 (2007) 114007. 

[5] C. Boehm and P.Fayet, Nucl. Phys. B 683 (2004) 219. 

[6] Y Kami, M. Schmitt, T.M.P. Tait, Phys. Rev. D 78 (2008) 115002. 



44 



STUDY OF THE 77 -> e+e ^ DECAY 

MALGORZATA HODANA*™ and PAWEL MOSKAL*^ 
for the WASA-at-COSY Collaboration 

1.1 Introduction 

Since the r\ meson is a short-lived, neutral particle, it is not possible to investigate its structure 
via the classical method of particle scattering. To learn about its quark wave function, one 
studies the decay processes of this meson, in which a pair of photons is produced, at least one of 
them being virtual. The virtual photons have a non-zero mass and convert into lepton-antilcpton 
pairs. The squared four-momentum transferred by the virtual photon corresponds to the squared 
invariant mass of the created lepton-antilepton pair. Therefore, information about the quarks' 
spatial distribution inside the meson can be achieved from the lepton-antilepton invariant mass 
distributions by comparison of empirical results with predictions, based on the assumption that 
the meson is a point-like particle. The latter can be obtained from the theory of Quantum 
Electrodynamics. The deviation from the expected behaviour in the leptonic mass spectrum 
expose the inner structure of the meson. This deviation is characterized by a form factor. It is 
currently not possible to precisely predict the dependence of the form factor on the four-momentum 
transferred by the virtual photon in the framework of Quantum Chromodynamics. Therefore, to 
perform calculations, assumptions about the dynamics of the investigated decay are needed. 

The knowledge of the form factors is also important in studies of the muon anomalous magnetic 
moment, = (g^ — 2)/2, which is the most precise test of the Standard Model and, as well, may 
be an excellent probe of new physics. The theoretical error of calculation of a M is dominated by 
hadronic corrections and therefore limited by the accuracy of their determination. At present, 
the discrepancy between the a M prediction based on the Standard Model [1] and its experimental 
value is (28.7 ±8.0) • 1CT 10 (3.6er) [2]. 

The aim of the data analysis is the investigation of the electromagnetic structure of the rj meson 
by determining the transition form factor using the i] — > e + e~j decay mode. The probability of 
creation of a diclectron pair in considered decay is proportional to the probability of emission of 
a virtual photon with a time-like four-momentum. The square of this four-momentum vector is 
equal to the square of the mass of created e + e~ pair. By studying the probability of given decay as 
a function of the dilepton pair mass, one obtains information about the hadron-photon transition 
and hence about the electromagnetic structure of decaying neutral meson [3]. 

1.2 Current status of the analysis 

Data were collected using the pd — > z Her\ reaction at proton beam momentum of 1.69 GeV/c. 
The experiment was performed using the WASA-at-COSY detector [4] in November 2008. Data 
collected during 4 weeks, yielded approximately 10 million -q mesons tagged by the 3 He ions 
measured in the Forward Detector. 

The event selection in the Central Detector aims at choosing the decay channel of interest. We 
demand (i) that two tracks corresponding to oppositely charged particles are reconstructed and, 
(ii) that at least one neutral particle was registered and, (iii) that the signals are correlated in 
time with signals observed in the Forward Detector within a 12 ns and a 31 ns window for charged 
and neutral tracks, respectively. The energy deposited by a neutral particle {E~ t ) is demanded to 
drop linearly with increasing opening angle with the nearest charged particle (f2) starting from 
100 McV for Q = 0° and having E y = MeV for Q = 180°. Additionally, in the r\ centre of mass 
system, there has to be only one such neutral candidate, forming with a lepton pair a A</> angle 
in the range from 60° to 300°. The identification of electrons is achieved using energy- momentum 
plot as shown in Fig. 1.1 (right). An additional restriction is imposed on the missing mass of the 
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pd — > Xe + e~"f in the range of 2.66 — 2.87 GeV/c 2 , futher reducing background coming from rj 
decay channels with pions. The efficiency of supression of the rj — > 7r + 7r~7 and rj — > 7r + 7r~7r° can 
be seen in Fig. 1.2. 

Despite of asking in the analysis for events with two oppositely charged particles, there is a 
large background coming from the rj — ¥ decay channel (the one with the highest branching 
ratio). Since photons may undergo conversion at the beam pipe, the rj meson decaying to two 
photons, may give the same signature at the end as the decay of interest. One can see that in 
Fig. 1.2 (left) which shows the M e + e ~ spectra from experiment and simulations before applying 
cut on photon's conversion. In order to suppress background from rj — > 77 decay, where one 
of photons undergoes conversion into e + e~ pair in the beam pipe, the radius of point of closest 
approach (Rca) is used as presented in Fig. 1.1 (left). 




Figure 1.1: Left: experimental distribution of the radius of the point of the closest approach 
between two helices as a function of the e + e~ mass (taken from the leptons' three momenta 
calculated at the beam pipe). The x-axis is zoomed in to see events from the pd — > 77 reaction. 
The region outside the black lines is discarded from further analysis. Right: experimental spectrum 
of the energy deposited in the Electromagnetic Calorimeter as a function of particles' momenta. 
Electrons arrange themselves on the band, along the momentum- £^ ep diagonal (between black 
and red lines). 




Figure 1.2: The M e + e - spectra from the simulations, and the experiment (after prompt pion 
background subtraction), before applying the cut on the photon conversion (left) and with this 
cut (right). The spectra are not acceptance corrected. 

The M e + e - spectrum after selection of rj — > e + e~7 decay channel as described above is shown in 
Fig. 1.2 (right). Approximately 700 events of r/ — > e + e~j has been reconstructed. The extraction 
of the transition form factor and the systematic studies are in progress. 
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ANALYSIS OF THE DOUBLE DALITZ DECAY 

rj — > e + e~e + e~ 

PATRICK WURM* 
for the WASA-at-COSY Collaboration 



1.1 Introduction 

The Wide Angle Shower Apparatus (WASA) is a large-acceptance detector to study the decay 
channels of light mesons [1, 2]. It is operated at the Cooler Synchrotron COSY Jiilich. 

A large number of rj mesons is being produced in proton-deuteron and proton-proton collisions. 
The data permit the study of very rare 7/-decay channels, like the double Dalitz decay, where the 
r) meson decays via two virtual photons into two electron-positron pairs. 

So far, the only published value for the branching ratio is BIZ = (2.4±0.2 sta t±0.1 sys t) x 10~ 5 [3]. 
One objective of the WASA-at-COSY experiment is to determine a value for this branching ratio. 
A first analysis is based on apd-> "Tie 77 beam time from 2008, which comprises approximately 
10 7 77 mesons [4]. The analysis presented here is based on a total amount of 3 x 10 7 i] mesons, 
which includes additional 2 x 10 7 i] mesons, which have been produced using the same reaction in 
2009. 



1.2 Analysis chain 

In the production reaction pd — »■ Tic 77 the 77 meson is tagged by the Tie nuclei, measured in 
forward tracking detectors and an arrangement of hodoscopes. Tie can be identified with the 
AE — E method, where the energy deposit in a thin plastic scintillator is plotted versus the energy 
deposit in the stopping layer of the hodoscope. The Tie-missing mass distribution is used as a 
monitoring spectrum to count the number of 77 events after every cut. 

The central part of the detector is used to reconstruct and identify the 77-decay products. First, 
it is checked if the event has two or more positively and negatively charged particles. The charge 
of a particle is determined by the curvature of the track helix. The helix is measured by a Mini 
Drift Chamber, embedded in a solenoid magnet. 

Due to the large amount of pions, an effective method to distinguish electrons (positrons) and 
pions must be found. WASA provides particle identification with the momentum information, 
which is measured by the Mini Drift Chamber, combined with the energy deposit in a Plastic 
Scintillator Barrel and in an Electromagnetic Calorimeter. Fig. 1.1 shows the electron (positron) 
and pion bands. 

The decision, whether the particle is an electron (positron) or a pion, is based on an ar- 
tificial neural network implemented in ROOT [5]. After the particle identification, an event 
candidate must have at least two electrons and two positrons. Out of these tracks all possible 
pair-combinations are build, which can be the decay product of one of the two virtual photons. 
In the following, each combination is evaluated to be a double Dalitz decay event. 

To get rid of pions which have been misidentified as electrons (positrons), the opening angle 
between the supposed electron and the positron from the same virtual photon is checked. In case 
of an electron-positron pair it has to peak at small values. 

Another source of background are events, where a photon converts into an electron-positron 
pair. In case of the single Dalitz decay rj — > e+e~7, there is the same number of electrons and 
positrons in the final state as in case of the double Dalitz decay. To suppress such events, the 
distance between the point of closest approach of the positively and negatively charged particle and 
the beam center is analyzed as a function of the invariant mass of the lepton-pair, calculated from 
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Figure 1.1: Energy deposit in the central part of WAS A against the momentum of a particle times 
its charge. Left: Plastic Scintillator Barrel. Right: Electromagnetic Calorimeter. 



their four- vectors at the beam pipe. For conversion events, the invariant mass is zero, since the 
lepton-pair comes from a real photon. The radius of the closest approach peaks around 30 mm, 
since this is the radius of the beam pipe. The left panel of Fig. 1.2 shows the distribution for 
simulated conversion events and for signal events, and the right panel shows the distribution from 
the data. 
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Figure 1.2: Radius of the closest approach of the positively and negatively charged particle against 
the invariant mass calculated from their four-vectors at the beam pipe. Left: Simulation. Right: 
2009 data. The dashed line in denotes the used cut. 



After all these cuts the signal to background ratio for r\ events is ~ 2/1, where the main 
background channels are r\ — > 7r 7r 7r and rj — > 7r + 7r~7r°. 



1.3 Results 

Fig. 1.3 shows the development of the Tie-missing mass distribution for the 2009 data after each 
cut. The remaining background events are mainly from direct pion production. To get rid of those 
events the final Tie-missing mass distribution is fitted with a polynomial function for the direct 
pion background and a Gaussian function for the r\ signal. Fig. 1.4 shows the Tie-missing mass 
distribution from the 2008 and 2009 data after subtracting the direct pion background. After 
subtracting the background from ij decays besides the double Dalitz decay, 52 ± 13 events remain 
in the distribution. 

The combined result from the analysis presented in Ref. [4] and from the analysis, which is 
presented here, is BIZ = (3.0 ± 0.8 s tat ± 0.7 sys t( nO rm.)) x 10~ 5 (preliminary), which is in agreement 
with the value from Ref. [3]. The systematic error comes from the branching ratio calculation 
based on two different normalization channels, e.g. rj — > e + e~7 and rj — > 7r + 7r~e + e~. 
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Figure 1.3: Development of the ^e-missing 
mass distribution (2009 data). 
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Figure 1.4: Final 3 He-missing mass distribution 
after direct pion background subtraction (2008 
and 2009 data). 



1.4 Outlook 

As a next step the systematic error, which is approximately on the same order of magnitude as 
the statistical uncertainty, has to be studied in more detail. 

The statistics can be improved with pp —> pprj reactions, where WASA-at-COSY has already 
measured > 10 s t) mesons. On the one hand the production cross section is in this reaction around 
twenty times larger compared with the pd ^e 77 production reaction, but on the other hand 
there is also a larger contribution from direct pions. 
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DALITZ DECAY OF THE w MESON 



FARHA ANJUM KHAN* 
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1.1 Introduction 

Our motivation is to measure the electromagnetic transition form factor of the lo meson Dalitz 
decay lo — > 7r°e + e _ . The differential decay width at the vertex of the w-pion-transition can be 
described as a product of the differential decay width for point like - particles (a fuction of invariant 
mass of dilepton) and the transition form factor, where the form factor is the description of the 
electromagnetic structure arising at the transition point. The form factor is a function of q 2 where 
q is the transferred four momentum of the dilcptons equivalant to their invariant mass. In this 
Dalitz decay, the w vector meson decays into a pseudoscalar meson tt° and a dilepton. According 
to the VMD (Vector Meson Dominance) assumption this decay occurs via a virtual p meson, 
m w w m p = 0.77 GcV, thus a resonance at m* = \/\q 2 ) = m p . Experimentally, the form factor 
is determined by comparing the dilepton invariant mass spectrum with the point-like Quantum 
Electrodynamics prediction. 

The issue is that the transition form factor for the lo meson does not agree with VMD pre- 
dictions compared to other meson decays involving dileptons; for example, ir° — > 7°e+e~ and 
r) — > 7°1 + 1~ arc consistent with VMD predictions. However, even extensions of VMD do not seem 
to work for the lo meson. There is agreement when the decays of narrow light vector mesons into 
pseudoscalar mesons and dileptons are calculated to leading order by treating pseudoscalar and 
vector mesons on equal footing. These theoretical efforts can attempt to go beyond VMD in a sys- 
tematic way and are in need of further experimental input [1]. Experiments have been performed 
and results for the lo meson form factor measured from fi + p~ are found in the 1986 review by 
L.G.Landsberg [2] which shows disagreement with calculations. Recently, the NA60 collaboration 
has claimed to have confirmed the old data [3]. It would be useful to have data with e + e~ as an 
additional experimental approach giving access to smaller virtual photon masses including the full 
reconstruction of the decaying meson. 

Two sets of experiments have been performed with the WASA detector at COSY using two 
different reaction mechanism. The intention is to compare the quality of the data between pd 
and pp reactions, in the sense of a feasibility and background study for lo — >c + c~7r° decays. For 
p + d — »■ 3 Hc + lo 20 TB of data have been recorded during 12 days of beamtime at two different 
beam kinetic energies 1.5 GeV and 1.45 GeV. The excess energies are 88 MeV and 63 MeV, 
respectively. This reaction has a smaller cross section and smaller background than the p+p 
reaction. The cross section at 1.45 GeV is 83.6±1.5±2.2 nb [4]. The two beam energies have been 
selected to improve the background subtraction. The pilot experiment for the p + p — > p +p + lo 
reaction has been performed for 11 days at 2.063 GeV beam kinetic energy, 60 MeV excess energy 
and 18TB of data have been recorded. This reaction has a bigger cross section and significantly 
more background than p+d. The cross section for this reaction at 2.063 GeV of beam kinetic 
energy is 5.7±0.6±0.8±0.9^ib [5]. The bigger background makes the online and offline event 
selection more challenging. 

1.2 Preliminary analysis 

The first analysis is being performed using the p d induced reaction. A preselection, i.e. a subset 
of acquired data is used for the analysis. The preselection of the p+d reaction was based on the 
3 He selection by extracting those events where 3 He was formed. Particles emitted in the forward 
direction i.e. 3 He, can be detected in the forward detection system that includes tracking detectors 
and range hodoscopes using the AE-E method. The decay particles i.e. e + ,e _ ,7 can be identified 
in the central part of the detector. The central detector consists of a mini drift chamber operated 
in the field of a superconducting solenoid, a plastic scintillator barrel and an electromagnetic 
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calorimeter. We aim to analyse u> — > j*ir° — s>c + c~7r° but start with the real photon case w — > 771- 
which is one of our reference channels. For uj — » jtt — > 777, the tt° is reconstructed via the two 
decay photons, oj mesons can be tagged via the missing mass deduced from the 3 He detected 
in forward direction M w = (E p + Ed — Eue) 2 ~ (P P — Phc) 2 - The uj is also fully reconstructed 
with the invariant mass M 7r o 7 = (E„o + E 7 ) 2 — (P n a + P 7 ) 2 . A peak is expected at the uj mass 

0. 782 GeV/c 2 in both spectra. A two-dimensional spectrum is shown on the left panel of Figure 

1. In this histogram one can see an enhancement at the uj mass indicated by the black lines. The 
background comes from multi pion production. The projection of the spectra is shown on the 
right panel of Figure 1, one can see a peak at the uj mass (red line) for both energies of the p+d 
data set. In this histogram one can see an enhancement at the uj mass indicated by black lines. 
The black dashed line corresponds to 1.5 GeV/c 2 and the blue line correspond 1.45 GeV/c 2 beam 
kinetic energy. One can see the uj peak at the same point for both energies but with different 
background features. This will allow to have a better understanding of the background. 




Missing Mass ( 3 He) GeV/c 2 Missing Mass ( 3 He) GeV/c : 



Figure 1.1: The left panel shows the invariant mass vs missing mass plot. The intersection of the two 
black lines indicate the location of uj peak. The figure on the right panel is the projection for invariant 
masses above 0.5 GeV/c 2 . This figure shows the peak at 0.782 GeV/c 2 (red line) for both energies. The 
black dashes correspond to 1.5 GeV/c 2 and the blue histogram corresponds to the 1.45 GeV/c 2 data set, 
showing the two phase space regions and the different position of the peak relative to the background. 

For fast analysis we need preselected data for p+p reaction as well. The preselection has been 
started and data is saved into four different streams, using dedicated experiment triggers and 
extracting those events where at least two proton tracks were found. Central detector conditions 
for the four streams are: 

streaml: at least 3 neutral tracks in central detector (777°) 

stream2: at least 2 charged tracks and at least 2 neutral tracks in central detector 
(e + e _ 7r°, and other decays with similar topology) 

streamS: at least 2 charged in central detector 

stream^: exactly 2 neutral tracks in central detector (background studies) 

The goal of the analysis is to find the event candidates in the p+d reaction and compare with 
the p+p reaction and then decide which reaction can be used for the production beamtime. 

Bibliography 

[1] C. Terschluesen, St.Leupold, arXiv:1003.1030vl [hep-ph]. 
[2] L.G. Landsberg, Phys. Rept, 128, 301, (1985). 

[3] NA60 Collaboration (R. Arnaldi et al.), Phys. Lett. B, 677, 260, (2009). 

[4] K. Schonning et al. CELSIUS/WASA Collaboration, Phys. Rev. C, 79, 044002, (2009). 

[5] S. Barsov et al., Eur. Phys. J. A, 31, 95, (2007). 



53 



STUDY OF THE -> rje + e~ DECAY AT KLOE 

JAROSLAW ZDEBIK* 
for the KLOE-2 Collaboration 

1.1 Introduction 

The analysis of <f> — > rje + e~ decay is interesting from several points of view. The structure of cj> 
and rj mesons and underlying quark dynamics in the transition region can be extracted from e + e~ 
invariant mass spectrum. By comparing the experimentally measured spectrum of the lepton pair 
with QED calculations for pointlike particles, it is possible to determine transition form factor in 
the time-like region of momentum transfer [1]. 

The only one measurement of the form factor comes from SND collaboration is not in good 
agreement with predictions from the Vector Meson Dominance (VMD) framework [1, 2]. The form 
factor is often parametrized in one-pole approximation: 

W) = 1^7X3- (1-D 

where q = M ee and A is a free parameter. 

The theoretical calculation for A is 1.0 GeV (VMD), which the value measured by SND is 
0.5±0.1 GeV [3]. Also is interesting to compare experimental results with another theoretical 
models [4]. In this paper preliminary results of the investigation of the <p V e+e ~ decay at 
KLOE are presented. 

Additional goal of this work is also to understand the signature of this channel in the KLOE 
detector, since it is the main and nonreductive background for dark matter particles searches [5, 6]. 

1.2 The KLOE detector 

The KLOE (Klong Experiment) detector is installed at the interaction point of the electron 
and positron beams of the DA$NE (Double Annular ^-factory for Nice Experiments) collider 
operating in the Laboratori Nazionali di Frascati (LNF). 

DA$NE is an e + e~ collider running at a center of mass energy of ~ 1020 MeV (the mass of 
the <f> meson). The positron and electron beams collide at an angle of 7r-(25 mrad), producing 
<j> mesons nearly at rest. At the interaction point (IP) the beam pipe has the shape of a sphere 
which is made of a beryllium- aluminium alloy with 10 cm diameter and 50 fjxa thickness [6]. 

This detector was fully constructed by the end of the year 1998 [7]. It consists of two main 
subsystems: an electromagnetic calorimeter and a large drift chamber. The drift chamber and the 
calorimeter are inside a superconducting coil which produces a 0.57 T magnetic field parallel to 
the beam axis. 

Energy and time resolutions for calorimeter are ue/E — 5.7%/ \J~E (GeV) and 
a t — 57 ps/y/E (GeV) 100 ps, respectively. For the drift chamber, the spatial resolutions are 
<j xy ~ 150 /um and a z ~ 2 mm. The momentum resolution is a(p±)/p± w 0.4%. Vertices are 
reconstructed with a spatial resolution of ~ 3 mm. 
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1.3 Event Selection 

The analysis of the <f> — > rje + e~ decay with subsequent, r) — > ir + ir~ir a , has been performed on 
1.52 fb _1 of the KLOE dataset. The signal Monte Carlo (MC) simulation has been produced with 
dT((f> — > 7] e + e~~)/dq weighted according to Vector Meson Dominance model [1], using the form 
factor parametrization from the SND experiment [3]. Data-MC corrections for cluster energies 
and tracking efficiency, evaluated with radiative Bhabha events and <f> — > pn samples respectively, 
have been applied [6]. 

The first step of the analysis was preselection of events, that have to satisfy the following 
criteria: 

1. two positive and two negative tracks with point of closest approach to the beam line inside 
a cylinder around the interaction point (IP), with transverse radius R=4 cm and length 
Z=20 cm; 

2. two energy clusters in calorimeter with E > 7 MeV not associated to any track, in an angular 
acceptance |cos6> 7 | < 0.92 and in the expected time window for a photon (|T 7 — i? 7 /c| < 
MIN(5cr t ,2ns)); 

3. best 7r + 7r~77 match to the r\ mass with the pion hypothesis to assign ^ tracks; the other 
two tracks are then assigned to e ± ; 

4. loose cuts on r\ and tt° invariant masses (495 < M 7r+7r - 77 < 600 MeV, 70 < M 77 < 200 MeV). 

After this preselection, a clear peak corresponding to <fi — > rje + e~ events is observed in the 
distribution of the recoil mass to the e + e~ pair (Fig. 1.1). The second peak at ~ 590 MeV 
is due to 4> —¥ KsKl, K$ — > 7r + 7r~ events with a wrong mass assignment. Events in the 
535 < M reco ii(ee) < 560 MeV window are retained for further analysis. 
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Figure 1.1: Recoiling mass against the e + e pair for the data sample and MC events after prese- 
lection. 

The next selection level was performed taking into account the information from the signal 
channel MC simulation. After the preselection, the remaining background contamination mainly 
comes from <fi — ¥ 777 reaction and events which have more than two charged pions in the final state. 

The background, due to <f> — > 777 events with photon conversion on beam pipe (BP) or drift 
chamber walls (DCW), is rejected by tracing back the tracks of the e + , e~ candidates and 
reconstructing their invariant mass (M ee ) and distance (D ee ) at the BP/DCW surfaces. As 
both quantities are small in case of photon conversions, <j) — ¥ rpf background is removed by 
rejecting events with: M ee (BP) < 10 MeV and D ee (BP) < 2 cm, M ee (DCW) <_80 MeV and 
D ee (DCW) < 10 cm. The second relevant background, originated from <p —> KK decays sur- 
viving analysis cuts, has more than two charged pions in the final state and is suppressed using 
time-of-flight (ToF) to the calorimeter. When an energy cluster is connected to a track, the ar- 
rival time to the calorimeter is evaluated both using the calorimeter timing (Truster) and the track 
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trajectory (T trac k = £track//3c). The AT = T trac k — r c i us t er variable is then evaluated for both 
electron (AT e ) and pion (AT ff ) mass hypotheses. Events with an e + , e~ candidate outside a 3cr's 
window on the AT e variables are rejected. The distribution ol invariant mass of e~e + pair after 
all analysis cuts is reported in Fig. 1.2. 
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Figure 1.2: Distribution of the e + e invariant mass spectra for the process <fi — > rje + e with 
rj — > Tr + TT~ir° decay chain. 

After all selection cuts the contamination coming from background events is less than 5% in the 
final sample and more than 14000 events of <fi — > rie + e~ reaction we reconstructed. The statistics 
is almost two orders of magnitude larger than in any previous measurment. 



1.4 Preliminary Results 

A fit to the invariant mass distribution corrected with acceptance and background subtracted was 
done using parametrizantion from Ref. [1]. 



dT{(j)^r]e + e-) 
dq 2 

a \F H (g 2 )\ 2 
3tt q 2 



/ Am 2 ( 2m 2 
1 s- I 1 + — s- 



m 2 - m 2 



m 2 , — m 2 



(1.2) 



Free parameters of the fit are A (reported in Eq. 1) and an overall normalization factor. The 
fit result is shown in Fig. 1.3. 
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Figure 1.3: Fit to the e + e~ invariant mass distribution corrected with acceptance and background 
subtracted. 

The achieved statistical error from the fit for A is about 0.005 GeV. The evaluation of systematic 
errors and implementation of the smearing matrix into fit procedure is in progress. 
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ELECTROMAGNETIC TRANSITION 
FORM FACTORS 
OF PSEUDOSCALAR AND VECTOR MESONS 

CARL A TERSCHL USEN*^ and STEFAN LEUPOLD* 

1.1 Introduction 

The aim of this contribution is the description of two kinds of decays, 

• light vector meson (p, to, K*, <fr) into a light pseudoscalar meson (n, K, 77, rj') and a dilepton, 

• pseudoscalar meson into a photon and a dilepton. 

Both processes might contain intermediate vector mesons. 

The most general transition form factor for a decay including a virtual vector meson with mass 
m is given as 

2 94 

Tfi Q q 

Fgcn.(q) = 90 — o + ( X ~ 90) + 91 — o + 92 —7 + ■ ■ ■ (1.1) 

m z - q z tji to 4 

with g 2 denoting the square of the invariant mass of the dilepton. In general, this form factor is an 
infinite series. Therefore, the parameters cannot be fixed. Progress can be made if it is possible to 
introduce a power counting scheme which orders the coupling constants go, gi, 92, ■ ■ ■ according 
to their importance. Therewith, it is possible to perform calculations up to a given order with a 
finite number of parameters. Additionally, one will be able to improve the results systematically 
by adding the next order. 

Due to the running coupling constant in QCD, perturbation theory can only be used for high 
but not for low energies. A possible solution are effective theories which take hadrons instead 
of quarks as relevant degrees of freedom. The effective theory called Chiral perturbation theory 
(ChPT) [1] takes the light pseudoscalar mesons as relevant degrees of freedom and treats all other 
mesons, in particular the vector mesons as heavy. Therefore, it is not applicable for the energy 
range of the hadronic resonances. 

We are using a new counting scheme [2] which treats the masses of both light vector and pseu- 
doscalar mesons as soft, i.e. of the order of a typical momentum q. For decays, all involved 
momenta are smaller than the mass of the decaying particle and, hence, also of the order of q, 

m v , to p , ~ q. (1.2) 

In ChPT, the range of applicability, i.e. the range for q is limited (on tree level) by the not- 
considered mesons, in practice by m v , and (for loops) by the scale 47r/, where / denotes the 
pion decay constant. In the scheme of [2], where vector mesons are included and where two- 
particle reducible diagrams (rescattering processes) are resummed, it is suggestive that the range 
of applicability can be pushed to larger energies. Next-to-leading-order calculations are necessary 
to assess this proposition. At present, phcnomenological consequences are worked out in leading 
order and compared to data. 

1.2 Decay of a light vector meson into a pseudoscalar meson 
and a dilepton 

Using the counting scheme (1.2), one can determine the leading-order Lagrangian for the decay of 
a vector meson into a pseudoscalar meson and a (real or virtual) photon, 

* Uppsala University 
tcarla.terschluesenaphysics.uu.se 



58 



C. Terschliisen 



PrimcNet Workshop 2011, Jiilich 



16/ 

eymy 



tr{V^Q}d^A u . 



(1.3) 



Thereby, describes the vector mesons in antisymmetric tensor representation and $ the pseu- 
doscalar mesons, 



:= 



V2p+ 



-P° 



\/2tt- 
V2K- 



V2ir+ 



— TV 



7^ m 
V2K° 



V2K+ U 

V2K+ \ 
V2K° 

) 



(1.4) 



(1.5) 



Furthermore, xo = diag {m 2 , m 2 , m%} , Q — diag {2/3, — 1/3, — 1/3} and A v denotes the photon 
field. Note that strictly speaking the non-Goldstone boson r\\ was not included in the scheme of 
[2]- 

The first two terms in (1.3) describe the decay of a vector meson into a pseudoscalar meson and 
a virtual vector meson, the last term the decay of the virtual vector meson into a photon. The 
decay of the photon into a dilepton is described by usual QED. Obviously, the counting scheme 
(1.2) allows in leading order only for decays via virtual mesons. Additionally, the parameters h A 
and b A can be fixed by comparing the calculated results of two-body decays into real photons 
to experimental data. There are no additional parameters needed to describe the decays into 
dileptons so that our calculations for such decays a real predictions. 

If the decay can happen only via one type of virtual vector meson with mass m, the standard 
vector-meson dominance (VMD) form factor equals 



Fvmd(q) = 



rrr 



m 2 — q 2 



Our form factor calculated with (1.3) gets an additional non-VMD term, 

^,2 



F(g) = go 



m 2 — q 2 



+ (1 - <?o) • 



(1.6) 



(1.7) 



On the left-hand side of Fig. 1.1, the uj — > 7r° transition form factor is plotted [3] (solid 
line) in comparison to the standard VMD form factor (dashed line) and data taken by the NA60 
collaboration for the decay into a dimuon [4]. Obviously, the standard VMD form factor fails to 
describe the data whereas our calculation misses only the last three data points. In the middle of 
Fig. 1.1, the </>—»• 77 transition form factor is plotted [3] (solid line). Hereby, the physical r\ meson 
is approximated by the octet state 77s- Again, it is compared to the standard VMD form factor 
(dashed line) and data taken at the VEPP-2M detector for the decay into a dielectron [5]. Due 
to the relatively large error bars there is no assessment possible which form factor describes the 
data better. In the near future better data are expected form KLOE 1 . 



1.3 Decay of a light pseudoscalar meson into a photon and 
a dilepton 

As a more phenomenological approach for the decay of a light pseudoscalar meson into two (real 
or virtual) photons, we take both our Lagrangian (1.3) describing the decay via two virtual vector 
mesons and the Wess-Zumino-Witten Lagrangian [1] 

Avzw - e " VaP tr {<3 2 *} 9,A a dvAp + 0($ 2 ) (1.8) 

1 See the contribution by Jaroslaw Zdcbik in these proceedings. 



•59 



PrimcNct Workshop 2011, Jiilich 



C. Terschliisen 



which describes the direct decay into two photons. Additionally, the r\-rj mixing is given by 

rj = cos Oris — sm Orji A rj' = sin 9r]$ + cos 9 r\\ (1.9) 

with the mixing angle 9 w —19.5°. 

The rj — > 7 transition form factor is plotted on the right-hand side of Fig. 1.1 (solid line) 
in comparison to the standard VMD form factor 2 (dashed line) and data taken by the NA60 
collaboration for the decay into a dimuon [4] . Here, the form factors arc on top of each other and 
describe the data equally well. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 0.6 

q [GeV] q [GeV] q [GeV] 

Figure 1.1: Left-hand side: to — > tt form factor compared to dimuon data taken by the NA60 
collaboration [4]. Middle: (j> — >• 77 form factor compared to dielectron data taken at the VEPP-2M 
detector [5]. Right-hand side: 77 — > 7 form factor compared to dimuon data taken by the NA60 
collaboration [4]. The first two results are published in [3]. 



1.4 Summary and outlook 

Using a new counting scheme (1.2) we were able to determine the leading-order Lagrangian for 
the decays V — > Pl + l~ and for the vector part of the decay P — > The form factors for 

the transitions u> — > ir°, <f> — > r\ and r\ — > 7 describe the experimental data very well, especially 
the w — > tt transition is much better described than with standard VMD. Additionally, we were 
able to describe the partial decay widths very well [3]. As a next step, next-to-leading order 
calculations have to be performed. 
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WITH WASA-AT-COSY 

LENA HEIJKENSKJOLD* AND SIDDESH SAWANT* 
for the WASA-at-COSY Collaboration 

1.1 Introduction 

The WASA-at-COSY collaboration has recently performed experiments dedicated to studies of 
the uo meson decays. One of the decay channels that will be under detailed investigations is the 
most probable decay of the uo meson into tt + tt^tt°. The goal of this experiment is the investigation 
of the 10 decay mechanism to three pions. This information may be gained from a comparison of 
the density distribution in Dalitz plot with appropriate theoretical predictions. To reach this goal 
a high statistics experimental Dalitz plot distribution is necessary. Measurements of the Dalitz 
plot for co — »■ 7T + 7r~7r° decay have been done before in experiments on the determination of the 
spin and parity of the uo meson [1]. However, the statistics of existing measurements is limited to 
about 4600 uj — > 3tt events. The WASA-at-COSY experiment has produced a higher amount of 
fully reconstructed lo — > 37r events using proton-proton and proton-deuteron collisions. 

The decay mechanism was considered within Vector Meson Dominance model as proceeding via 
the pn intermediate state [2]. More recent predictions base on a counting scheme for flavor-SU(3) 
systems of Goldstone bosons and light vector mesons [3]. The importance of the irn final state 
interaction is also considered [4]. However, the expected deviations from Vector Meson Dominance 
model are small. 

1.2 WASA-at-COSY 

The experiments have been performed with the WAS A detector setup at COSY [5]. COSY (COoler 
SYnchrotron) is a synchrotron and storage ring which can provide a beam of (un)polarized protons 
or deuterons with a momentum range of 600-3700 MeV/c. 

WASA (Wide Angle Shower Apparatus) is a nearly 4-7T detector optimised for studying pro- 
duction and decay of light mesons. A unique pellet generator provides frozen pellets of hydrogen 
and deuterium as a fixed target. The central part of WASA can detect both charged and neutral 
particles while the forward detector is designed to detect scattered projectiles and charged recoil 
particles like protons, deuterons and He nuclei. 

1.3 The current analyses 

In spring 2011, WASA-at-COSY collected data of co production in proton beam on proton target 
(p-p) collisions (p + p — > p + p + cj) at T p =2.063 GeV and in proton beam on deuteron target (p-d) 
collisions (p + d — > 3 Hc + u;) at T p =1.45 GeV and T p =1.5 GeV. The corresponding excess energies 
are 60 MeV for the p-p reaction and 63 MeV and 88 MeV for the p-d reaction, respectively. These 
two reaction have complementary advantages and drawbacks, which are summarised in Tab. 1.1. 

1.3.1 First preliminary results 

Preliminary analyses have been performed on both data sets. In these analyses the protons or the 
3 He particles are identified in the forward part of WASA using the AE/E technique. The decay 
products of the uo are detected in the central part of WASA. For charged pions, the charge and 
momentum are reconstructed using information from the central mini drift chamber placed in the 
solenoid magnetic field. The n° is reconstructed using the energy deposit and angle in the central 
calorimeter of the two photons. 

In presented analysis 11% for the p-d data and 15% for the p-p data of the detected events have 
been analysed. On the left panel of Fig. 1.1, which shows the missing mass of the 3 He detected 
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Table 1.1: The table displays the advantages and drawbacks of the two production reactions. 



Reaction 


Advantage 


Drawback 


P-P 


Higher production cross section, 

at T p = 2.063 GeV er tot = 5.7 ^b [6]. 


Fast protons, majority punches 
through forward part of detector. 
Requires elaborate trigger. 


p-d 


Simple triggering on 3 He. 

3 He stops in forward part of detector. 

Advantageous for tagging. 


Lower production cross section, 

at T p = 1.45 GeV a to t = 83.6 nb [7]. 



for p + d — s- 3 He + tt + + ir~ + ir°, a peak at the value of to mass is clearly visible. In the middle 
plot of Fig. 1.1 the three pion invariant mass versus the missing mass of two protons detected for 
p + p— >p + p + 7r + + 7r~ + 7T° reaction is presented. By putting a cut shown in the figure on the 
invariant mass of three pions around the u meson mass, the missing mass of two protons (right 
plot in Fig. 1.1) does not show a clear peak but only a shoulder at the value of cj mass. This 
suggests a small uj meson signal hidden under the three pion background. 
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Figure 1.1: Left: The missing mass deduced from the 3 He in the p-d data at T p = 1.5 GeV. 
Middle: The invariant mass of three pions versus the missing mass from the protons in the p-p 
data. Black lines show the cut on the invariant mass of 7r + 7r~7r°. Right: The missing mass from 
the protons in the p-p data after the cut on the invariant mass of the three pions. 



1.3.2 Outlook 

In order to extract the u signal, the next steps in the analyses would be fine tuning of energy 
calibration of all detectors. Then the kinematic fitting method will be applied to improve overall 
resolution. These steps should allow to obtain the best signal to background ratio and properly 
select the to signal. 
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U)/<j) — > 37T AND T}' TJ7T7T 
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1.1 Introduction 

The advent of high-statistics measurements of hadronic three-body decays calls for equally accu- 
rate theoretical descriptions. In this context rescattering effects between decay products play a 
very important role if one aspires to perform precision amplitude analyses. While perturbative 
approaches, such as chiral perturbation theory or non-relativistic approaches, implement final- 
state interactions up to a certain order in a small power-counting parameter and are restricted 
to a low-energy regime, dispersive approaches resum hadronic rescattering up to all orders and 
thus allow for an extension to higher energies. While the groundwork for dispersive analyses of 
hadronic three-body decays has been laid long ago [1], high-accuracy parameterizations of phase- 
shifts and sufficient computational power, which are necessary for such an enterprise, have only 
recently become available. 

We apply the dispersive approach to two specific examples. uj/<f> — > 'Stt is the simplest imag- 
inable system in which three-body effects in hadronic rescattering can be studied, due to the fact 
that the decaying particles are of vector type. If one neglects contributions from F- and higher 
partial waves, the system is completely determined by P-wave interactions. In the case of <j> Sir 
the phase space is sufficiently large to study three-body effects on the p-resonance. With its large 
existing (<f>: KLOE/CMD-2 [2, 3]) and upcoming (w: WASA-at-COSY) data base u/4> -> 3tt 
provides an ideal testing ground for the approach. 

77' — > r]nn on the other hand is one of the very few processes in which effects from nri scattering 
might be studied, in the hope of gathering information on 71-77 phase shifts. Additionally, the 
neutral channel r/' — > r]n°n° exhibits a cusp effect that may allow for a determination of tttt 
S-wave scattering lengths [4]. Finally, recent and upcoming experimental efforts from various 
collaborations (BES-III, WASA-at-COSY, ELSA, CB-AT-MAMI-C) call for a renewed analysis. 

1.2 The framework 

The goal of the dispersive approach is to relate the decay amplitudes to a set of integral equations 
derived from the fundamental principles of unitarity and maximal analyticity, which are then solved 
numerically. We follow a Khuri-Treiman type dispersion relation ansatz [1], where the dispersive 
result for two-body scattering is linked to the three-body decay by analytic continuation in the 
mass of the decaying particle and the center-of-mass energy. 

1.2.1 Decay amplitude for uj/4> — >■ 7r°7r + 7r _ 

In the following we resort to the standard Mandelstam representation for three-body decays, 
s = (p-k- + fV+)> t = (p n - +?V)) 2 , u = (7^+ +p. n °) 2 with s + t + u = 3s = M v + 3M%, where 
My = M w /0 is the mass of the decaying vector meson. The process is of odd intrinsic parity and 
we may decompose the amplitude A4(s, t, u) according to 

M(s,t,u) = e )iva pn IJ -p^ + p'Z-pP J r {s,t,u) , (1.1) 

where n M is the polarization vector of the decaying u /<j). Due to Bose symmetry only odd angular 
momenta contribute in the partial- wave expansion of IF{s, t, u). We restrict our analysis to P-wave 
contributions only, since the effects of F- and higher partial waves are expected to be significantly 
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suppressed. The P-wave part of the amplitude can be decomposed into a sum of single-variable 
amplitudes T(s) (compare [5]), 

T(s,t,u) = J*(s) + F{t) + F{u) . (1.2) 
Imposing unitarity and elastic final-state interactions provides us with the discontinuity of J"(s) [6], 

discF(s) = 2i[F(s)+J r (s)]e- lSl( - s hm(5 1 (s)) , (1.3) 

where J r (s)+J r (s) is the P-wave projection of J"(s, t, u) and 8\{s) the irir P-wave phase shift taken 
from phenomenology [7, 8]. The inhomogeneity -F(s) is given as 

f(s) = 3((1 - z 2 ) ?) , (z n f) =\fdz z n f (^- s + z < s ) 



2 J-i V 2 



«(*) = yjl- ^V((Mv + M n y - 7) pv - M w )2 - s ) . (1.4) 

To solve for the amplitude we use a product ansatz J-"(s) = fl(s)<p(s) where f2(s) is the well-known 
Muskhelishvili-Omnes (MO) solution to the homogeneous equation J"(s) = [9]. For J"(s) we 
thus obtain 

-7-/ \ , J * /°° JXf)sirW)) ) m . (a f°° ds' 5(s>) \ 

The subtraction constant a is fixed by the partial decay width. The number of subtractions is 
chosen such that the dispersion integral is guaranteed to converge. The issue of performing the 
double integration correctly in view of the analytic continuation in My and s will not be discussed 
here. We refer to [10] and references therein. 

1.2.2 rj' — > 777777 

Similarly to the co/(f) — >• 37r case we can decompose the full amplitude into isospin components, 

M(s,t,u) = MZ*{s)+Ml r] {t) + MZ r, (u) , (1.6) 

with the amplitudes of the 7T7t S-wave and the ttij S-wave Mq V . We neglect the exotic iri] 

P-wave throughout. The set of integral equations can be obtained analogously as before, 



ds' sin£J' r (s / ).Mo 7r ( s/ ) 



f s 3 rc 

MV(s) = fTO{a„ + „„, + + - / 4m; s , 3 |nr(( , )|(s , _ s _ u) 

^ w - gi»f, 4 £ g |$^' ( rfi'2) }■ <, = (m,+^) 2 . m 

We fix the 7r?7 phase shift i5^ n using a simplified unitarization model. The subtraction constants are 
matched to Resonance Chiral Theory (RChT, see [11]), using a narrow- width-resonance approxi- 
mation for the Omnees function O(s) — > M|/(M| — s). This is by no means a perfect method to 
fix the subtraction constants, since obviously the rnr S-wave is not well-described by a resonance 
approximation. An alternative would be to match to large- Nc chiral perturbation theory. In 
this case, however, it is not clear how to properly treat resonance contributions in the low-energy 
constants. 



1.3 Preliminary numerical results 

The integral equations are solved numerically by an iterative procedure. We start out with an 
arbitrary input function and calculate the inhomogeneity J-(s). The single- variable amplitude is 
then calculated from this inhomogeneity, which is in turn used as input for the following iteration 
step. We repeat the procedure until convergence is reached up to satisfactory accuracy. 
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1.3.1 (f>/<*> — > 37r amplitudes and Dalitz plot 

To illustrate the convergence of the iteration procedure the <j> — > Ztt amplitude T(s) is plotted 
after each iteration step in Fig. 1.1. We have chosen the MO function as the starting point for 
the iteration. It should be mentioned that the final result of the iteration is independent of the 
specific choice of the starting point; however, choosing the MO function allows us to quantify 
crossed-channel rescattering effects in a plausible way. We find that these effects are sizable 
and from the point of view of precision analyses certainly not negligible. As seen in Fig. 1.1, 
convergence of the amplitude is reached after three iterations; for uj — > 37r (not shown here), this 
is the case after two iterations. 

For a three-body decay, the Dalitz plot distribution is of particular interest. The Dalitz plot is 
a two-dimensional scatter plot in terms of two kinematic variables; in the case of (j)/w 37r it is 
common practice to use x = (t — u)/(2My) and y = [(My — M^) 2 — s]/(2My). The normalized 
Dalitz plot is shown in Fig. 1.2; we note that it is independent of the subtraction constant, so that 
Fig. 1.2 is free from any input aside from the irir P-wave phase shift, which is well-established up 
to roughly 1.2 GcV. 

The uj — ¥ 37r Dalitz plot is a relatively smooth distribution, which rises from the center to its 
border. The outer corners show a maximum increase of roughly 20% with respect to the center. 
In contrast, the cj> 3 n Dalitz plot exhibits significantly more structure. From its center the 
Dalitz plot rises towards the bands created by the p resonance, and then steeply falls off towards 
the outer corners. The profile of the — > 3ir Dalitz plot seems to agree qualitatively with the one 
shown in [2]; for an in-depth analysis, closer comparison to the data is certainly required. 

Finally we wish to study the effects of crossed-channel rescattering on the Dalitz plot distribu- 
tion. For that purpose, we follow two different approaches to fix the subtraction constant. First we 
assume that the subtraction constant is given by some independent method, and we are interested 
in what bearings the crossed-channel effects have on both the overall shape of the Dalitz plot 




-200 -150 -100 -50 50 100 150 200 -300 -200 -100 100 200 300 



x [MeV] x [MeV] 

Figure 1.2: Normalized Dalitz plot for u> — > 3tt (left) and 4> — > 3tt (right). 
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Figure 1.3: l-^fuiil 2 / 1-^=0 1 2 ^ or u ~~ ^ ^ 7r O e ft) an d 4> ~ * 3tt (right). The subtraction constant is 
fixed so that it reproduces the decay width before and after the iteration. 

and the partial decay width. We fix the subtraction constant from the experimental decay width 
for the case T = 0. The effect on the decay width in both processes amounts to roughly 20%. 
However, while the partial width of u — > 3tt sees an increase, it suffers a decrease in <f) — > 37r. The 
same qualitative behaviour shows up in the Dalitz plot distribution. 

In our second approach, we match the subtraction constant to reproduce the experimental 
decay width in both cases, that is with and without crossed-channel rescattering effects included. 
We then study the effects on the profile of the Dalitz plot. The result of this approach is shown for 
both u) — > 3tt and <fi — > 3ir in Fig. 1.3. One observers that the effects of crossed-channel rescattering 
are to a large degree absorbed in the partial width. The overall behavior is very similar: in both 
processes one observes an increase from the outer corners to the center of the Dalitz plot. It is 
noteworthy that the p resonance bands in (f> — > 3n are left relatively untouched. 

For a rough error estimate we have varied the cutoff in both the integral of the MO function 
and the dispersion integral. Further we have varied between two sets of phenomenological repre- 
sentations of the 7T7T P-wave phase shifts [7, 8] . We found that these effects are small and certainly 
not comparable to the sizable crossed-channel effects. 

1.3.2 r]' — > TjiziT Dalitz plot 

The rj — > 777T7T amplitude exhibits a similar convergence behavior as (j) — > 37r. Crossed-channel 
effects are sizable and convergence is reached only after three iterations. It is also interesting to 
study the effects of crossed-channel rescattering on the Dalitz plot, more specifically the Dalitz 
plot parameters. Slightly different conventions for the Dalitz plot variables in rj — > r/mr as opposed 
to u>/(f> — > Sir are used, namely 

= (M v + 2M„)(M ri , -M v f -a _ = V3(t - u) 

V 2M V 'M W (M V > - M n - 2M„) 1 X 2M V >(M V > — M v — 2M„) ' [ '' 

where s and t are defined analogously as before. The amplitude squared can then be expanded 
around the center of the Dalitz plot in terms of x and y according to 

\M(x,y)\ 2 = \N\ 2 {l + ay + by 2 + dx 2 } , (1.9) 

neglecting terms of cubic order and higher, where Af is the overall normalization and a, b and 
d are the Dalitz plot parameters. A term oc x is forbidden by C-parity. Table 1.1 shows the 
parameters as determined by our procedure. The following comments are in order: we matched 





RChT [11] 


M = 


full 


VES Collaboration [12] 


a 


-0.119 


- 0.239 


- 0.478 


-0.127+ 0.016 + 0.008 


b 


+ 0.001 


- 0.025 


+ 0.053 


-0.106+ 0.028 + 0.014 


d 


-0.056 


- 0.023 


- 0.043 


-0.082 + 0.017+ 0.008 



Table 1.1: Determination of the rj — > -qnir Dalitz plot parameters by RChT, the pure MO solution 
{M. = 0), the full dispersive result and the experimental measurement by the VES collaboration. 
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the dispersive representation to RChT in the limit fl(s) — > Mg/(Mg — s). This is not fulfilled well 
by the irir S-wave, which manifests itself in large modifications of the parameter a from the RChT 
result to the MO solution. From that vantage point, if the RChT result is already in marginal 
agreement with experiment, one cannot expect the outcome of the full analysis to be close to 
the measured parameters, especially if the crossed-channel rescattering corrections are as large 
as evidenced in the table. Certainly the matching procedure and thus the determination of the 
subtraction constants is still the weakest link in our analysis and requires further improvement. 

Finally, while we do observe effects of nrj phase shifts on the Dalitz plot parameter d, these are 
not large enough to extract precision information about nr] scattering, especially if the subtraction 
constants are not under control. 
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HARD PION CHIRAL PERTURBATION THEORY: 

WHAT IS IT AND IS IT RELEVANT FOR 77' DECAYS? 

JOHAN BIJNENS* 

1.1 Introduction 

In this talk I will try to convince you that we can give predictions from chiral symmetry also for 
cases where not all pions are soft. This is something I called hard pion Chiral Perturbation Theory 
(HPChPT) and there have been a few recent papers using this [1, 2, 3, 4, 5]. 

I will first give a short introduction to effective field theory (EFT) and remind you of the 
underlying principles of Chiral Perturbation Theory (ChPT). I will remind you of the fact that in 
ChPT with baryons and other heavy particles a power-counting has been achieved by consistently 
absorbing the heavy mass dependence into the low-energy-constants (LECs). 

The arguments will then be generalized to the case of processes with high energy or hard pions. 
The arguments also apply to cases where we can treat the strange quark mass as small as well. 

After that I will show applications to K — » -kit, to scmilcptonic decays of pseudo-scalar mesons 
or to more general vector form-factors and to charmonium decays to two pseudo-scalars. 

Unfortunately, there seem to be no rf decays where the present method are applicable. 

1.2 Effective field theory and ChPT 

The underlying idea of EFT is a general theme in science, restrict yourself to the relevant degrees 
of freedom. So, in cases where there exists an energy or mass gap we keep only the lower degrees 
of freedom. Lorcntz-invariance and quantum mechanics imply that we are restricted to a field 
theory but we should build the most general one with our chosen degrees of freedom. We have no 
predictability left since the most general Lagrangian will have an infinite number of parameters. 
This can be cured if we find an ordering principle, power-counting, for the importance of terms. 

ChPT is "exploring the consequences of the chiral symmetry of QCD and its spontaneous 
breaking using effective field theory techniques" and was introduced as an EFT in [6, 7]. The 
degrees of freedom are the Goldstonc bosons from the spontaneous breakdown of the global chiral 
SU(n) l x SU {ti)r to the diagonal vector subgroup SU (n)y . That the Goldstone boson interactions 
vanish at zero momentum allows to construct a consistent power- counting [6] . 

The basic form of ChPT has since been extended to baryons, mesons and baryons containing a 
heavy quark, vector mesons, structure functions and related quantities as well as beyond the pure 
strong interaction by including weak and electromagnetic internal interactions. Many models of 
alternative Higgs sectors also use the same technology. 

1.3 Power- counting and one large scale 

In purely mesonic ChPT the power-counting is essentially dimensional counting and this works 
since all the lines in all diagrams have "small" momenta. Already when discussing baryons, this 
lead to problems because now there is a large scale, the baryon mass. However, by setting the 
baryon momentum ps — Mbv + k with v the baryon four-velocity, a consistent power-counting 
can be achieved. This works "obviously" since the heavy line goes through the entire diagram 
and all momenta apart from Mg« are soft as indicated by the thick line in Fig. 1.1(a). The same 
arguments apply to ChPT for mesons containing a heavy quark. It works because the "soft" stuff 
can be expanded in "soft/Ms" and the remaining Mb dependence can be absorbed in LECs. 
For vector meson ChPT there is a problem since they can decay. A typical diagram is shown in 
Fig. 1.1(b). However it was argued that the non-analytic dependence on the light quark mass 
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Figure 1.1: (a) A typical baryon ChPT diagram with the baryon going through the entire diagram, 
(b) An example of a diagram in vector meson ChPT with no continuous vector meson line. 




Figure 1.2: The process of cutting the soft lines and reproducing the non-analytic dependence by 
the diagram on the right. The hard lines are thick, soft lines are shown thin. 

could still be obtained, see the discussions in [8] . Again, the underlying idea is that the large My 
allows to expand in "soft/My" and the remaining My dependence is absorbed in the LECs. 

1.4 Several large scales or HPChPT 

In [1] the authors applied heavy Kaon ChPT to at the endpoint, i.e. the pion is soft, this 
works as in usual ChPT. They also applied it to the region for small q 2 where the pion has a large 
momentum and gave arguments based on partial integrations why this would give a correct chiral 
logarithm. The argument was generalized in [2, 3, 4, 5]. The underlying idea is similar to the 
previous section. The "heavy /fast /hard" dependence on the soft stuff can always be expanded 
and the remaining dependence goes into the LECs. That this might be possible follows also 
from current algebra. Non-analyticities in the light masses come from the soft lines and soft 
pion couplings are restricted by current algebra via \iiR q ^ (ir k (q)a\0\(3) = — j^(a\ [Q§,0] \/3) . 
Nothing prevents hard pions to be in the states a or (3, so by heavily using current algebra one 
can get the light quark mass non-analytic dependence 

A field theoretic argument is: (1) Take a diagram with a given external and internal momentum 
configuration. (2) Identify the soft lines and cut them. (3) The resulting part is analytic in the soft 
stuff, so it can be described by an effective Lagrangian with coupling constants dependent on the 
external given momenta (Weinberg's folklore theorem [6]). (4) The non-analytic dependence on 
the soft stuff is reproduced by loops in the latter Lagrangian. The process is depicted in Fig. 1.2. 

The remaining problem is that we have no power-counting. The Lagrangian that reproduces 
the non-analyticities is fully general. In the HPChPT papers it was shown that for the processes 
at hand, all higher order terms can be reduced to those with the fewest derivatives thus allowing 
the light quark mass chiral logarithm to be predicted. The underlying arguments were tested by 
comparing to a two-loop calculation [4] and by explicitly keeping some higher order terms [2, 3, 
4,5]. 

1.5 Applications 

We have applied the method to K — > irir decays [2] where we treat the Kaon as heavy and look 
for the dependence on the pion mass M 2 . The result is, up to linear in M 2 and higher order: 

Ag L °/A%° = 1 + (3/8)1 , A% LO /A%° = 1 + (15/8)1 , A = -M 2 la(M 2 /fx 2 )/(Ww 2 F 2 ). 

The scalar and vector form- factors of the pion are known to two-loops in ChPT [9]. HPChPT 
predicts at large t [4] with Fy(t,0) and Fs(t,0), the form-factors at large t in the chiral limit, 
completely free: 

F v (t,M 2 ) = F v (t,0)(l + A), F s (t,M 2 )=F s (t,0)(l + (5/2)A). 
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Figure 1.3: The CLEO data on D —> ir and D — >• K. Form- factors as measured to the the left and 
corrected with the chiral logarithms to the right. Note the improved agreement. 

The full two-loop result expanded for large t should have this form and it does with for e.g. Fy 
F v (t,Q) = 1 + (t/(167r 2 F 2 )) (5/18- 16tt 2 ^' + i7r/6 - (1/6) ln(i/> 2 )) . 

The first application was semileptonic form- factors in K^. We extended this to -> 
D, 7r, K, ^-decays in [3, 4]. This allowed to test our results experimentally. The form factor f+(t) 
as measured by CLEO [10] in D —> it and D — >• K decays are different by about the amount 
expected from the chiral logarithms as shown in Fig. 1.3. One puzzling observation [4] was that in 
the limit of a hard pseudo-scalar in the final state the correction was always the same for / + and 
/_. This is in fact due to the LEET relation [11] which shows that there is only one form- factor 
in this limit and it is nice to see that our calculation respects this without it being used as input. 

A last application was to Xco.2 ~> ■K 1 KK,nn [5]. Here it was found that there were no chiral 
logarithms to the order considered. Comparing with the known experimental results indeed shows 
SU(3) breaking to be somewhat smaller than in e.g. Fk/F^. Details can be found in [5]. 
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HUNTING RESONANCE POLES 

PERE MASJUAN* 



1.1 Introduction 

The non-perturbative regime of QCD is characterized by the presence of physical resonances, 
complex poles of the amplitude in the transferred energy at higher (instead of the physical one) 
complex Ricmann sheets. From the experimental point of view, one can obtain information about 
the spectral function of the amplitude through the Minkowsky region (q 2 > 0) and also about its 
low energy region through the experimental data on the Euclidean region (q 2 < 0). 

In reference [1], the particular case of the tttt Vector Form Factor (VFF) was analyzed using 
the available Euclidean data and its first and the second derivatives were determined at q 2 = 
with Pade Approximants (PA) centered at the origin trough a fit procedure to that data [1]. In 
such a way, the vector quadratic radius (r 2 ) v and the curvature c v were extracted from the fit 
and, as a consequence, a value for the low-energy constant L 9 = (6.84 ± 0.07) • 10~ 3 was obtained, 
[1, 2]. Similar ideas are applied in reference [3] to study the pion Transition Form Factor (TFF) 
and its first (a n ) and second (b n ) derivatives with better precision compared to the current world 
average result [4]. The implications of these parameters on the tt° contribution to the hadronic 
light-by-light scattering to the muon anomalous magnetic moment are also studied showing the 
relevance of this kind of model-independent methods. 

Despite the nice convergence and the systematical treatment of the errors, this procedure does 
not allow us to obtain properties of the amplitude above the threshold, such as in the case of the 
tttt VFF, the p-meson pole position, or trough the tttt phase shift the cr-pole. The reason is simple: 
the convergence of a sequence of Pade Approximants centered at the origin of energies (q 2 = 0) is 
limited by the presence of the tt — tt production brunch cut: the PA sequence converge everywhere 
except on the cut [2]. Still, the mathematical Pade Theory, through the Montessus de Ballore 
theorem, allow us to produce a model independent determination of resonance poles when certain 
conditions are fulfilled. The most important one is to center our PA sequence above the branch-cut 
singularity (beyond the first production threshold) instead of at origin of energies (q 2 = 0). This 
small modification also provides the opportunity to use Minkowskian data in our study instead of 
the Euclidean one. The relevance of this model-independent method to extract resonance poles 
is clear since does not depend on a particular lagrangian realization or modelization on how to 
extrapolate from the data on the real energy axis into the complex plane. 

The Montessus de Ballore's theorem states that the sequence of one-pole Pade Approximants 
Pf around x , 

N-l i , N 

P x (X,X ) = a k (x-X ) + : aW+T? ~ 7. (1-1) 

fe=0 [ 0) 

converges to F(x) in any compact subset of the disk excluding the pole x p , i.e, limjv^oo P 1 J (x, x ) = 
F(x). 

As an extra consequence of this theorem, one finds that the PA pole xpa = x + aN converges 
to x p for N — > oo. Since experiments provide us with values of Fj at different Xj instead of the 
derivatives of our function, we can use the rational functions P^ as fitting functions to the data 
in a similar way as in Ref. [1, 3] In this way, as N grows Pf gives us an estimation of the series 
of derivatives and the x p pole position. 

For a physical amplitude, the function F(x) (without a right-hand cut) is analytic from x = — oo 
up to the first production threshold x t h and within the disk B Xth (0). In the nn VFF case, the 
threshold is found to be at x t h — 4m^., where is the mass of the pion. One can use Pade 
Approximants in a safe way by using the Montessus' theorem and center the approximants at 
xq + i0 + over the brunch cut between the first and the second production threshold. In the tttt 
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Rcf. 


M p (MeV) 


r p (McV) 


[7] 


762.5 ±2 


142 ±7 


[8] 


754 ± 18 


148 ±20 


[9] 


763.0 ±0.2 


139.0 ±0.5 


[10] 


764.1 ± 2.7±2l5 


148.2 ± 1.9±|;^ 


this work 


763.7 ± 1.2 


144 ±3 



Table 1.1: Comparison of different results for the determination of the mass and width of the 
p-meson. 

VFF that would correspond to the range between xth = Am 2 and $th — 4to^ (assuming small 
multipion channels). In such a way we unfold the Second Ricmann sheet due to the analytical 
extension of the function F(x) from the first Riemann sheet at x + i0 + into the second one. 

In the case of resonant amplitudes, a single pole appears in the neighborhood of the real x 
axis in the second Ricmann sheet which is related to a hadronic state, a resonance. Once we 
have unfold the Second Riemann sheet by locating our approximants over the brunch cut, the 
application of the Montessus' theorem in the disk of convergence (which is the region defined 
between the thresholds) is straightforward and allows us to locate the position of the resonant 
pole if lies inside that region. If that is the case, our Pade approximant sequence determine 
systematically its position in a model independent way. 

1.2 Application of the method in a physical case 

We would like to use the method to analyze the final compilation of ALEPH tttt VFF data for 
the squared modulus \F 7T7T (q 2 )\ 2 , [5], and the / = J = 1 tttt scattering phase-shift 5 %v , identical 
to the 7T7T VFF phase-shift in the elastic region Am 2 < q 2 < Am 2 K (if multipion channels are 
neglected), which conforms the range of applicability of our P^ sequence (more details can be 
found in Ref. [6]). For N > 3 the fit x 2 already lies within the 68% confidence level (CL) and 
becomes statistically acceptable. At this point one needs to reach the typical fitting compromise. 
On one hand the experimental errors have an statistical origin and the contribution of this error 
increase as one considers higher order , with larger number of parameters. On the other hand, 
the systematic theoretical error decreases as N increases and the PA converges to the actual VFF 
and the phase-shift. For the p-mass determination we have taken N — 6 as our best estimate 
as the new parameters of PA with N > 7 turn out to be all compatible with zero, introducing 
no more information with respect to the previous Pf. Furthermore, the model studied in Rcf. [6] 
shows that the theoretical errors for mass and width results are smaller than 10~ 2 — 10 _3 MeV 
for N > 6, being negligible compared to the 0(1 MeV) experimental errors. This yields to the 
determination: 

M p = 763.7 ± 1.2MeV, T p = 144 ± 3MeV , (1.2) 

which is found to be in good agreement with previous determinations using more elaborated 
and complex procedures and with similar size of uncertainties as shown in Table 1.1. 

Similar results are found when analyzing the ir — ir phase-shift data for the tr-pole position, 
yclding in this case 

M a = 451 ± 21MeV, Y a — 619 ± 18MeV , (1.3) 
again in good agreement with previous determinations with more complex methods [4]. 

1.3 Conclusions 

We have develop a new model-independent method for extracting resonance poles (mass and 
width) from physical amplitudes based on the well defined mathematical theory of Pade Approx- 
imants which makes use of the Montessus' theorem to unfold de Second Riemann sheet of the 
desired amplitud to locate the resonant pole in the complex plane. This method systematize the 
algorithm of extracting the desired resonance pole (or low energy properties of these amplitudes as 
commented in the introduction) by building up a convergence sequence of those PA which at the 
very same time provides with a simple way to estimate the systematic errors. Our method has, 
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however, a larger application since does not rely on a particular lagrangian or in a modelization 
on how to extrapolate from the data on the real energy axis into the complex plane (such as other 
methods discussed in Ref. [11]). 

In the particular case presented here we have analyzed on one hand the experimental Euclidean 
data to obtain the first and the second derivatives at q 2 = for both the 7T-VFF and for the 
7T° TFF. And on the other hand, also Minkowskian data have been analyzed for 7T7T-VFF and 
also 7T7r-scattering data by a Pf* Pade Approximant sequence centered between the first and the 
second production thresholds (in such a way that we can unfold the second Riemann sheet of 
our amplitude). We have obtained a prediction of the p- meson and a- meson pole positions with 
compatible accuracy compared to other determinations but with a simpler and systematic method. 
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EXPLORING LOW-ENERGY QCD 

USING A GENERALIZED LINEAR SIGMA MODEL 

AMIR H. FARIBORZ* 

The decays of rf can provide very useful information on the physics of scalar mesons. For 
example, the decay rf — > -qmr involves intermediate scalar states /o(600) (or sigma), /o(980) and 
ao(980) as shown in Fig. 1.1. Therefore, at least in principle, such decays can shed light on 
various aspects of scalar meson. However, in practice, this objective is overshadowed by various 
uncertainties such as different final state interactions as well as various underlying mixing among 
scalar states. Nevertheless, one would hope to find useful approximations that could provide some 
insights into the physics of scalar mesons. The present work summarizes one such approach within 
the context of a generalized linear sigma model. 




Figure 1.1: Contribution of scalar mesons to decay rf — > i]irir. 

Below 1 GeV, the known scalar states are listed in Fig. 1.2 [1]: the light and broad isosinglet /o(600) 
or sigma, followed by isodoublet if* (800) or kappa meson and the two nearly degenerate states, 
isosinglet / (980) and isotriplet ao(980). Above 1 GeV, the known scalar states are also listed in 
Fig. 1.2: The isosinglet state / (1370) followed by the isodoublet #£(1430), the isotriplet a (1450) 
and the two isosinglet states /o(1500) and /o(1710). There have been numerous works addressing 
the physics of these illusive states and at times conflicting point of views have surfaced. Of course, 
getting into the history of these states and various approaches by different investigators is beyond 
the scope of this short summary and the reader is encouraged to refer to the available literature 
for details. However, in a nutshell, the scalar states below 1 GeV seem to be close to four-quark 
states [2] (with some deviations) and those above 1 GeV seem be close to two-quark states (with 
some deviations). To generate such deviations from pure four-quark and pure two-quark pictures, 
it seems natural to investigate a possible underlying mixing among four- and two-quark states. 
Such a mixing was first investigated within a nonlinear chiral Lagrangian framework in ref. [3] and 
it was shown that it provides a natural understanding of some of the puzzling aspects of if,* (1430) 
and ao(1450). The same type of mixing was further extended to the case of isosinglet scalar states 
below and above 1 GeV (which can also mix with scalar glueballs) in [4] , in which again consistent 
results were obtained and estimates on the quark and glue substructure of the isosinglet scalars 
were made. Such mixing patterns were further studied within the framework of a generalized 
linear sigma model which is the topic of the present summary. In this approach, the pseudoscalar 
mesons below and above 1 GeV (also given in Fig. 1.2) are brought into the picture. The main 
motivations for using a generalized linear sigma model are: First, it provides an independent check 
of the results obtained for scalar mixings in the nonlinear chiral Lagrangian approach mentioned 
above, and tests whether the results obtained are robust, or are artifacts of the framework used. 
Second, as scalar mesons are probed in various pseudoscalar interactions (such as, for example, 
7T7T, ttK, and ttt] scatterings and decays such as rf — > r/WTT, • • • ) it seems important to see what 
roles (if any) two- and four-quark mixings play in the physics of pseudoscalars mesons, for which, 
at least for pseudoscalars below 1 GeV, the phenomenology is fairly established, and as a result, 
allowing pseudoscalars to participate in these mixing patterns can provide a reliable probe into 
the understanding of the mixing. Third, some of the heavier pseudoscalar states above 1 GeV are 
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Figure 1.2: Scalars (left) and pseudoscalars (right) below 2 GeV together with their "bare" 
mixed) underlying two- and four-quark mixing nonets (middle). 



un- 



known to be dubious and speculated to be of a non quark- antiquark nature, therefore, extending 
the mixing mechanism to the pseudoscalar sector may provide useful insight into the nature of 
some of the heavier pseudoscalars. A generalized linear sigma model [which is formulated in terms 
of two scalar meson nonets (a two-quark nonet and a four-quark nonet) and two pseudoscalar 
meson nonets (a two-quark nonet and a four-quark nonet)] seems to be an appropriate framework 
for exploring the mixing patterns in both scalar as well as pseudoscalar sectors. The idea of 
mixing within this framework is investigated in Ref. [5] (and references therein) and the outcome 
(within the approximation used in [5]) is symbolically summarized in Fig. 1.2. In this approach, 
it is found that (as expected) the underlying four-quark scalar nonet is lighter than the two-quark 
scalar nonet, and the situation in reversed for the two pseudoscalar nonets (i.e. the two-quark 
pseudoscalar nonet tends to become lighter than the four-quark pseudoscalar nonet). This of 
course agrees with the conventional phenomenology for light pseudoscalars where their properties 
are understood based on a quark-antiquark picture. Further details can be seen in Ref. [5]. The 
notation is clarified in Fig. 1.2 in which the chiral nonet M is formed out of quark-antiquark scalar 
nonet S and quark-antiquark pseudoscalar nonet <f> (similarly, the chiral nonet M' is formed out 
of four-quark nonets S" and <p'). The general structure of the Lagrangian in terms of M and M' 
is: 

C = -Itt (9 M M^M f ) - ±Tr {d^M'd^M*) - V (M, M 1 ) - V SB , (1.1) 

where Vo(M, M') stands for a function made from SU(3)lX SU(3)r, (but not necessarily U(1)a) 
invariants formed out of M and M' and Vsb is the flavor symmetry breaking term. In dealing with 
this Lagrangian, two approaches have been considered: Approach 1 is based on the underlying 
chiral symmetry and the resulting generating equations [6] . In this approach no specific choice for 
the chiral invariant part of Vq is made, and only the axial anomaly, the Vs b and the condensates 
are modeled. The details of this approach can be found in Ref. [6] and are not discussed here. The 
main conclusions are: (i) The light pseudoscalars, as expected from conventional phenomenology, 
remain dominantly close to quark-antiquark states; (ii) The kappa meson tends to become a 
dominantly four-quark state; and (iii) This approach does not provide any information on ao 
and fo systems. To study ao and /o systems, specific choices for Vq should be made, and that 
leads to the second approach which is summarized here. In approach 2, a specific choice for the 
chiral invariant part of Vq is made [5] (in addition to modeling the axial anomaly, the Vsb and 
the condensates). The main conclusions are: (i) Again, light pseudoscalars have the tendency 
of becoming quark-antiquark states; (ii) Light scalars tend to become mainly four-quark states; 
and (iii) Predictions for various low-energy processes such as irir, irK, irr/ scatterings, and decays 
such as rj' decays, semileptonic decays of D s , ■ ■ ■ can be made. Here, we only give a summary of 
approach 2. Obviously, there are infinite number of terms that can be written down for Vq. Up 
to dimension 4, there are twenty one SU(3)l x SU(3)r invariant terms in Vq which can be made 
out of M and M' . To work with such a high number of terms, it seems reasonable to consider an 
approximation scheme in which Vq is organized in terms of the number of quark and antiquark 
lines. In the first approximation attempt: 1) We only consider terms with eight or fewer quark 
or antiquark lines; 2) We mock up the U(1)a anomaly exactly; 3) We ignore terms that are not 
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consistent with OZI rule; and 4) We consider the leading flavor symmetry breaking term that 
mocks up the quark mass term in the fundamental QCD Lagrangian. These approximations lead 
to the potential: 

V = - c 2 Tr(MM f ) + c\ Tv(MM^MM^) + d 2 Tv(M'M^) + e%(e abc e def M%M*M'f + h.c.) 



+ c 3 



, /detM\ /Tr(MM't) 
7lln ( v deTMtJ +(1 - 7l)ln U(M-Mt) 



+ k 1 [Tr(AM)+h.c.}. (1.2) 



where matrix A = diag. (Ai, A 2 , A 3 ) induces the flavor symmetry breaking and the two- and 
four-quark condensates of fields S and S' obey (S^) = 5 b a a a and (S' a ) — 5 h a j3 a . We assume isospin 
limit which implies: A\ = A 2 ^ A 3 and ot\ = a 2 ^ ot 3 and f3\ = (5 2 7^ /?3- At the level of the 
approximate potential of Eq. (1.2), altogether, there are 12 unknown parameters c 2 , c%, d 2 , eg, c 3 , 
71, ai, 0:3, pi, /?3, Ai, A3 that need to be determined using 12 inputs. We take the following eight 
experimental inputs: m[a (980)] = 984.7 ± 1.2 MeV; m[a (1450)] = 1474± 19MeV; to[tt(1300)] = 
1300 ± 100 MeV; = 137 MeV; F n = 131 MeV; = 20 ^ 30; det(M^) = det(Af2) cx P- ; and 

Tr(M^) = Tr(M2) 0X P, together with four minimum conditions: = (||^) = (^r) = 

(t^ 5- ) = ^' These inputs and minimum conditions allow a determination of the Lagrangian 
parameters and the detailed numerical analysis of this study in given in [5] . The main uncertainties 
are clearly on m[7r(1300)] as well as the ratio and therefore, the results have some variations 
that reflect these two uncertainties. The work o [5] shows that the underlying mixing among two- 
and four-quark pseudoscalars does not change the well-established picture of the light pseudoscalars 
as dominantly being quark-antiquark states (the situation is of course reversed for the heavy 
pseudoscalars). The present level of this investigation predicts four etas: Two of the predicted 
etas are below 1 GeV and are close to the physical states 7?(547) and ?i(958). The two heavier etas 
above 1 GeV may be identified with two of the several physical eta states above 1 GeV. In the 
work of [5] various scenarios for this identification is studied in detail and the closest agreement 
corresponds to identifying the two predicted etas above 1 GeV with ri(1295) and 77(1760) . The 
situation for scalar mesons is the opposite of pseudoscalars: Those below 1 GeV come out close to 
four-quark states and those above 1 GeV tend to become more of two-quark states. These natural 
tendencies for scalar mesons are consistent with those studied within nonlinear chiral Lagrangians 
described above. Therefore, in summary, it seems that the generalized linear sigma model discussed 
here is an appropriate framework for learning about scalar and pseudoscalar mesons, and already 
has provided useful insights into their quark substructure. Further extensions will involve inclusion 
of higher order terms in the potential [with higher number (more than eight) of underlying quark 
and antiquark lines]. Also inclusion of glueballs seems to be quite important. Various predictions 
of the model for several low-energy processes are being currently looked at and will be reported 
elsewhere. Particularly, the present approach is hoped to provide useful insights into the physics 
of 7/ and the eta's above 1 GeV. 
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1.1 Introduction 

The r\ meson provides an excellent laboratory for the study of matter in the low-energy regime. 
Because all strong and electromagnetic decays of the r\ are forbidden in the first order, rare 
processes are accessible which would otherwise be highly suppressed. The measurement of the 
decays of the ij can test fundamental symmetries, probe the structure and dynamics of hadrons, 
and investigate possible new physics. 

A major goal of the WASA-at-COSY experiment is to make precision measurements of r\ 
decays. r\ mesons are produced by impinging a proton beam onto a proton or deuteron target. 
With the WASA detector, both forward-scattered hadrons and wide-angle decay products of the 
r] can be measured. 

WASA has recorded 3 x 10 7 r\ mesons in the pd — > 3 Her/ reaction. In the pp —>■ ppr\ reaction, over 
10 8 r\ mesons have been produced with a trigger biased towards charged decays. Proton-deuteron 
reactions feature a favorable ratio between r\ production and background processes, but a lower 
rate of r\ production of about 10r;s _1 . Proton-proton reactions have a factor of ten higher rate of 77 
production but significantly higher background. The two production methods are complementary. 

1.2 Current Work on Selected Decay Channels 

In the following, the current status of the analysis of selected decay channels will be described. 

The decay r\ — > tt + it~j provides the possibility to search for a contribution of the box-anomaly 
term from the Wess-Zumino-Witten (WZW) Lagrangian [1, 2]. The predicted decay width cal- 
culated with chiral perturbation theory is significantly lower than the measured value, partially 
because the phase space available for this decay is constrained by the mass of the tt and is far 
from the chiral limit. Attempts have been made to include final-state interactions as unitarized 
extensions to the WZW Lagrangian. In order to test these predictions, both the decay rate and 
the dynamics of the final state particles can be investigated. 

In the WASA measurement, a total of 13960 ± 140 event candidates for this reaction were 
recontructed from pd — > 3 Herj data and the shape of the relevant kinematic distributions have 
been determined [3] . The data deviates from calculations using the simplest gauge-invariant matrix 
element. For a more realistic description, a form factor as described in [4] has been applied to 
the prediction of the simplest matrix element. This allows for a model independent description of 
the distribution with a single free parameter, a. This parameter can be compared to theoretical 
predictions. 

Another aspect of this decay is the measurement of the branching ratio. A value about 3-cr 
lower than the former experimental average was recently measured by the CLEO collaboration and 
confirmed by the KLOE collaboration [5, 6]. As a next step for this channel, WASA will measure 
both the branching ratio and the photon-energy dependence in one analysis from high-statistics 
proton-proton data. 

The decay r\ — > tt + tt~ e + e~ is a rare decay with a branching ratio of (2.68±0.11) x 10~ 4 [7]. One 
goal of the analysis is to confirm the experimental value of the branching ratio, which has only been 
measured once with high statistics [8]. This decay is interesting because it may provide a novel 
test for CP-violation outside of the standard model. It has been shown that the asymmetry of the 
angle between the electron and tt decay planes is sensitive to a contribution of the CP-violating 
E-l transition to the decay amplitude [9, 10]. Unlike the reaction 77 — > 7r + 7r~, CP- violation in this 
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system is not suppressed by experimental limit on the electric dipole moment of the neutron. The 
theoretical upper limit for this asymmetry is about 1%. 

So far (230 ± 22) signal event candidates have been identified in proton-deuteron data, shown 
in Figure 1.1a. This number is consistent with a branching ratio near the measured and predicted 
values. Systematic studies are currently being performed in order to extract the branching ratio 
with the lowest possible error. The decay is also being analyzed in proton-proton data where 
a larger sample of event candidates is visible. The goal is to measure the branching ratio with 
proton-deuteron data and to subsequently use the proton-proton data to confirm the branching 
ratio and measure the decay plane asymmetry. 

Because it proceeds via two virtual photons, the decay r\ — > e + e~e + e~ has a very small 
branching ratio, with a QED prediction of 2.6 x 10~ 5 . It has been measured only once and the 
experimental value was found to be within errors of the QED calculation [11]. 

WASA has measured this decay in the pd — > 3 Her] reaction. Two independent analyses report 
(50 ± 14) event candidates. The event sample from one analysis is shown in Figure 1.1b. The 
branching ratio returned by these two analyses is consistent, allowing WASA to quote a preliminary 
value of (3. Q±0. 8 s tat±0 -7 sys . (norm)) x 10 -5 which is in agreement with theoretical and experimental 
values. Systematic studies as well as consistency checks between the two analyses are currently 
being performed. 



Data 
MC Total 

MC n -* n + n _ [Tr°-> e*e"y] 
MC r\ —> ji + jt"e + e" 
MC pd^ 3 Heft*;i 7t° 




0.56 0.58 0.6 0.62 
Missing Mass 3 He [GeV. c 2 ] 



(a) 



i0 ^_ — • — data 

;o l_ sum i) channels (w/ signal) 

r,-»e*e-e*e- 

sumr|-background channels 
10 ^- (w/o signal) 

;0 ~ preliminary 



0.55 0.6 
'He-missing mass [GeV/c z ] 



(b) 



Figure 1.1: Missing Mass of 3 He for events satisfying conditions for (a) rj —> tt + tt e + e and (b) 



rj — > e + e e + e 



The decay of rj into two electrons is a fourth order electromagnetic process proceeding through 
two virtual photons and additionally suppressed due to helicity conservation to the order of 10~ 9 . 
This extremely low branching ratio makes the reaction very sensitive to any processes that could 
enhance the decay width, including possible effects outside the standard model [12, 13, 14]. Mea- 
surements of the closely related 7r° — > e + e~ decay have found branching ratios above theoretical 
predictions. Measurements in the rj system should be even more sensitive to such effects due to 
the lower branching ratio predicted by the standard model. 

The current experimental upper limit for this decay is 7.7 x 10~ 5 [7] and was measured by 
the CELSIUS- WASA collaboration. Using two weeks of data in proton-proton interactions from 
WASA-at-COSY, a value one order of magnitude lower has been calculated. When considering 
the full statistics available in the pp — > pprj reaction, it is expected that additional sensitivity can 
be gained corresponding to at least another order of magnitude improvement. Additionally, a 
more intelligent trigger has been developed and is being tuned to enable data collection at higher 
luminosities. 



1.3 Conclusion 

The WASA-at-COSY collaboration is in the processes of analyzing the large amount of data taken 
on r\ decays while simultaneously gearing up for extended production runs to increase statistics 
for rare processes. Results for some decay channels have been released while several more are at 
the stage of careful systematic checks. 

Additionally, the meson program at WASA-at-COSY is expanding. A pilot run intended to 
produce large numbers of tt° mesons in order to measure the very rare ir° — » e + e~ decay is currently 
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being analyzed. Also, the first data on w decays was taken in proton-deuteron and proton-proton 
reactions in Spring of 2011. 
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1.1 Introduction 

At low energies the coupling constant of Quantum Chromodynamics (QCD) becomes large and 
renders perturbation theory useless. The phenomenology of the strong interactions in this domain 
remains one of the major challenges of modern particle physics. Experimental study of decays of 
light mesons in general and ir° , r\ and rj pseudoscalars in particular is a unique way to explore low 
energy QCD. By comparing data obtained at different facilities like KLOE, CLEO, BES, MAMI, 
COSY, and CLAS with a predictions of effective field theory at low energy, Chiral Perturbation 
Theory (ChPT), one can get an insight into the non-perturbative regime of strong interactions and 
provide important information for a firmer foundation of hadronic physics rooted in the standard 
model. Until recently the quest for high statistics accurate experimental data was not fulfilled. 

In this talk we give a brief overview of experimental data on a decay of light mesons collected 
with the CLAS setup at JLAB either with unprecedented high statistics or comparable with the 
worlds highest statistics measurements depending on the decay mode. Below we outline the quality 
and statistics of available data at CLAS in a different decay channels of light mesons, rather than 
present results of finalized analyses that are in progress. 

1.2 Light mesons in CLAS 

The tagged beam of photons produced from 6 GeV CEBAF electron accelerator is used to study 
photo-production of light mesons from liquid hydrogen target. Until now the main focus of the 
experiments in CLAS related to light mesons was the study of different aspects of photo-production 
cross sections aiming to understand mechanism of their production. However as it is shown below 
the wealth of data allows to make a studies exploiting different decay modes detached from the 
production vertex and considering decaying mesons as a laboratory by themselves. 

Experimental data are presented by accenting on photo-production reactions j+p — > p^°(r], v') 
collected in the following decay modes: 

• Dalitz decays ir°(ri,r]') — > e + e~-f 

• Radiative decays 77(77') K + n~-f 

• Hadronic decays 17(17') — > 7r + 7r~7r° 

• Hadronic decay 77' — > tt + tt^t] 

1.2.1 Dalitz decays 

The branching ratios of radiative decay of pseudoscalar mesons ir° and 77 were measured and are 
recorded in PDG, however there is only an upper limit quoted for r{ . Below in Fig. 1.1 (left panel) 
we present distribution of the invariant mass of e + e~j measured in photoproduction reaction with 
CLAS. As one can there are clear peaks of all three light pseudoscalar mesons ir°, 77 and eta'. One 
of the less studied channel of the 77 — > e + e~-f decay allows to extract form factor of 77 in the 
timelike region by studying distribution of the invariant mass of e+e~ pairs presented in Fig. 1.1 
(right panel). In both panels in Fig. 1.1 number of events are presented without corrections due 
to the CLAS acceptance and efficiency. 
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P(i" '1 'l') I Invariant mass of e + e for e + ey 




M(eVy) [GeV] M(eV) [GeV] 



Figure 1.1: (Color online). Left panel: the e + e~7 invariant mass distribution in the reaction 
7 + p — > pe + e~7, the vertical lines show positions of pseudoscalar mesons 7T , r], r)'. Right panel: 
invariant mass of e + e~ lepton pairs for events under the 77 peak M(e + e~j) = 0.548 ± 0.01 GeV. 



1.2.2 Radiative decays 

Another very interesting decay of rj and rj' mesons is their radiative decay 77(77') — > tt + it~j. Recent 
theoretical paper [1] demonstrates that the photon energy distribution from radiative decay of 77 
and 77' mesons provides very important test of theoretical models. In Fig. 1.2 (left panel) missing 
mass of the proton is presented for exclusive reaction 7 + p —¥ pir + Tr~j. In both 77 and 77' peaks 
a statistics is by more than an order of magnitude higher than existing world data (see [1] and 
references therein). In Fig. 1.2 (middle (right) panel) distribution of the energy of the photon 
from the decay of 77(77' in the center-of-mass frame of the parent is are presented. 
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Figure 1.2: Left panel: distribution of missing mass of the proton in exclusive reaction 7 + p — > 
p7r + 7r _ 7. Middle panel: the distribution of the 7 energy in the center-of-mass of 77. Right panel: 
same as the middle panel plotted for the 77' decay. All three plots show raw number of events 
before the acceptance and efficiency corrections. 



1.2.3 Hadronic decays 

In this section we present experimental data for the reaction 7 + p — > pir + ir~ X(ir°(r])), where 
7T° (77) arc identified via missing mass of proton, ir + and ir~ . In Fig. 1.3 (left panel) a distribution 
of missing mass of the proton is presented showing clear peaks of 77 and ui mesons with ~2M and 
~20M events respectively as well as hints of 77' and <f> mesons. To see peaks of 77' and 4> better 
in Fig. 1.3 (right panel) we plot the same distribution for the range above lo meson. As one can 
see there we observe one of rare decays 77' — > ir + ir~Tr° (Br=3.6 x 10 -3 ) and OZI violating decay 
(f> —> 7r + 7T~7r° (Br=15.3%) observed for the first time in the photo-production experiments. The 
photo-production of the latter can be studied to check the influence of meson-baryon interference 
that might occur in the charge and neutral KK decay modes and is absent in this channel. 

In Fig. 1.4 (left panel) we present distribution of missing mass of the proton in the reaction 
7 +p — >■ p7r + 7r~77, where 77 is reconstructed in the missing mass of proton, tt + and tt~ . As one can 
see there is a clear peak of 77' with ~340K events which is almost an order of magnitude higher 
than recent BES [2] data. This will allow to extract coefficients of polynomial distributions of 
Dalitz plot projections with much higher statistical precision. On the other hand there appears a 
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Figure 1.3: Left panel: distribution of missing mass of the proton in the reaction 7+p — > pir + ir tt°. 
Right panel: the same for the range of invariant mass above u> meson production. 

hint of G-parity decay mode of <fi — >• -K + n~r). In In Fig. 1.4 (right panel) the same distribution is 
plotted above the rj' peak. Here we see clear peak of cj) meson. This decay mode has never been 
observed and PDG quotes only an upper limit Br< 1.8 x 10~ 5 . Theoretical model [3] predicts 
even smaller upper limit on the order of ~ 3 x 10~ 7 . 




Figure 1.4: Left panel: distribution of missing mass of the proton for the reaction 7+p — > p-K + n r\. 
Right panel: the same for the range of invariant mass above 7/ meson production. 
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1.1 Introduction 

In the low energy regime of Quantum Chromodynamics, the interaction between quarks and gluons 
cannot be treated perturbatively and so far the understanding of the processes governed by the 
strong forces is unsatisfactory. Therefore, it is essential to carry out measurements involving the 
production and decay of hadrons and to interpret them in the framework of effective field theories 
experiencing recently an enormous development in applications to the description of meson decays 
and production. In this contribution we concentrate on the r\ meson. The progress in understand- 
ing of the production processes of the r\ meson will strongly rely on the precise determination of 
spin and isospin observables. So far these observables have been determined only for few excess 
energies and with low statistics [1-5]. 

1.2 Partial waves 

For an unambiguous understanding of the production process relative magnitudes from the partial 
waves contributions must be well established. This may be achieved by the measurement of the 
analysing power which would enable to perform the partial wave decomposition with an accuracy 
by far better than resulting from the measurements of the distributions of the spin averaged cross 
sections [6, 7]. Analysing power A y may be understood as a measure of the relative deviation 
between the differential cross section for the experiment with and without polarized beam (a and 
0o respectively): 

a(0, ip) = a (0, ip) ■ [1 + A v {6) ■ P • cos(p)] (1-1) 
where P denotes the beam polarization. 

1.3 Studies of A y with the WAS A detector at COSY 

Using the WASA-at-COSY detector [8] we intend to determine the energy and the angular de- 
pendence of A y (Q,6) and the total and differential cross sections for the pp — > ppr\ reaction in the 
excess energy range from the threshold up to 100 MeV. In November 2010 first measurements for 
Q = 15 MeV and Q = 72 MeV have been conducted [9]. From Table 1.1 we can see the beam 
parameters and the expected number of events for each excess energy. 



Q [MeV/c] 


P [MeV/c] 


<7tot\pfo] 


Acceptance 






15 


2026 


10 H 


0.55 


99770 


81861 


72 


2188 


5 • 10 a 


0.63 


447739 


375580 



Table 1.1: Estimate of the number of produced rj mesons 



Protons from the pp —> pprj reaction are registered in the forward part of the detector and 
photons from the r] meson decay are detected in the Electromagnetic Calorimeter of the central 
part. Simultaneously to the pp — > pprj reaction elastically scattered protons were registered. The 
pp — > pp reaction will be used for monitoring of the polarization degree, luminosity and the detector 
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performance. In the case of the pp —> pp reaction one proton is registered in the Forward Detector 
and the other in the Central Detector. 

Within the geometrical acceptance of the Forward Detector (from 3° to 18°), the Central 
Detector covers proton scattering angle from 60° to 84°. In the center of mass system that 
corresponds to a scattering angle in the range of 30° to 46° as seen in Fig. 1.1a. 

In Fig. 1.1b the scattering angle of one proton is depicted as a function of the scattering angle 
of the second proton. The first proton is registered in the Central Detector, the second one in 
the Forward Detector. Colored lines are predictions for different beam momenta. The blue line 
corresponds to a beam momentum of 2 GeV/c, close to the experimental value. 
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(a) The scattering angle in the centre of mass 
system is presented as a function of the scatter- 
ing angle in the laboratory frame. 
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(b) The distribution of the scattering angle of 
proton registered in the Central Detector as a 
function of the scattering angle of the proton 
registered in the Forward Detector. For more 
details see the text. 



Figure 1.1: The protons angular distributions in the Monte Carlo simulations of the pp — > pp 
reaction. 



1.4 Extraction of A y from the experiment 

In the first step we can obtain the value of the polarization, P, from the pp — > pp reaction using 
the analysing power already measured by the EDDA experiment [10]. In the second step, we will 
calculate the geometrical averages, N±, of the number of rj meson produced in direction (6, ip) 
during the spin down and spin up modes, defined as follows: 



A_ 


= A 


JnIn^ = a 


/ K 


N L 
6 L L± 


N+ 




= A, 


1 K 


jvi 

CrL^ 



(1.2) 



Ni and Nl can be determined according to the Madison conventional]. Knowing the angles 6 
and <p of the outgoing r\ meson we can calculate the analysing power for the three particle final 
state using the following formula: 



A (0) N +(0^)-N-(6,y) 
y[ ' ~ Pcostp ' N + (9,<p)+N_(9,<p) [ ' 

where the polarization P was calculated in the first step. 

An A y for the proton-proton elastic scattering at 2026 MeV and 2188 MeV beam momenta, 
based on the results of the EDDA experiment is shown in Fig. 1.2. We can conclude that values 
of A y for Q=15 MeV and Q=72 MeV are within the range of [0.32-0.38]. 
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Figure 1.2: Distribution of A y for pp — > pp reaction as a function of the protons scattering angle, 
9, in the centre of mass system. Square points denotes data based on EDDA experiment [12] and 
superimposed line indicates results of Rcf. [13]. 
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HIGH PRECISION rj MESON MASS 
DETERMINATION AT COSY-ANKE 



PAUL GOSLAWSKI* for the AN KE- Collaboration 



Recent measurements on the rj meson mass performed at different experimental facilities (i.e. 
CERN-NA48, COSY-GEM, CESR-CLEO, DA$NE-KLOE, MAMI-Crystall Ball) resulted in very 
precise data but differ by up to more than eight standard deviations, i.e. 0.5 MeV/c 2 [1]. In 
order to clarify this situation a high precision measurement using the ANKE spectrometer at the 
COoler SYnchrotron COSY has been realized. 

Using the two-body reaction dp — > 3 Her? at low excess energies the rj mass can be determined 
only from pure kinematics by the determination of the production threshold. Therefore, twelve 
data points at fixed excess energies in the range of Q = 1 — 12 MeV were investigated. The final 
state momentum pf of the 3 He -particles in the Center of Mass System (CMS) 



Pf (s) = 



(toshc + rn v Y 



(rasHe - m v Y 



(1.1) 



measured with the ANKE spectrometer, is very sensitive on the rj mass and the total energy y/s, 
where the latter one is completely defined in a fixed target experiment by the masses of the initial 
particles and the momentum of the deuteron beam pd- 



2m, 



pi 



(1.2) 



For a precise determination of the production threshold both quantities, the final state momenta 
of the 3 He -particles and the corresponding deuteron beam momenta have to be measured with 
highest accuracy. Fitting the dependency of the final state momentum pf on the beam momentum 
Pd and the rj mass, pf = Pf{pd, rn v ), for the twelve points the mass can be extracted as fitting 
parameter (see Fig. 1.3). 



The beam momentum for each fixed excess energy was determined using a method developed at 
VEPP-2M at Novosibirsk at an electron-positron machine [2] using the spin dynamics of a polarized 
beam. Here the spin precession frequency of a relativistic particle is disturbed by an artificial spin 
resonance induced by a horizontal rf-magnetic field of a solenoid leading to a depolarization of the 
polarized accelerator beam. The depolarizi ng reson ance frequency f r depends on the kinematical 
7-factor (i.e. the beam momentum p = myj 7 2 — 1) and the beam revolution frequency /o via the 
resonance condition: 

/ r = (fc + 7G)/o, (1-3) 

where k is an integer and G the gyromagnetic anomaly of the beam particle. By measuring these 
two frequencies the beam momentum of a polarized beam can be determined with a precision 
below Ap/p < 10~ 4 . For the first time this method was used at COSY with a vector polarized 
deuteron beam and the momenta in the threshold range of 3.1 — 3.2 GeV/c were determined with 
an accuracy of Ap/p < 6 • 1CP 5 , i.e. with 170 keV/c [3]. 

The correct final state momenta for the twelve different energies of the 3 He -nuclei of the reaction 
dp — > 3 Her7 can only be extracted fulfilling two conditions: a clear identification of the reaction of 
interest and a precise detector calibration. At ANKE the produced 3 He -nuclei can be identified 
by the energy loss and time of flight information. By this, the background consisting mainly of 
protons and dcutcrons of the dp elastic scattering and the deuteron break-up, can be suppressed 
effectively. For a two-body reaction at a fixed center of mass energy t/s the final state momenta in 
the CMS are distributed on a momentum sphere with constant radius pf, which can be visualized 
by plotting the transversal versus the longitudinal reconstructed momentum, as shown in Fig 1.1a. 

*Institut fur Kcrnphysik, Westfalischc-Wilhelms Univcrsitat Miinstcr, 48149 Miinstcr, Germany, 
paul . goslawskiOuni-muenster . de 
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Figure 1.1: a) The momentum loci for the 3 He?7 and 3 He7r channels. For the 3 He?7 channel 
ANKE covers near threshold the full solid angle, while for the 3 He ir channel only forward scattered 
3 He -nuclei are detected, b) Final state momentum pt = p/(cos??) (red) at an excess energy of 
Q ps 1 MeV, the background description (blue) and the extracted 3 FIc?7 signal (green). 



According to Eq. 1.1, one expects a centered momentum locus with a fixed radius Pf = (p 2 . +Py + 
pi) 1 / 2 , indicated in Fig. 1.1a as dashed line. Using the feature that the ANKE facility has full 
geometrical acceptance for the reaction dp — > 3 Hery near threshold up to 15 MeV the detector 
calibration can be verified and improved by studying the kinematics of this two-body reaction. The 
main idea to verify the calibration is that the momentum sphere has to be completely symmetric 
in p x , p y and p z (or ?? and <fi) so that the final state momentum pt should be constant in all 
directions. By a careful investigation of the momentum dependency on the cosine of the polar 
angle ?? and the azimuthal angle <f> 

Pf = pf (cos??) (1.4) 
Pf = Pf{4>) (1-5) 

the shape of the momentum sphere or locus can be studied in more detail. Deviations from this 
symmetric shape will indicate the need of an improvements of the calibration. Therefore the 3 Hc ?; 
signal has to be extracted background free. The background left after cutting on 3 He -nuclei (see 
Fig. 1.1a originating mainly from the multi pion production can be subtracted by data taken below 
the rj threshold at an excess energy of Q « —5 MeV, but analyzed as if they were taken above. 
In [4] the successful applicability of this approach on missing mass spectra is shown, but it is also 
applicable to final state momentum spectra as shown in Fig. 1.1b for different cos?? bins for a data 
point close to threshold, i.e. Q sa 1 MeV. Similarly the final state momentum dependency on </> 
can be studied. The background free 3 He?7 distribution allows to extract the mean final state 
momentum pj(cos??) and Pf{<t>) as shown in Fig. 1.2. In contrast to Monte-Carlo simulations on 
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Figure 1.2: Final state momentum dependency: pf — p/(cos??) and Pf{4>) for data (red) simula- 
tions without momentum smearing (black line) and with momentum smearing (black). 

p/(cosi?) and Pf(4>) without momentum smearing (black line) the obtained data points indicated 
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as red crosses show a dependency of the final state momentum in cost? and (f>. The shape of the 
momentum sphere is stretched to values cos-d — >• ±1 and shows an oscillation in tp. This behavior is 
caused by a kinematic effect due to different momentum resolutions of the ANKE detector for p x , 
p y and p z . Assuming that the p x , p y and p z distributions are gaussian distributed with different 
widths a, it is possible to reproduce the final state momentum dependency on cost? and <f> with 
Monte-Carlo simulations shown as black crosses. 

For the final state momentum determination the plots shown in Fig. 1.2 are of high importance 
because of following three reasons 

i) Improvement of the calibration. Asymmetric shapes in cos # and (j) can be corrected by 
minor changes of the ANKE calibration parameters. These changes are so small that typical 
calibration quantities like missing masses of different reactions show no variations, because 
they are not sensitive enough. 

ii) Extraction of the correct momentum resolution in (p x ,p y ,p z ). Assuming gaussian distribu- 
tions with different widths we found for ANKE (a x , a y , a z ) = (3.2, 7.8, 16.4) MeV/c. 

iii) Correction of the reconstructed final state momentum. Because of the finite momentum 
resolution the extracted momentum (red) are larger compared to the original one (black 
line) . This has to be considered in the determination of the twelve final state momenta. 

Currently the final state momentum analysis is still in progress but already now the twelve mo- 
menta in the range of 30 — 100 MeV/c can be determined with a precision below 400 keV/c. Fitting 
the pf — pj(pd, rn v ) dependency, as shown in Fig. 1.3 the rj mass is determined preliminary to 

m v = (547.869 ± 0.007 ± 0.040) MeV/c 2 . (1.6) 

The accuracy, which will be achieved at ANKE, will be comparable and competitive to the precision 
achieved at other recent experiments. 




Figure 1.3: Determination of the r\ mass by fitting the dependency pf — Pf(pd, Tn v ). 
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1.1 COSY-11 detection set 

The reported experiment has been performed 
in the Research Centre Jiilich at the cooler syn- 
chrotron COSY [1] by means of the COSY-11 
detector system [2, 3] presented in Figure 1.1. 

The collision of a proton from the COSY 
beam with a proton cluster target may cause 
an rj meson creation. In that case all outgoing 
nucleons have been registered by the COSY-11 
detectors, whereas for the rf meson identifica- 
tion the missing mass technique was applied. 

The cluster target width was reduced by 
a factor of 5 compared to the normal opera- 
tion mode in order to decrease the systematic 
uncertanities during the determination of the 
total width of the 7/ meson. The dimensions 
of the target stream were determined using a 
wire probe installed for this purpose. 




Figure 1.1: COSY-11 detector setup. From left 
to right: quadrupole (yellow) and dipole (orange) 
magnets of COSY, two drift chambers (silver) 
and two scintillator detectors (blue-black). 



1.2 Total width of the 7/ meson 

The mass of unregistered meson was determined via the missing mass technique, while the total 
width was directly derived from the mass distributions established at five different beam mo- 
menta [4, 5]. Based on the sample of more than 2300 reconstructed pp —> ppr{ events the deter- 
mined total width of the rf meson amounts to T rf = 0.226 ± 0.017(stat.) ± 0.014(syst.) MeV/c 2 , 
which is the most precise measurement until now [6, 7, 8]. On the plots in Figure 1.2 the exper- 
imental data are presented as black points, while the histograms correspond to the sum of the 
Monte Carlo generated signal for T v > = 0.226 McV/c 2 and a normalised second order polyno- 
mial obtained as a fit to the signal- free background region for another energy [4, 5]. The plot in 
the bottom right corner presents the x 2 distribution as a function of the r^' determined by the 
comparison of the missing mass spectra from the measurement with Monte Carlo simulations. 
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Figure 1.2: The missing mass spectra for the pp 
system equal to 4.8, 2.8, 1.7, 1.4, and 0.8 McV (from left to right, top to bottom). The rf meson 
signal is clearly visible. The experimental data are presented as black points, while in each plot the 
histogram corresponds to the sum of the Monte Carlo generated signal for = 0.226 MeV and 
the shifted and normalised second order polynomial obtained as a fit to the signal-free background 
region for another energy. Bottom right: \ 2 as a function of T v i in the range up to \ 2 = Xmin + 1- 
The minimum corresponds to the determined T v > value. 
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1.1 Introduction 

More than fifty years ago Abashian, Booth and Crowe [1] observed in the reaction pd -^ 3 HcX 
that the missing mass of the detected 3 He ejectile corresponding kinematically to the production 
of a pion pair, i.e. X = -kit, contains a peculiar enhancement starting right at threshold at low 
7T7r-masses. In subsequent bubble-chamber [2, 3] and single-arm magnetic spectrometer measure- 
ments [4, 5, 6, 7, 8, 9, 10, 11] this enhancement was found in fusion reactions to d, 3 He and 4 He, 
where an isoscalar pion pair was produced, however, not in the fusion reaction leading to tritium, 
where an isovector pion pair was produced. From these results it was concluded that this effect, 
which later-on was named ABC effect after the initials of the first authors, only appears in the 
double-pionic fusion in combination with the production of an isoscalar pion pair. 

Since these experiments were undertaken at beam energies, which energetically allow the mu- 
tual excitation of two colliding nucleons into their first excited state, the A(1232)P 33 , theoretical 
attempts were undertaken to understand the ABC effect by the excitation of a AA system via 
i-channel meson exchange or variations of it [12, 13, 14, 15, 16, 17]. These calculations, indeed, 
predicted a low-mass enhancement, however, also a high-mass enhancement. Surprisingly enough, 
the inclusive single-arm measurements seemed to support such a high-mass enhancement. And it 
was only recently that the first exclusive and kinematically complete measurements of solid statis- 
tics [18, 19, 20, 21, 22, 23] over practically the full phase space revealed that such a high-mass 
enhancement is not present. 



1.2 Recent Experiments 

In order to understand the two-pion production process in pp, pd and dd collisions in more detail, 
systematic experimental studies were initiated at CELSIUS and continued later-on at COSY. These 
exclusive and kinematically complete high-statistics measurements have been carried out at the 
detector installations PROMICE/WASA [24, 25, 26], CELSIUS/WASA [27, 28, 19, 29, 30, 31, 32], 
COSY-TOF [33, 34], WASA-at-COSY [35, 23, 36] and COSY-ANKE [37] covering the energy 
range from the two-pion production threshold up to yfs = 2.5 GcV (corresponding to a nuclcon 
beam energy of 1.4 GeV). The bulk of these experiments has been performed with the WAS A 
detector setup, which has an angular coverage of nearly 47r and features windowless hydrogen 
and deuterium pellet target systems [27, 35] as well as a very reliable particle identification for 
7, ir,p, n, d, 3 He and 4 He. Measurements at WASA generally are performed with detection of all 
ejectiles of an event (with the exception of spectator nucleons) leading thus to kinematic fits with 
overconstraints in the data analysis. For the investigation of the pn initiated reactions the quasifree 
pd process with a spectator proton resulting from the target deuteron is utilized. By exploiting 
the Fermi motion of the active nucleon in the deuteron target that way the energy dependence of 
pn and pp initiated reactions can be measured over a range of more than 100 MeV by use of a 
single beam energy. 

*Physikalisches Institut der Universitat Tubingen, Germany 
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[GeV] 

Figure 1.1: Total cross sections obtained from this experiment on pd — > dir°ir° + p spe ctator for 
the beam energies T p = 1.0 GeV (triangles), 1.2 GeV (dots) and 1.4 GeV (squares) normalized 
independently. Shown are the total cross section data after acceptance, efficiency and Fermi mo- 
tion corrections. The hatched area indicates systematic uncertainties. The drawn lines represent 
the expected cross sections for the Roper excitation process (dotted) and the i-channel AA con- 
tribution (dashed) as well as a calculation for a s-channel resonance with m = 2.37 GeV and 
r = 68 MeV (solid). From Ref. [23] 

1.3 Results 

From the overwhelming amount of high-quality data for total (integral) and differential cross 
sections in the various two-pion production channels the following features emerge for the energy 
region from threshold up to y/s — 2.6 GeV: 

• The isovector two-pion production [38] initiated by pp collisions can be quantitatively 
understood by t-channel meson exchange with 

— excitation of one of the colliding nucleons into the Roper resonance, N*(144Q)Pn and 
its subsequent decay into the Nttt: system either directly {Na channel) or via the A7r 
intermediate state, 

— mutual excitation of the two colliding nucleons into the A state each with subsequent 
decay into Nn, or 

— excitation of a higher lying A state, favorably the A(1600)P33 state with subsequent 
decay into Nrnr via An. Such an excitation, however, plays a dominant role only in 
the nnn + Tt + channel possessing an isotensor pion pair. 

Whereas the Roper process dominates at lower energies close to threshold, the AA process 
dominates at higher energies close to y/s = 2mA = 2.46 GeV. In particular the isovector 
double-pionic fusion reaction pp — > rf7r + 7r° is well described in its energy dependence up to 
at least 3 GeV by the latter process [30]. 

• The isoscalar two-pion production reaction pn — > dn^ir - measured as quasifree process 
in pd collisions [19, 23] - exhibits a pronounced low-mass enhancement in the 7T7r-invariant 
mass distribution as expected in view of the ABC history on double-pionic fusion reactions. 
However, what was not expected is a strict correlation of this ABC effect with the appearance 
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of a resonance-like structure in the total cross section - see Fig. 1 - with m = 2.37 GeV and 
r = 70 MeV, i.e. roughly 100 McV below 2mA and three times narrower than expected from 
the t-channel AA process. From the angular distributions the spin-parity of this structure 
is derived with the result that I(J P ) = 0(3 + ). 

Note that this reaction constitutes the only purely isoscalar two-pion production channel in 
NN collisions. Hence exotic processes in the isoscalar channel are likely to be seen best in 
this particular reaction without being diluted by isovector contributions. 

• The isoscalar double-pionic fusion reaction dd — > 4 He7T7r also exhibits the ABC effect in strict 
correlation with a resonance- like structure in the total cross section at the same energy above 
the two-pion threshold [21, 22]. The increased width of this structure can be explained by 
the Fermi motion of the active isoscalar pn pair in target, projectile and fused nucleus. 
The energy dependence of the reaction pd — > 3 He7T7r is under investigation. Preliminary 
results [39] indicate a similar trend as observed for the 4 He. If true than this isoscalar pn 
resonance structure is obviously strong enough to survive even in helium nuclei. 
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1.1 Introduction 

The ABC effect named after Abashian, Booth and Crowe [1] denotes a low-mass enhancement 
in the 7T7r-invariant mass distribution resulting from two-pion production processes. More than 
fifty years ago Abashian, Booth and Crowe were the first to observe this enhancement in the 
reaction pd — > 3 HcX by measuring the momentum spectrum of the detected 3 He ejectile. This 
spectrum translates into a 3 He- missing mass spectrum, which corresponds to the 7T7r-invariant mass 
spectrum for X = tttt. In subsequent bubble-chamber [2, 3] and single- arm magnetic spectrometer 
measurements [4, 5, 6, 7, 8, 9, 10, 11] this enhancement was observed in doublc-pionic fusion 
reactions leading to d, 3 He and 4 He, if an isoscalar pion pair was produced. However, such an 
enhancement was not seen in fusion reactions leading to deutcron and triton, if an isovector pion 
pair was produced. These results led to the conclusion that this effect only appears in the double- 
pionic fusion in combination with the production of an isoscalar pion pair. 

Recent exclusive and kinematically complete measurements [12, 13, 14] of the pn — > dTr°ir° 
reaction from threshold up to i/s = 2.5 GeV revealed a strict correlation between the appearance 
of the low-mass enhancement in the 7T7r-invariant mass spectrum (ABC effect) and a narrow 
resonance-like structure in the total cross section of this reaction. This surprising structure peaks 
at m = 2.37 GeV, i.e. roughly 100 MeV below 2mA, where the conventional excitation of a 
AA system peaks. Early attempts to theoretically understand the ABC phenomenon favored 
exactly this conventional i-channcl process for its understanding [15, 16, 17, 18, 19, 20]. The 
recent high-statistics data do not only show that the reaction peaks far below the conventional 
AA process, but that also the narrow width r = 70 MeV of the observed resonance- like structure 
is three times less than expected from the conventional process. 

In fact, all pp- induced, i.e. isovector two-pion production processes including the double- 
pionic fusion reaction pp — > d-K + -K a [21] show neither a pronounced ABC effect nor a narrow 
resonance- like structure in their total cross sections [22]. They rather are well understood by 
conventional ^-channel nucleon excitations. From this follows that any unconventional process 
in connection with the ABC effect has to be associated not only with the isoscalar 7T7r system 
but also with the isoscalar NN system. This hypothesis is supported by the finding that also 
in the double-pionic fusion to 4 He the appearance of the ABC effect is strictly connected with a 
resonance-like structure in the total cross section [23, 24]. 

1.2 New Experiments and Results 

Since the only other doublc-pionic fusion reaction exhibiting an ABC effect is the one leading 
to 3 He, this reaction is ideally suited to test this hypothesis. Hence corresponding exclusive 
and kinematically complete measurements have been performed at the WASA-at-COSY facil- 
ity [25]. These measurements covering the energy range from threshold up to the A A region 
have been carried via the reactions pd -^ 3 Hc7r°7r° at T p = 1.0 GeV and via the quasifree process 
dd -^- 3 }Ieir TT + n spec t a tor, which covers a wide range of energies due to the nucleons' Fermi motion 
in the target. The experiments are complemented by the recent CELSIUS/ WAS A measurement 
at T p = 0.875 GeV [26]. The data analysis is still in progress. As a preliminary result we show 
the energy dependence of the total cross section in Fig.l, where also the results from previous 
experiments at CELSIUS [26, 27] and COSY [28] are plotted. Again we find a resonance-like 
behavior of the total cross section connected with the appearance of the ABC effect in the tttt 
-invariant mass spectrum. 
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Figure 1.1: Preliminary results for the total cross section for the double-pionic fusion to 
3 He obtained from this experiment on pd — ^ 3 He7r°7r° at T p = 1.0 GeV (yellow square) and 
dd — » 3 He7r°7r + n specta tor for the beam energy Td — 1.4 GeV (black dots) in comparison with 
previous results from PROMICE/WASA [27], CELSIUS /WAS A [26] and COSY-MOMO [28]. 
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1.1 Introduction 

Two-pion production in NN collisions gives access to the investigation of nucleon excitations in 
course of the production process. In particular, it provides the unique possibility to study the 
mutual excitation of the two colliding nucleons into one of their excited states. The simplest such 
process is the excitation of a AA system by i-channel meson exchange in the collision process. 

In case that the two-pion production process leads to a fused nucleus in the final state the 
so-called ABC effect - a peculiar low-mass enhancement in the 7T7r-invariant mass spectrum - 
shows up, if the pion pair is of isoscalar nature [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16]. 
Recently it has been demonstrated [17, 18] that the ABC effect in the pn — > dir°ir° reaction is 
strictly correlated with the appearance of a pronounced resonance- like structure in the total cross 
section. This structure with maximum at y/s = 2.37 GcV peaks roughly 100 MeV below the 
maximum of the conventional ^-channel AA process. Moreover it has a width of only 70 MeV, 
i.e. about three times less than the conventional process in its energy dependence. 

In a recent work is has been claimed that a kind of ABC effect is also present in the pp — > ppir°Tr° 
reaction, when the two emitted protons have very small relative momenta representing a kind of 
quasi-bound 2 He system [19]. 

In order to clarify the situation and to contribute to a detailed understanding of the two-pion 
production process, we have studied all different two-pion production channels in pp collisions 
constituting the isovector part of the NN and NNirir systems. 

1.2 Experiments and Results 

In order to complement other data sets [19, 20, 21, 22, 23, 24, 25] exclusive and kinematically com- 
plete measurements of the reactions pp — > pp7T°Tr° [26, 27, 29], pp — > ppir + ir~ , pp — > pnir + Tr° [27], 
pp — > nmr + TT + [30] and pp — > dn + ir° [28] were carried out at the CELSIUS storage ring using 
the 47r WASA detector setup including the pellet target system [31]. The pp — ► ppTT°n° and 
pp — > ppir + 7r~ reactions were evaluated from close to threshold up to T p = 1.36 GeV, the other 
channels only at T p = 1.1 GeV. 

Forward going charged pions, protons, neutrons and deuterons were detected in the forward 
detector and identified by the AE-E technique. Charged pions as well as gammas from w° decay 
have been detected in the central part of WASA detector, which contains a mini drift chamber 
surrounded by thin superconducting magnet and the electromagnetic calorimeter. With the ex- 
ception of neutrons the four-momenta of all emitted particles have been measured giving rise to 
1-6 overconstraints for the subsequent kinematical fit. 

1.2.1 (7T7r)/ = o : pp — > pp7T°7T° 

In Fig. 1, left, the total reaction cross section is shown. It keeps rising from threshold up 
to T p w 1 GeV, where the total cross section levels off and proceeds only slowly rising until 
1.2 GeV. Thereafter it continues steeply rising. Whereas the low-energy region is due to the 
Roper resonance [26] , the renewed rise at higher energies can be associated with the dominance of 
the AA excitation. In Fig. 1, right, the 7r°7r°-invariant mass spectrum is shown for T p = 1.2 GeV, 
the energy which is close to the maximum cross section in the pn — > d7r°7r° reaction. In fact 
we also observe a slight low-mass enhancement, however, in our case it is well accounted for by 
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the conventional calculations. The solid lines in Fig.l show calculations based on the Valencia 
model [32] , where Roper and A A processes have been modified [29] . 



.£5 




T p [GeV] M A „ [MeV/c 2 ] 

Figure 1.1: Energy dependence of the total cross section (left) and 7r°7r°-invariant mass spectrum 
at T p = 1.2 GeV (right) of the 7r°7r° reaction. Solid lines show calculations with the modified 
Valencia model, see Ref. [29] 

1.2.2 (7T7t)j = i : pp — > dir+ir 

In the 7r + 7r°-invariant mass spectrum of this fusion reaction no ABC effect is observed [28]. Due to 
Bose symmetry the isovector ir + ir° system must be in relative p-wave. Both differential distribu- 
tions and the energy dependence of the total cross section are well described by the conventional 
AA process. Note that the Roper process is strongly suppressed in this isovector ttit channel even 
at energies close to threshold. 

1.2.3 (7T7r)/ = 2 : pp — y nn7r + 7r + 

This isotensor tttt channel is special, since single baryon excitations can only contribute, if they 
originate from the excitation of a higher-lying A state. Hence it was expected that the AA process 
provides the main contributions to this channel. However, as demonstrated in Ref. [30], a large 
contribution has to come also from the excitation of a higher-lying A state, favorably the A(1600). 



1.3 Conclusions 

Summarizing we find all data for the two-pion production initiated by isovector pp collisions 
to be well accounted for by conventional i-channel processes involving single or double nucleon 
excitations. This includes also double-pionic fusion processes to d and quasi-bound 2 He. This 
situation is fundamentally different from that in the two-pion production initiated by isoscalar 
TV TV collisions, where the ABC effect has been shown to be correlated with a narrow resonance- like 
structure in the total cross section. 
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1.1 Introduction 

The High Acceptance Di-Electron Spectrometer (HADES) installed at GSI Hclmholtzzcntrum fur 
Schwerionenforschung in Darmstadt was designed to investigate dielectron production in heavy-ion 
collisions in the 1-2 AGeV beam kinetic energy range. The main goal of the HADES experiment 
is to study the properties of hadrons inside the hot and dense nuclear medium via their dielectron 
decays [1, 2, 3]. 

One specific issue of heavy- ion reactions in that regime is the important role played by the bary- 
onic resonances, which propagate and regenerate due to the long life-time of the dense hadronic 
matter phase. The A(1232) resonance is the most copiously produced but with increasing incident 
energy, higher lying resonances also contribute to pion production. A detailed description of the 
resonance excitation and coupling to the pseudoscalar and vector mesons is important for the 
interpretation of the dilcpton spectra measured by HADES. Baryonic resonances are indeed im- 
portant sources of dileptons through two mechanisms: their Dalitz decays (e.g. A/ N* — > Ne + e~ ) 
and the mesonic decay with subsequent dielectron production. 

Simultaneous measurement of one and two-pion production channels can bring information on 
the baryonic resonances excitation. Two-pion production, in particular, is an important subject 
since it connects tttt dynamics with baryon and baryon-baryon degrees of freedom. In our energy 
regime, AA excitation is the leading process [8, 9]. We report on one- and two-pion production 
in pp and np reactions investigated with HADES in exclusive measurements for the beam kinetic 
energies 1.25 and 2.2 GeV. 



1.2 HADES experiment 

The HADES detector [4], as shown in Fig. 1.1, consists of 6 identical sectors covering the full 
azimuthal range and polar angles from 15° to 85° with respect to the beam direction. Each sector 
contains: A Ring Imaging CHcrcnkov (RICH) detector used for electron identification; two sets of 
Mini-Drift Chambers (MDC) with 4 modules per sector placed in front and behind the magnetic 
field to determine momenta of charged particles; the Time Of Flight detectors (TOF/TOFINO) 
and the Pre-Shower detector improving the electron identification. For reaction time measurement 
in heavy-ion reactions, a diamond START detector is located in front of the target. An (e + ,c~) 
invariant mass resolution at the u peak of ~ 2.7% and a momentum resolution for protons of ~ 
3% can be achieved. The first level trigger is obtained by a fast multiplicity signal coming from 
the TOF/TOFINO wall, combined with a signal from the START detector, while the second level 
trigger is made by using the informations from the RICH and Pre-Shower to select the lepton 
candidates. 

In the pp experiment at 1.25 and 2.2 GeV kinetic energy, the intensity was about 10 7 parti- 
cles/second, and a target filled with liquid hydrogen was used. The START detector was not used 
in this run because of the too high rate of secondaries produced by the START detector itself. 
Thus a specific algorithm has been developed to calculate the time of flight and then to identify 
the charged particles. 

*IPN Orsay, sudolSIpno.in2p3.fr 
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Figure 1.1: (Color online) Schematic view of the HADES detector as implemented in the simula- 
tion. Simulated particle tracks are shown as well. 

1.3 Results 

As mentioned in the introduction, it is important to constrain the models predicting the cross 
section and the production mechanisms of the A resonance. In order to do that the hadronic 
channels (pp — > np7r + and pp — > pp7T°) have been measured and studied in parallel to the leptonic 
channels. 

The pp — > np7r + channel is studied using the reconstruction of the undetected neutron. The 
reaction was selected first by the charged particle identification based on momentum and recon- 
structed time of flight, then a (p,7r + ) missing mass cut was imposed around the neutron mass. 
This cut efficiently suppresses the background coming from misidentified protons and two-7r con- 
tributions. 

The pp7r° reaction was extracted from the (p,p) missing mass spectrum by subtracting the 
background under the ir° peak after suppression of the pp elastic and two-pion contribution. 

The absolute normalization of the data has been achieved by a simultaneous measurement of 
elastic scattering for which the cross section is known. Error bars include statistical and systematic 
errors. 

Experimental data have been compared with the original and modified version of the resonant 
model [5] using model A and model B (for details see Table 1.1). A comparison to different spec- 
tra has been made. In particular the neutron angular distribution in the center of mass system 
measured in the pp — >pmr + channel, has been investigated, since it is sensitive to the angular 
distribution of A resonance production (dominant process is pp — » nA ++ ). The acceptance cor- 
rected experimental spectra are presented in Fig. 1.2 (black points) together with slightly modify 
resonance model "model A" (for details of the modification, see Tab. 1.1). The distribution shows 
strong forward/backward peaking, as expected due to the peripheral character of the A resonance. 
The original Teis model (black solid line) shows underestimation at cos9 n ~ 0, which is improved by 
changing cutt-off parameter (A w ) (black dashed line). In addition in Fig. 1.2 comparison between 
HADES and WASA measurement is presented. WASA data were taken from [7] and normalized 
to the cross-section obtained by HADES. One can see reasonable quantitative agreement, keeping 
in mind, that the WASA data are presented without systematic errors. 

Fig. 1.3 exhibits the projection on the (p,7r°) (Left), (p,7r + ) (Middle) and (n,7r + ) (Right) 
invariant masses for T p = 1.25 and 2.2 GeV. The data are corrected for reconstruction efficiencies. 

The prominent peak of Mj n „(p, 7r + ) and Mi„„(p,7r ) around 1.23 GeV/c 2 confirms that most 
of the pions are produced via A decay, which is consistent with the resonance model . However, 
the data present some deviations with respect to the model A, which is close to the original Teis's 
model ([5]). This, together with the above mentionned failure in the description of the angular 
distribution motivates some changes in the model. A better description of the yield and shape 
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Figure 1.2: Angular distribution of neutron in centre of mass system after acceptance correction 
for the pp— s-pii7r + reaction. (Left) Comparison of the HADES data to simulation based on [5] and 
model A. Both simulation curves are normalized to reproduce the integrated experimental yield. 
(Right) Comparison of HADES and WAS A data [7]. 



pp -> pp7i° pp -> np7i + 




M inv (p,7i°) [GeV/c 2 ] M inv (p, K + ) [GeV/c 2 ] M inv (n,;t+) [GeV/c 2 ] 



Figure 1.3: (Color online) Invariant mass distributions measured of (p,7r°) (Left), (p,7r + ) (Middle) 
and (n,7r + ) (Right) measured in pp —> pp7r° and pp — > np7r + at 1.25 GeV (top row) and 2.2 GeV 
(bottom row). HADES data (full dots) are compared within the experimental acceptance to the 
predictions of the modified resonance model (see text). 

of the missing mass spectra is indeed obtained after some modifications which are included in 
the model B (see Tab. 1.1 for details). Due to the dominance of the Delta(1232) resonance, the 
two channels are correlated by isospin symmetry: cr(pp — > nA ++ ) ~ 3(pp — > pA + ), leading to 
"■(PP np7r + ) <~ 5 (pp — » pp7r°). N* contribution seems also reasonable and the invariant mass 
distributions are rather well reproduced by the pp — > nA ++ and pp — > pA + simulations. At 
2.2 GeV, the contribution of higher lying resonances is clearly visible at high invariant masses. 

Furthermore, as another important platform for studying resonances properties, double-pion 
production in nucleon-nucleon collisions is also of interest. Preliminary results obtained for the 
analysis of 2 channels: (1) pp — ¥ ppTr +r K~ and (2) pn — > pn tt + tt~ have been presented in comparison 
to model predictions [5, 6] including double-A and N(1440) excitation. Both channels were selected 
by identifying all charged hadrons involved in the reaction and applying a missing mass cut 
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Tp/Model 


Model A 


Model B 


1.25 [GeV] 


[5], (1) 


Model A, (2) 


2.2 [GcV] 


[5], (1) 


Model A, (3) 



Table 1.1: Summary of the modifications introduced to the resonance model [5] in order to better 
fit the data: (1) pp and pn final state interactions have been implemented using the Jost function 
and adjustment of the cut-off parameter (A ff ) entering the ttNA and ir NN vertex form factor 
has been changed from 0.63 to 0.75 GeV, (2) the A production angular distribution has been 
further adjusted to reproduce the measured neutron angular distribution, (3) the cross section of 
the higher lying resonances has been increased and a non-resonant contribution has been added. 

(Mp p7r + W - in case (1) and M JOT + W - in case of (2)). It turns out that the invariant mass (M n + 7r -) 
and the opening angle in CM (cos<S„.-t-„.- ) of the pion pair are the most sensitive distributions to 
the different model contributions. Fig. 1.4 exhibits an experimental distributions of the invariant 
mass and opening angle of tt + it~ in comparison to the pure phase space (PHPS) calculations. 
The PHPS calculations has been normalized to the experimental yield. It is seen, that both 
experimental distributions deviate from PHPS calculations. Some enhancement can be observed 
at low 7T + 7r~ invariant masses, which is not present in the simulation with the PHPS only. This 
peak at small invariant masses appears in the models [8, 9] as being due to the AA excitation and 
to the decay channel of the Ropper resonance N(1440) into NA (N* ->N A). 




0.2 0.3 0.4 0.5 0.6 -1 -0.5 0.5 1 



M^. [GeV/c 2 ] cos8 n+n . 

Figure 1.4: Distribution of the tt + it~ invariant mass M n + V - (left) and the 7r + 7T~ opening angle 
in centre of mass 5 n + w - (right) for pp — ► ppir + ir~ reaction at beam energy Tkin — 1-25 GeV are 
shown in comparison to the phase space calculation. Both spectra are inside HADES acceptance 
and efficiency. The simulation is normalized to reproduce the measured experimental yield. Only 
statistical errors are shown. 



1.3.1 Conclusion 

The aim of measuring one and two-pion production in elementary reactions with HADES are three 
fold: (1) to study the detector response, and cross-check the analysis of dilepton channels, (2) test 
the ingredients of the transport models used for the dielectron production, which are based on 
resonance models and (3) to obtain a parametrization of meson and baryon resonance production 
for dielectron channels analysis. 

HADES provides high statistics data in one and two-pion production channels, which comple- 
ments measurements of previous experiments. For the one pion channel, which is dominated by A 
resonance, the data have been compared to the resonance model usually used in transport models, 
and modifications have been proposed to better reproduce the data. The ongoing analysis of the 
two-pion production and comparison to the models will bring constrain on double-A production 
as well as N(1440) production and decay mechanism. 
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1.1 Introduction 

The High Acceptance Di-Electron Spectrometer (HADES) [1] operates at the GSI Helmholtzzen- 
trum fur Schwerionenforschung in Darmstadt, Germany. One of the main physics goals of HADES 
is to investigate spectral modifications of light vector mesons in strongly interacting matter via 
dilepton (e + e~) decay channel. The HADES collaboration has set up an experimental program 
to measure dilepton spectra in p+p and p+n collisions [2]. In year 2007 electron-positron pair 
production has been measured in elementary collisions of two protons with 3.5 GeV projectile 
kinetic energy. One of basic observables in this measurement is inclusive e + e~ mass distribution , 
shown in Fig. la after Combinatorial Background (CB) subtraction. The expected e + e~ produc- 
tion channels are given by Dalitz decays of 7r°, T), to mesons and A(1232) resonance. In the high 
mass region 2-body decays of p/ui play a role. Indeed, the experimental data can be described by 
PLUTO [3] simulation (see Fig. la) assuming aforementioned components, but not in the mass 
region below vector meson pole {Mf n £ £ (0.5 — 0.7)). However, only Dalitz decay channel for 
baryonic resonances is included for A(1232) here. 



^10" 1 t i ' | ' ' i | ' i i | i ' i | i ' i | ' ' ' 4 




Fig. la) Efficiency corrected inclusive invariant mass 
distribution of dielectrons inside the geometrical ac- 
ceptance of HADES. e+e — - signal pairs (black 
dots), after combinatorial background subtraction for 
p(3.5GeV)+p collisions with coresponding simulation 
coctail components. The data are normalized to the 
measured p+p elastic events. 




Fig. lb) Efficiency corrected invariant mass distribu- 
tion of dielectrons for exclusive ppe + e~ channel inside 
the geometrical acceptance of HADES after combinato- 
rial background subtraction. The data are normalized 
to p+p elastic events and compared to simulation cock- 
tail assuming sources given in the legend. 



In order to study contribution of higher lying baryonic resonances ppe + e~, ppn , pnir + , 
ppcu(r]) — > 7T + 7r~7r° exclusive channels have been analysed. The ppe + e~ channel has been se- 
lected by condition on pe + e~ missing mass ~ 0.938 [MeV/c 2 ] and the respective e + e~ invariant 
mass spectrum is shown in Fig. lb. The mass region below vector meson pole is not described 
by channels including ui, rj decays which are constrained by hadronic channels mentioned above 
(see [7]). This means that Dalitz decays of baryon resonances must contribute in this area. In 
order to estimate production cross sections of the resonances which decay in the pmr + and ppir 
final states, the results have been compared to Monte-Carlo calculations based on the resonance 
model [4] assuming incoherent sum of various resonances. 

* Jagiellonian University, Cracow; adybczakSgmail . com 
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To convert resonances cross sections into e + e~ yield one should remember that the respective 
partial decay widths are mass dependent. Two models of T^^pg+g- (M) of ([5], [6] are applied in 
the simulation. 

1.2 Analysis of pn7r + , ppn final states 

In both final states only events with two positive tracks were selected. Protons and positive 
pions were distinguished using measured momentum and particle velocity of reconstructed tracks. 
Channels were selected via condition on missing mass of neutral particles (n and 7r°, respectively). 
Signal was obtained as a number of counts above fitted continuous background. In the case of ppir 
discrimination on clastic scattering and multipion production was required additionally to enhance 
purity of the signal. Resonances cross section were estimated by simultaneous fit of invariant 

masses (Mf* v , ^Tnv ; ^fnv) and angular distributions {cos&q M , cosO'^m ; cos ^cm) m severa l 
bins of the N — n invariant mass. Iterative procedure was applied: In the first step a flat angular 
distribution has been assumed for resonances production. Since experimental distributions show a 
clear forward-backward peaking, in the next step a modified distribution of the form da/dt ~ A/t a 
has been used. The t, A, a are t-Mandelstam variable and fit parameters depending on resonance 
mass, respectively. Furthermore, the isospin relations between cross sections of various resonance 
have been preserved. 

1.2.1 Acceptance corrected pnir + , pp(tt°) results 

In Fig. 2a, the results obtained after acceptance/efficiency corrections assuming aforementioned 
model are presented. The pnn + channel allows to distinguish doubly and singly charged reso- 
nances. 
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Fig. 2a) Acceptance and efficiency corrected spectra of 
pir + , nir + invariant mass with corresponding simulation 
components. 



Fig. 2b) Acceptance and efficiency corrected spectra of 
p7r + , nn + angular distribution in the centrum of mass 
of with corresponding components obtained from fitting 
procedure described in the text. 



1.3 Conclusions and outlook 

The applied resonance model of baryon resonance production can fully describe production of pions 
in the two ppir and pnir + channels. Obtained cross sections for the various resonances fulfill all 
isospin relations. The derived cross section for A(1232) production is in a good agreement with 
observations of other experiments. Estimated cross sections for the total tt° and ir + yields are 
also in a good agreement with previous measurements. The obtained cross sections for resonance 
production will be further used to calculate dielectron production in the ppe + e~ channel. 
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Fig. 3a) Acceptance and efficiency corrected spectra of 
pit invariant mass with corresponding simulation com- 
ponents. 




Fig. 3b) Acceptance and efficiency corrected spectra 
of pn° angular distribution in the centrum of mass of 
resonance with corresponding simulation components. 
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1.1 Introduction 

The High Acceptance Dielectron Spectrometer (HADES) is a second-generation experimental 
setup located at the heavy-ion synchrotron SIS18/GSI Darmstadt. Its powerful particle iden- 
tification capabilities and excellent momentum resolution allows for simultaneous measurements 
of hadronic (p, ^ , K^) and electromagnetic (e + e~) final states over a wide angular range 
(18° < 9 < 85°) in proton, deuteron and HI induced reactions [1]. One of the main goals 
of the HADES experiments is to achieve a detailed understanding of dielectron emission from 
hadronic systems at moderate (1 — 4 AGeV) bombarding energies. A particular feature of this 
energy domain is the important role of baryonic sources as N + N bremsstrahlung and resonance 
Dalitz decays (i.e A(N*) — > Ne + e~) which are poorly known. The results obtained on electron 
pair production in elementary N + N collisions represent a doorway for the understanding of 
pair production in HI collisions where additional in-medium effects, as for example hadron mass 
modifications, are expected to occur [2, 3, 4]. In order to shed more light on the aforementioned 
processes, HADES performs a combined study of pion and dielectron production in p+p and d+p 
reactions. The results from pion production are given in the contributions of M. Gumberidze [5] 
and A. Dybczak [6]. In this report, a survey on important results on e+e~ pair production is 
given. 

1.2 e + e~ results from p + p and quasi-free p + n collisions 

Results from the pioneering DLS experiment at Bevelac on inclusive e+e~ production in p + p 
and p + d reactions indicate a strong isospin dependence of the pair production in the 1 — 4 AGeV 
region [7]. The cross section ratio of pair production is significantly larger than 2 which was 
observed at E^ in > 2 GeV and could only be qualitatively described by model calculations. 

Fig. 1 shows inclusive e + e~ invariant mass distributions measured in p+p (left-top) and quasi- 
free p + n (left-bottom) collisions at 1.25 GeV [8] inside the HADES acceptance. Quasi-free n + p 
reactions were selected from d + p — >• e+e~ + X reaction on trigger level by the detection of a fast 
spectator protons from the deuterium break-up in a dedicated forward hodoscope wall at small 
angles w.r.t. beam direction. As one can see, a large (at least one order of magnitude) excess of the 
pair yield is visible for the p+n reaction channel at higher invariant masses. The data are compared 
to simulated pair distributions calculated with the Pluto event generator [9] assuming a dominant 
role of the tt° — > e + e~ +7 and A(1232) — > e + e~ +N Dalitz decays. For both reaction channels the 
measured yield in the tt° Dalitz decay region is very well reproduced taking into account the inclu- 
sive 7T° production cross section from the resonance model assuming A(1232) dominance [8]. This 
assumption is also corroborated directly by the HADES measurements of the dominant exclusive 
channels pp — » ppir and pp — > pnn + (see [5]). To model A(1232)+(°) — > p(n)e + e~ the differential 
partial decay width dT Ne+e-(M/±) / 'dM e+e from [10] has been utilized, assuming a constant mag- 
netic transition form factor Gm = 3.02 ± 0.03 fixed by photo-production experiments (the other 
two form factors, electric and Coulomb, are assumed to be zero) [8]. A possible modification of 
the magnetic transition form factor due to intermediate vector mesons was simulated assuming 
the two-component quark model [11] and is visualized by the shaded area. As expected, it leads 
to a slightly enhanced yield for high pair masses. Whereas the p + p spectrum is very well de- 
scribed by both contributions, the description of the p + n channel is not satisfactory even if the 
rj — > e+e~ + 7 decay is included. The latter one is fixed by precise data from WASA at CELSIUS 
on near threshold production in pn — > dn and pn — > pnn [12]. 
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Figure 1.1: Differential cross sections for dielectron production as a function of the invariant mass 
as measured in p + p, quasi-free p + n reactions at 1.25 GeV (left) and p + p at 3.5 GeV (right). 
The distributions are compared to model calculations assuming various hadronic sources (for more 
details see the text). 



Next, the obtained results are compared to predictions of the OBE model [13], where resonant 
(A) and non-resonant ("quasi-free" N + N bremsstrahlung) contributions are added coherently. 
The result of the simulation is shown in Fig. 1 as solid black line. The yield calculated in this 
approach overestimates (underestimates) measured distribution in the p+p (p + n) reactions. One 
should add here that more recent calculations of [14] based on similar effective Lagrangian are 
successful in description of the p + p reaction channel and also come closer to the data measure 
in the p + n reaction channel. For the latter case, the implementation of the pion electromagnetic 
form factor in graphs including exchange of charged pions turns out to be decisively to produce 
e + e~ yield enhancement visible in the p + n reaction channel at higher masses. 

The right part of Fig. 1 shows the e + e~ invariant mass distributions measured in p+p collisions 
at 3.5 GeV [15]. In contrast to lower energy, the spectrum extends to much higher masses due 
to the rapidly changing excitation function of dielectron production. A clear ui peak with mass 
resolution AM/M ~ 2% can be seen in e + e~ channel, for the first time at such low energy. Below 
the vector meson mass region dielectron production is dominated by the n and 7r° Dalitz decays. 
The comparison of the experimentally measured invariant mass distribution to a incoherent sum 
of known hadronic sources is also shown in this figure. Calculations have been performed with 
PYTHIA + PLUTO codes [16, 15]. The simulated cocktail describes the data reasonably well 
except for the mass range around 0.55 GeV/c 2 where the yield is underestimated. The latter 
deviation is not too surprising since in this energy region contribution from the resonances (in our 
calculations only from A(1232)) can be expected. Furthermore, the virtuality of the photon reaches 
quite high values and one can expect shape modifications due to unknown mass dependency of the 
A(1232) transition form factor. Moreover, contributions of higher lying (A, TV*) resonances has 
not been included, yet. The respective contributions can be better constrained by means of the 
exclusive dielectron and hadronic channels: ppe + e~ , ppir , pnir + , ppn and ppco inside the HADES 
acceptance. The strategy is to fix the cross sections for resonance (A, 7V*( + '°)) production from the 
hadronic final states and convert it to ppe + e~ yields by means of dT Ne+e-{M res ) / dM e+e - given 
by available calculations [17, 18]. Many of these hadronic channels have already been analyzed. 
The channels with pion production is discussed in the contribution of A. Dybczak [6]. 
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1.3 Summary and outlook 

HADES measured diclcctron production in proton, deuteron and HI induced reactions in the 1 — 4 
AGeV energy region. Diclcctron production in N + N collisions is a mandatory prerequisite for 
the understanding of pair production in HI collisions and already observed in-medium effects. The 
interesting feature of the pair production in low-energy domain is the important role of baryonic 
e+e~ sources , i.e. resonance Dalitz decays and N + N bremsstrahlung. The description of the 
respective pair production process is directly linked to the questions of how the virtual massive 
photon couples to baryons and how do the respective transition form factors look like. The results 
already obtained by HADES in p+p (1.25, 2.2, 3.5 GeV) andp+n (1.25 GeV) indicate a sensitivity 
to such effects. Further measurements are however mandatory to pin down respective effects. In 
this context, future experiments with pion beams planned with HADES at GSI can be extremely 
helpful. In contrast to proton induced reactions, resonances can be directly excited in the s channel 
and separated by means of excitation function and angular measurements [19]. 
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1.1 Introduction 

The HADES [1] (Acceptance Dielectron Spectrometer) collaboration demonstrated, using pp and 
quasifree pn reactions [3] that baryonic resonances play an important role in dielectron emission at 
SIS accelerator energies (1-2 AGeV). This is first due to their mesonic decays and the subsequent 
direct dilepton (e.g. p/u> — > e + e~) or Dalitz (also called conversion) decay (7r°/ry — ► r ye + e^ or 
co — > Ti e + er) modes of these mesons. Baryonic resonances are also expected to contribute directly 
to dilepton emission via their Dalitz decay modes (N* /A — > Ae+e~). (see Piotr Salabura's talk 
for more details) 

The Dalitz decay of baryonic resonances can be studied in a more precise way using the GSI pion 
beam. In n~ + p reactions, the production mechanisms are better controlled and all resonances 
are produced at fixed mass, in the s channel. In addition to electromagnetic channels, meson 
production channels (irN — > irirN, irN — > r]N, ttN — > p/ujN, irN — > KA ...) can also bring useful 
information to study resonance properties, via partial wave analysis (PWA). There is indeed a 
big interest from theoreticians, for new precise differential cross section data for various reaction 
channels and at various incident momenta. 

These experiments present some technical challenges. This secondary pion beam is produced 
on a beryllium targets in C+Bc or p+Be reactions. To achieve pion beam intensities of the 
order of 10 6 /s, primary carbon or proton beams need to reach an intensity close to space charge 
limit of the accelerator. To produce pions with momentum from 0.7 to 1.5 GeV/c, the reaction 
C+Be@2.0GeV should be chosen but to produce pions with momenta from 1.5 to 2.2 GeV/c 
with sufficient intensity, p+Be@3.5GeV has to be used. The acceptance of the beam line is 
Ap/p w ±4%. For exclusive studies momentum reconstruction is needed. Furthermore, the 
bacground due to secondary particles produced by the beam halo has to be suppresed. Therefore, 
an additional event by event pion in-beam tracking is required. This task is handled by two in- 
beam silicon detectors X1,X2. Trackers are two double sided, 2x128 strip, 300^m thick silicon 
detectors manufactured in high radiation hard technology. This system provides a momentum 
resolution of Ap/p ~ 0.1% which is sufficient for exclusive analysis studies. A diamond detector 
placed in front of the target complements the beam tracking and will be used as start detector. 

1.2 Results 

The poster shows simulations and feasibility studies for future experiments with the pion beams. 
First the total cross section to produce baryon resonances up to A°(1950) [4] where calculated. 
Next step consisted in Dalitz decay simulations of the baryonic resonances , where two electro- 
magnetic decay models of baryons were compared. The first one where they are treated as point 
like [5] and the other one using the eVMD (extended Vector Meson Dominance) model for the 
electromagnetic form factors [6] . In addition, the production and decay of p, to and r/ mesons is 
included in the simulations. Finally, the HADES spectrometer acceptance for dilepton pairs was 
determined. This allowed to produce inclusive invariant mass spectra inside HADES acceptance 
and estimate count rate per one week of experiment at two different energies (fig. 1 right). 

A second type of simulations was devoted to two pion production in tt~ +p reactions. Like in 
the case of dilepton analysis, cross sections and branching ratios were determined for bench mark 
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Figure 1.1: Left: Invariant mass spectra of -k + h~ pairs from simulations of 7r~ + p — > n7r + 7r~ 
reactions at 0.8 GeV/c and estimated count rates per 1 week of experiment. Right: Inclusive 
invariant mass spectra of e + e _ pairs from simulations of tt~ + p reactions at 0.8 GeV/c and 
estimated count rates per 1 week of experiment. 

channels. Next the HADES acceptance for 7r + 7r~ pairs and the inclusive Tr + n~ invariant mass 
spectra with estimated count rate per one week of experiment have been presented. 

1.3 Perspectives 

Pion beams at GSI are available for many years and used by other groups (e.g. FOPI coll.). After 
the heavy ion campaign scheduled for nearest future, the HADES collaboration is preparing a 
physics program for experiments with pion beams. The simulations will be further developed by 
including Dalitz decay of higher lying resonances and by studying exclusive channel analysis like 
7r~p — > ne + e~ , 7r~p — > n7r + 7r~, n~p nr]. The aim is to define in detail the physics case and 
beam momenta for these future experiments. 
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1.1 Introduction 

A data sample of integrated luminosity C = 242.5 pb -1 collected at yfs = 1 GeV with the KLOE 
detector at the DA^NE <ji-factory has been analyzed in order to study the e+e~ — > e + e~X reac- 
tions, X being the 77 meson or the 7r°7r° state produced in the scattering of two quasi-real photons. 
The measurement is done with final state e going along the beam pipe and escaping detection. In 
this case the virtual photons are quasi-real and the overall cross section can be factorized as the 77 
subprocess cross section times a 77 luminosity function, cF e + e -^ e + e - x — J CF it^x{w)'^^ L dw [1]. 



1.2 77 — > 7] 

Two independent analyses have been performed looking for 77 — > 77 production with 77 — > 7r + 7r~7r° 
and 77 — > 7r 7r°7r°. In both analyses the signal has been simulated with a Monte Carlo generator 
based on the complete matrix clement calculation which allows the full phase space generation [2] . 

77 ~~ ► V ~ ► 7r + 7r _ 7r°. Events are selected asking for two photons from a ir° decay and two 
tracks with opposite curvature coming from the interaction point, the 7^ mass being assigned to 
each track. A kinematic fit is performed using the 77 — > -K + n~n° decay hypothesis; the fit returns 
a x 2 value and events with \ 2 > 20 are rejected. Further cuts are applied to suppress processes 
with photons and e + e~ in the final state. The number of signal events is extracted from the fit to 
the squared missing mass distribution M^ iss using Monte Carlo distribution for both signal and 
backgrounds (Fig. 1.1). From the number of signal events given by the fit one obtains the cross 
section er(e+e~ — > e + e~r/, y/s = 1 GeV) = (41.7 ± 4.0 s t a t) pb, the systematic error being under 
evaluation. Work is in progress to extract the partial width T v ^ 71 . 

The main irreducible background is due to e + e _ — > 777 — > Tr + ir~ir°j events with the monochro- 
matic photon (Ej = 350 MeV) lost in the beam pipe. In order to normalize this process 
properly a dedicated analysis has been done using the same C = 242.5 pb -1 data sample, se- 
lecting events with three photons in the final state and performing a kinematic fit asking for 
energy and momentum conservation; improved variables have been used to fit to data distribu- 
tion of the energy of the monochromatic photon. The preliminary result for the cross section is 
a(e+e- -> 777, yfs = 1 GeV) = (0.866 ± 0.009 stat ± 0.093 syst ) nb [3]. 

77 7r 7r°7T . Events are selected with six photons from the interaction point (prompt 

photons) ; the photons are paired choosing the combination which minimizes the x 2 variable 

2 _ {m^-m lj f (m^-rriki) 2 (m OT - m st ) 2 . . 

x 77 — 9 + n 2 + 2 ' y 1 - 1 -) 

a ij a ki a st 

with rriij the two-photon invariant masses and their resolutions. A kinematic fit is performed 
constraining the six-photon invariant mass to the 77 meson mass and asking for the space-time 
relation t — |~F'|/c = 0tobe satisfied for each photon; a cut is then applied on the value of the 
kinematic fit \ 2 ■ Events with tracks in the drift chamber are rejected; a cut is applied on the 
energy of the most energetic photon, to reject 777 events. The squared missing mass distribu- 
tion shows an excess of events well reproduced by the e+e~ — > e + e~r] -4- e + e~7r 7r°7r° Monte 
Carlo. The fit gives normalization for the Monte Carlo signal and the backgrounds distributions 
(Fig. 1.2 left). The preliminary result for e + e~ — > e + e~rj production cross section at yfs = 1 
GeV is a(e + e~ -> e + e~r/, 1 GeV) = (37.0 ± 1.4 s ^ a ^ ± 2.2 S y S ^) pb, to be compared with the result 
obtained studying the 77 — > 7r + 7r~7r° decay channel. From 777 normalization a preliminary value 



*Univcrsita dcgli Studi Roma Tre and INFN, Sezione Roma Tre; plonghiaroma3.infn.it 



119 



PrimeNet Workshop 2011, Jiilich 



I. Prado Longhi 




Figure 1.1: 77 — > 7/ — > 7r + 7r~7r° analysis: squared missing mass distributions for data (points) and 
Monte Carlo normalized according to fit result (solid lines). Colour code: red = 777, light blue = sig- 
nal, blue = W7T°, green = e + e _ 7 (preliminary plot). 




Figure 1.2: 77 — > 77 — > 7r 7r°7r° analysis. Left: squared missing mass distributions for data (points) 
and Monte Carlo normalized according to fit result (solid lines). Colour code: red = 777, light 
blue = signal, blue = lott (preliminary plot). Right: preliminary KLOE result for er(e + e~ — > 777) 
at y/s = 1 GeV, (red point, with statistical uncertainty only) and SND results [4] at several values 
of yfs. 



for the cross section a(e + e — > 777, y/s = 1 GeV) = (0.875 ± 0.018 s t a t ± 0.035 S y S ^) nb is evaluated, 
in agreement with the measurement done in the 77 — > 7r + 7r~7r decay channel. This value is shown 
in Fig. 1.2 (right) together with the SND experimental results [4] as a function of y/s. 



1.3 77 7r°7r° 

The main goal of this analysis is to investigate the low 7r°7r° invariant mass region, just above 
the production threshold, where a contribution from the ct(600) scalar meson as a resonant in- 
termediate state is expected. The main background processes are annihilation reactions to states 
with four or more prompt photons: e + e~ — > KsKl, 777, W7r°, /07, ao7- In addition, the e + e _ — » 77 
process is also considered as a source of background due to possible splitting of the photons energy 
deposits in the calorimeter. Events with 4 prompt photons are selected. The photons are paired 
as in the 77 —> 7/ — >• 7r 7r 7r analysis, and events with > 4 are rejected: the effect of this 
selection is shown in Fig. 1.3 (left). Events with no tracks in the drift chamber are selected; cuts 
on photons energies and on the transverse momentum of the four-photon system are applied; a 
cut is applied on the ratio of the sum of the energies of four photons to the total energy deposit 
in the calorimeter, to reject KgK^ events with large amount of non-prompt energy released in 
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Figure 1.3: 77 — > tt q tt° analysis. Left: two-photon pairs invariant masses before (top) and after 
(bottom) rejection; selected events have both photons pairs invariant masses centered around 
7r mass value. Right: four-photon invariant mass spectrum for data (points with error bars) and 
backgrounds Monte Carlo. Colour code: light blue = KsKl, blue = ljtt , violet = /07, green = 7/7, 
red = 77 (preliminary plot). 



the detector. The four-photon invariant mass spectrum for the selected data sample is shown in 
Fig. 1.3 (right) together with the normalized Monte Carlo background. An excess of events in the 
low invariant mass region signals ir°ir° production processes. 

1.4 Conclusions 

The reported analyses show excess of events with respect to the annihilation processes in the region 
where the signal from 77 or 7r°7r° production via 77 interaction is expected. For the e + e~ — > e + e~rj 
cross section at y/s = 1 GeV two compatible measurements are obtained studying the r; — > 7r + 7r~7r 
and the r\ — > 7r 7r°7r° decay channels. For the 77 — > 7r°7r° process work is in progress to determine 
the signal efficiency and the 77 luminosity function in order to extract the cross section and 
compare it with the only available measurement in literature [5]. 

The forthcoming data taking with KLOE-2 detector, equipped with dedicated small angles 
taggers for e ± in the final state, will give precious additional information on 77 hadron production 
at low energy [6]. 
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1.1 Introduction 

Presently, there are several experiments dedicated to the study of low-energy pion-photon inter- 
actions. At COMPASS, pion Compton scattering is being investigated by means of the Primakoff 
effect in n~ Z — > tt~ Z"f with the aim to extract pion polarizabilities [1], while at KLOE the reac- 
tion 77 — > tt°tt° is being measured in e + e~ — > e + e~7r°7r° [2]. Measurements of the cross section 
for 77 — > tt°tt° at low energies are particularly valuable, since they provide an alternative way to 
access the pion polarizabilities and the absence of Born-term contributions in this channel renders 
the cross section sensitive to the influence of the a resonance. However, in both cases one needs 
theory input to reliably extract the information of interest, as the polarizabilities are defined at 
the threshold for pion Compton scattering, and studying the a even requires analytic continuation 
far into the complex plane. The necessary extrapolations should be performed respecting all avail- 
able constraints, i.e. analyticity, unitarity, crossing symmetry, and gauge invariance. Roy-Steiner 
equations for 77 — > tttt provide a framework that meets all these requirements [3]. 1 

The Roy equations for tttt scattering [5] are a coupled system of partial wave dispersion relations 
that respects analyticity, unitarity, and crossing symmetry of the scattering amplitude. In recent 
years, partial wave dispersion relations in combination with unitarity (and chiral symmetry) have 
been used for high-precision studies of low-energy processes, both in tttt [6, 7] and irK [8] scattering. 
An important application of tttt Roy equations in combination with Chiral Perturbation Theory 
(ChPT) was the precise prediction of the pole parameters of the a resonance [9] 

M a = 441^ 6 MeV, Y a = 544±£gMeV. (1.1) 

The reaction 77 — > tttt provides an alternative to tttt scattering for the excitation of the a. In 
particular, Roy-equation techniques in 77 — > tttt allow us to constrain the <r's two-photon width 
T ajl at a similar level of rigor as M a and r CT based on tttt Roy equations [3] . 



1.2 Roy equations for 7T7T scattering 

Roy equations for tttt scattering are obtained by starting from a twice-subtracted dispersion re- 
lation at fixed Mandelstam t, determining the i-dependent subtraction constants by means of 
crossing symmetry, and finally performing a partial wave expansion. This leads to a coupled sys- 
tem of integral equations for the tttt partial waves tj(s) with isospin I and angular momentum J 



I'=0J'=0jf2 



(1.2) 



where Kjj, are known kinematical kernel functions and the tttt scattering lengths — the only free 
parameters — appear in the subtraction term kj. Assuming elastic unitarity 

e 2iS '->^ - 1 / AM 2 

\mt I J {s)=a{s)\t I J {s)\\ #(*)= 2 . a{s) , a{ S ) = sjl-—^, (1.3) 

(1.2) translates into a coupled integral equation for the phase shifts Sj themselves. 

*Helmholtz-Institut fur Strahlcn- und Kernphysik (Theorie) and Bethe Center for Theoretical Physics, Univer- 
sitat Bonn, hoferichterahiskp.uni-bonn.de 

^Institute of Nuclear and Particle Physics and Department of Physics and Astronomy, Ohio University 
tCONICET - Departamento de Fi'sica, FCEyN, Universidad de Buenos Aires 
x The following presentation involves significant textual overlap with [4]. 
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1.3 Roy— Steiner equations for 77 — >• 7T7T 

Crossing symmetry in this case is less restrictive than for tttt scattering, as it couples 77 — > tttt 
to pion Compton scattering — > 771- , which we will consider as the s-channel process. Roy- 
Steiner equations are then most conveniently constructed based on hyperbolic dispersion relations 
[10]. The resulting system of integral equations couples the 77 — > nn partial waves hj ± (t) to the 
77T — > 77r partial waves fj±(s) (with photon helicitics ±), e.g. 

h I j,-(t)=Nj(t) + ± [ ds'J2Gjj l (t,s')Imf I J , t+ (s') + ^ f dt'J^K^it^lmh'j,^'), 

J 77_1 J 77 

Ml J _1 4M2 J 

(1.4) 

where Nj (t) includes the QED Born terms. Subtracting at t = 0, s = M%, the subtraction 
constants directly correspond to pion polarizabilities. In the once-subtracted case, one needs the 
dipole polarizabilities a\ ± while a second subtraction requires in addition knowledge of the 
quadrupole polarizabilities a.2 ± /3 2 . 

Elastic unitarity is also less restrictive than for tttt scattering, since the unitarity relation is 
linear in hj± 

Im^ i± (t) = c7(t)^ i± (t)tS(t)*. (1.5) 
Below inelastic thresholds the phase of hj± coincides with Sj ("Watson's theorem"). Assuming 
this phase to be known, the equations thus reduce to a Muskhelishvili-Omnes problem for hj± 
[11]. 



1.4 Muskhelishvili-Omnes solution and results for T ail 

To solve the equations for hj ± , we truncate the system at J = 2. Furthermore, we assume the 
amplitudes to be known above the matching point t m = (0.98 GeV) 2 . The solution can then be 
written down in terms of Omnes functions 

m) = ^] M JM-\ (1 . 6) 

We find that the solutions for different partial waves in general do not decouple, e.g. the equation 
for the S*-wave involves spectral integrals over the D- waves as well [3] . This is a new result of our 
dispersive treatment of 77 — > tttt based on Roy-Steiner equations. 

We approximate Im fj ± (s), which at low energies is dominated by multi-pion states, by a sum 
of resonances [12]. Above the matching point we use a Breit-Wigner description of the /2(1270), 
which dominates the cross section at higher energies. Within our formalism [3] we derive a sum rule 
for the 1 = 2 polarizabilities, which — in combination with ChPT results for dipole and neutral-pion 
quadrupole polarizabilities [13] — produces an improved prediction 

(a 2 - /3 2 ) 7r± = (15.3 ± 3.7) • lO^fm 5 (1.7) 

for the charged-pion quadrupole polarizability. This sum-rule result together with the ChPT 
values for the other polarizabilities [13] leads to the "ChPT" prediction for the total cross section 
of 77 — > 7r°7r° depicted in the left panel of Fig. 1.1. The result labeled "GMM" is found when we 
adopt the polarizability values of a recent fit of a two-channel Muskhelishvili-Omnes representation 
to 77 — > nn cross section data [12]. The uncertainty due to the tttt phases represented by the grey 
band is estimated by varying between two recent state-of-the-art analyses based on Roy and Roy- 
like equations [7, 14] . We see that especially for the twice-subtracted version the agreement with 
experiment in the low-energy region is very good. Given sufficiently accurate data on 77 — > 7r°7r° 
such an analysis could be used to extract pion polarizabilities. 

Since we have shown that the a lies within the domain of validity of our Roy-Steiner equa- 
tions [3], this formalism allows for a reliable analytic continuation to the a pole. The main result 
of our analysis is shown in the right panel of Fig. 1.1: there is a correlation between T (J ^ 1 and the 
1 = pion polarizabilities that follows from Roy-Steiner equations and input for the tttt phases 
alone. In combination with the ChPT-plus-sum-rule input for the polarizabilities, we obtain 

r CT77 = (1.7±0.4)keV. (1.8) 
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\ft I GeV (qi - /3i) J=0 / lO-'fm 3 



Figure 1.1: Total cross section for 77 — > tt°tt° for | cos0| < 0.8 1 (left) and T ail as a function of 
the 1 = pion polarizabilities (right). The black line refers to the unsubtractcd case and the 
colored lines to the twice-subtracted version with {a.2 — h) I=Q as indicated (in units of 10 _4 fm 5 ). 
The grey bands represent the uncertainty due to the 7T7t input. The cross corresponds to the 
twice-subtracted case plus ChPT input. 

1.5 Outlook 

The framework presented here should prove valuable for the interpretation of upcoming KLOE 
data for 77 — > 7r°7r°. Especially in combination with the ongoing efforts at COMPASS concerning 
77r~ — > 77r~, this will further improve knowledge of the pion polarizabilities and thus of r CT77 . 

Acknowledgements 

This research was supported by the DFG (SFB/TR 16), the program "Kurzstipendien fur Dok- 
torandlnnen" of the DAAD, the Bonn-Cologne Graduate School of Physics and Astronomy, the 
US Department of Energy (Office of Nuclear Physics), and CONICET. 

Bibliography 

[1] A. V. Guskov, Phys. Part. Nucl. Lett. 7 (2010) 192; Nucl. Phys. Proc. Suppl. 198 (2010) 112. 
[2] The KLOE2 Collaboration, arXiv: 1107.3782 [hep-ex]; I. Prado Longhi, these proceedings. 
[3] M. Hoferichter, D. R. Phillips and C. Schat, Eur. Phys. J. C 71 (2011) 1743. 
[4] M. Hoferichter, D. R. Phillips and C. Schat, arXiv: 1108.4776 [hep-ph]. 
[5] S. M. Roy, Phys. Lett. B 36 (1971) 353. 

[6] B. Ananthanarayan, G. Colangelo, J. Gasser and H. Leutwyler, Phys. Rept. 353 (2001) 207. 

[7] R. Garcia-Martfn, R. Kamihski, J. R. Pelaez, J. Ruiz de Elvira and F. J. Yndurain, Phys. 
Rev. D 83 (2011) 074004. 

[8] P. Biittiker, S. Descotes-Genon and B. Moussallam, Eur. Phys. J. C 33 (2004) 409. 

[9] I. Caprini, G. Colangelo and H. Leutwyler, Phys. Rev. Lett. 96 (2006) 132001. 

[10] G. E. Hite and F. Steiner, Nuovo Cim. A 18 (1973) 237. 

[11] N. Muskhelishvili, Singular Integral Equations, P. Noordhof, Groningen, 1953; R. Omnes, 
Nuovo Cim. 8 (1958) 316. 

[12] R. Garcia-Martm and B. Moussallam, Eur. Phys. J. C 70 (2010) 155. 

[13] J. Gasser, M. A. Ivanov and M. E. Sainio, Nucl. Phys. B 728 (2005) 31; 745 (2006) 84. 

[14] I. Caprini, G. Colangelo, H. Leutwyler, in preparation. 



124 



MEASUREMENT OF THE 
DOUBLE POLARIZATION OBSERVABLE G 
IN 7T° AND 77 PHOTOPRODUCTION 

AN NIK A THIEL* for the CBELS A/ TAPS- Collaboration 

1.1 Introduction 

The excitation spectrum of the nucleon consists of several strongly overlapping resonances with 
different strengths. To disentangle these resonances and to determine their contributions in the 
different final states, a partial wave analysis has to be performed. The model-independent de- 
termination of all amplitudes can be achieved by a "complete experiment" , for which a set of 8 
well chosen polarization observables is needed in each final state for single pseudoscalar meson 
production[l]. This set consists of at least 4 double polarization observables like the observable 
G, which has been measured with the CBELS A/TAPS experiment in Bonn. In the experiment 
linearly or circularly polarized photons and a longitudinally or transversally polarized target are 
provided, therefore several polarization observables are accessible. 

The first results for the double polarization observable G, which utilizes linearly polarized photons 
on a longitudinally polarized target, are shown here for tt° and r\ photoproduction off the proton. 

1.2 Experimental Setup 

The Crystal Barrel experiment is located in Bonn at the electron accelerator ELSA, which provides 
polarized or unpolarized electrons up to 3.5 GeV. Via bremsstrahlung an unpolarized, linearly 
polarized or circularly polarized photon beam can be produced, whose energies are determined in 
a tagging system. The polarized target is surrounded by the main calorimeter, the Crystal Barrel 
detector, while the forward direction is covered by the Mini- TAPS detector up to small angles. 
Each calorimeter is equipped with detectors for charge identification. 

This setup allows an angular coverage of nearly 47r with a high detection efficiency for photons 
and is therefore well suited for the measurement of reactions with neutral mesons like ir° and i] in 
the final state. 

1.3 Selection of the Reactions 

Two different polarization settings where used for this analysis with coherent edge at 950 McV 
and at 1150 MeV. With these two settings, which are shown in figure 1.1, left, the photon energy 
range of E 1 = 700 MeV— 1100 MeV can be covered with a high degree of photon polarization. 
For the selection of the two photon decay of the r\ and 7r°, first, the reactions were checked for two 
neutral particles and one charged for the recoiling proton. The time information of all particles of 
the reaction were used, to confirm that the observed particles originated from the same reaction. 
Also, the data was corrected for random time background. To ensure that the meson and the 
recoiling proton were detected in opposite directions in the center of mass system, the <f> and 6 
angle difference were analyzed. By treating the proton as a missing particle and plotting the 
missing mass of the system, a peak at the mass of the proton could be observed. In the invariant 
mass of the two photons, which is shown in figure 1.1, center, a clear signal for the 7r° and r\ can 
be seen above nearly no background. To choose the desired reaction, a cut on either meson was 
done. 

1.4 Determination of the Polarization Observables 

In the frozen spin target, butanol is used as target material, which is composed of hydrogen, 
carbon and oxygen. Because of that, reactions on free and on bound protons are contributing. 

*Hclmholtz-Institut fur Strahlcn- und Kcrnphysik, Universitat Bonn, thieiahiskp.uni-bonn.de 
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Figure 1.1: The two different settings of the coherent peaks for linearly polarized photons (left). 
Invariant mass of the two photons in logarithmic scale (center) with the peaks for the 7T° and the 
77 meson. Example fit to the 4>(tt°) distribution of the data (right) to extract the observables. 



The count rate for the reactions on the butanol target can be written as 

N(9, cp) = ( Nb + N f )-[l- ^±^1 • cos(20) + j^j^P* ■ V l ?G sin(2^)] (1.1) 

with Nf(8) (Nb(6)) being the count rate on the free (bound) protons. The observables, which 
can be measured are the beam asymmetries Sfc and £ / and the double polarization observable 
G, which can only be observed on the polarizable hydrogen protons. To successfully determine 
the double polarization observable G and the beam asymmetry £ on the butanol the following 
equation was used to fit the data: 

N{9, <p) = A ■ [1 - B ■ cos(2<£) + C ■ sin(2<£)] (1.2) 

An example fit is shown in figure 1.1, right. 



1.4.1 The Beam Asymmetry £ 

The fit parameter B allows the extraction of the beam asymmetry on the free (£/) and on the 
bound protons (£&)• By reducing the fraction of bound protons in the data set, the measured 
beam asymmetry converges to the beam asymmentry on the free proton and can be compared to 
previous measurements. In figure 1.2 the beam asymmetry for ir° and r\ photoproduction is shown, 
compared to previous measurements and results of the different partial wave analyses. The data 
agrees well with the measurements on the free proton, which indicates that the influence of the 
bound protons could successfully be reduced. 




COS cos u v 



Figure 1.2: The beam asymmetry £ (red) for ir°- (left) and ^-production (right). The values 
are compared with the measurements of the GRAAL collaboration[5] (black) and for the r\ with 
previous CBELSA/TAPS measurements [6] (blue). The lines give the solutions of the partial wave 
analyses: MAID [2] (black), SAID [3] (blue) and BnGa[4] (red). 



1.4.2 The Double Polarization Observable G 

From the fit parameter C the double polarization observable G can be extracted, modified by the 
dilution factor D = N f . To determine this factor, data sets with hydrogen and carbon targets 
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were compared. The comparison of the missing mass for the different materials in fig. 1.3, left, 
shows that the sum of the hydrogen and carbon measurements describes the butanol data. The 
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Figure 1.3: Missing mass distribution (left) and meson cos# (center) for tt° production on the 
different target materials: butanol (black), hydrogen (blue), carbon (red) and the sum of hydrogen 
and carbon (green). The dilution factor (right) was determined using these distributions. 

same applies to the cos# of the meson (fig. 1.3, center). With these distributions, it is possible to 
determine the dilution factor (fig. 1.3, right). This factor was determined for both mesons and all 
energies and is used to correct the values of the observable G. The double polarization observable, 
which is shown in fig. 1.4, was successfully extracted in the energy range for both mesons. 




Figure 1.4: The double polarization observable G for the 7r° (left) and the t\ (right). Colors of the 
predictions of the partial wave analyses as in figure 1.2. 



1.5 Summary 

The determination of polarization observables is crucial for the understanding of the excitation 
spectrum of the nucleon. For the beam asymmetry £ and the double polarization observable G 
data has been taken by the CBELS A/TAPS experiment in Bonn. These data will provide new 
information for the partial wave analysis, which will help to identify the contributing resonances 
to the different final states. Supported by the DFG (SFB/TR16). 
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THE AT* (1535) EXCITATION - 
THE ELECTROMAGNETIC TRANSITION 

FORM FACTORS 

G. RAMALHO* M.T. PEN A*^ and K. TSUSHIMA^ 

The A* (1535) resonance, an S\\ state with J p = A , is a very important nucleon excitation. 
It has a strong decay channel for r/N (« 50%), and has been studied by several models and 
frameworks [1, 2, 3, 4, 5, 6, 7, 8]. A notable propriety of the A* (1535) resonance is the relation 
between the two independent helicity amplitudes for the transverse A 1 / 2 and the longitudinal 
Si/2- Although it has been assumed till recently that the amplitude Si/ 2 was negligible compared 
to A1/2, recent data [9] revealed that these two amplitudes are correlated for the region Q 2 > 2 
GeV 2 , where Q 2 is the negative of the momentum transfer squared. The ratio between the two 
amplitudes for Q 2 > 2 GeV 2 can be expressed as [2] 

S 1/2 (Q)_-— ___A 1/2 (Q ), (1.1) 

i — o 2 
where Q = \fQ 2 , Ms is the TV* (1535) mass, M the nucleon mass and r = ( M ^+m) 2 • ^ ne * es * °^ 

the relation (1.1) using the A-yi 2 data [9, 10] and the model from Ref. [1], is presented in Figure 1.1. 








Q 2 (GeV 2 ) 

Figure 1.1: Si/ 2 helicity amplitude for the 7 A — > iV*(1535) reaction. The function F{Q 2 ) in the 
figure is F{Q 2 ) = — M 2 s m^q ' as defined by Eq. (1.1). Model for A X j 2 is from Ref. [1]. Data 
are from Refs. [9, 10]. 

Although small, the magnitude of Si/ 2 is about 10% of A 1 / 2 for very high Q 2 [2]. The results 
for S1/2 can be interpreted as a correlation between the valence quark degrees of freedom and the 
quark-antiquark excitations (or meson cloud excitations) , in the framework of the spectator quark 
model [1, 2, 11, 12]. 

In general at high Q 2 the valence quark degrees of freedom are dominant, since meson cloud 
effects are expected to be suppressed in QCD according to the dimensional counting rules [8]. At 
low Q 2 , meson cloud effects can be very important particularly for the lightest meson, the pion, 
due to chiral perturbation theory [13]. An an example of the insufficiency of the valence quark 
degrees of freedom alone to explain the experimental data at low Q 2 , is the *yN — > A reaction. 
Including only the valence quark effects it is impossible to describe the data and misses at least 
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30% of the magnetic transition form factor strength for the jN — >• A reaction, in a non-relativistic 
framework near Q 2 = [8, 14, 15]. This is interpreted as a manifestation for the necessity of 
including the explicit meson cloud effects [8, 14, 15, 16, 17, 18, 19]. 

A definite decomposition into the valence quark and meson cloud effects is however impossi- 
ble [13, 20]. The decomposition can be done only for a particular model. 

The application of a valence quark model to describe an electromagnetic transition between 
two baryon states has the advantage of accounting the quark structure (baryon wave functions 
and quark electromagnetic current), and can provide an upper limit for the magnitude of the 
quark core contribution for the electromagnetic form factors. [Recall the case for the — > A 
reaction] . Recent models developed for the treatment of meson-baryon reactions using the meson 
and baryon effective degrees of freedom, take the quark core contributions in consideration [15, 21]. 
In addition the quark model estimates can be tested in the high Q 2 regime, where meson cloud 
effects fall down. Another regime where meson cloud effects can be suppressed and test the quark 
model estimates are the lattice QCD simulations with large pion masses, if some quark models 
can be extended to that regime. 

To study the A r *(1535) system we use the covariant spectator quark model [1, 11, 12, 22]. 
In the model a baryon is described as a 3 constituent-quark system that can be re-arranged in 
a quark-diquark structure, where the quark can couple with photon in terms of a vector meson 
dominance (VMD) process and the diquark is treated as an on-shell (spectator) particle [11, 12]. 
The model has been developed, and calibrated by the nucleon [11] and A [16, 17, 18] systems, 
and later extended to study the other resonances such as the Roper and A(1600) [23, 24]. Using 
the VMD parametrization for the quark current, the model was also extended successfully to 
the lattice QCD regime for > 400 MeV, particularly for the elastic reaction [19] and for the 
transitions from the nucleon to the A(1232) and A r *(1440) systems [18, 19, 23]. The results show 
that the covariant spectator quark model has a good control of the valence quark model degrees 
of freedom. 

In the covariant spectator quark formalism the N* (1535) is described as a quark-diquark system 
with negative parity, a relative orbital angular momentum 1 and a quark core with spin 1/2 [1]. 
In addition we assume that the diquark has no internal P-state structure and can be treated as 
pointlike [1]. Furthermore, we note that the A r *(1535) wave function can only be reliably applied 

for Q 2 > |q|g = 0.23 GeV 2 , where |q| = M |m^ is the photon three-momentum in the A r *(1535) 
rest frame. The results for the transition form factors F* and J 7 ^* are presented in Figure 1.2 for 
Q 2 > 1 GeV 2 . From the figure, one can conclude that the model is reasonable to explain the 
results for Fj* for Q 2 > 1 GeV 2 . As for F£ one can notice that the model results differ in sign 
from the data, and also that it overestimates in absolute values of the data, where the data show 
F 2 * ss for Q 2 > 1.5 GeV 2 . The result F| = is equivalent to the relation given by Eq. (1.1). In 
the spectator quark model the deviation from the data can be interpreted as the missing of the 
meson cloud effects in the calculation [2] . If this interpretation is correct we can expect a negative 
contribution from the meson cloud for F% around Q 2 = 2 GeV 2 with the magnitude similar to the 
present model. 

Although one cannot compare directly the valence quark effects and/or meson cloud effects, 
from different formalisms, it is encouraging to note that our results for F% are close to the estimate 



129 



PrimcNct Workshop 2011, Jiilich 



G. Ramalho 



of the EBAC/Jlab group [6] for the quark core contributions, using a framework based in effective 
baryon-meson interactions. It is also worth to note that a calculation of the form factors with 
the assumption that the iV*(1535) is a dynamically generated resonance, and therefore takes into 
account only the meson cloud effects, gives only negative contributions for _F 2 * [4]- 

Based on the results from different works, one can conclude that valence quark effects are very 
important [4], or are the dominant contribution [1, 5, 6, 7, 8, 25] for the TV* (1535) transition form 
factors, although meson cloud is certainly necessary to describe the data accurately. 
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THE STRUCTURE AROUND W=1680 MeV IN 
77-PHOTOPRODUCTION OFF THE NEUTRON 

DOMINIK WERTHMULLER*^ and LILIAN WITTHAUER* 
for the A2, Crystal Ball and TAPS collaborations 

1.1 Motivation 

Previous experiments of the GRAAL [1, 2], LNS-Scndai [3] and CBELS A/TAPS [4, 5] collabora- 
tions have reported a narrow structure at W w 1680 MeV in the total cross section of quasi-free 
77-photoproduction off the neutron. This structure is not seen on the proton and its nature is up 
to now still unknown. As the width of the structure is close to the experimental resolution of the 
different experiments the true width should be unusually small compared to known resonances 
in this energy region. Therefore the structure can hardly be described using a single well-know 
broad nucleon resonance. Theoretical explanations thus include interferences of resonances [6] or 
contributions from intermediate meson loops [7] . Another solution is provided by the chiral quark 
soliton model [8] which predicts an anti-decuplet nucleon resonance that couples strongly to the 
formation channel 7n as well as to the final state rjN. 

We report preliminary results from two different experiments performed at the tagged-photon 
beam facility at MAMI. In one experiment a liquid deuterium target was used as a neutron target 
whereas in the other one a liquid 3 He target was installed. 




Tagger Crystal Ball TAPS 

Figure 1.1: Scheme of the experimental setup. 



1.2 Experimental setup 

The measurements were performed at the electron accelerator facility MAMI [9, 10] in Mainz. A 
schematic view of the experimental setup is shown in figure 1.1. A photon beam was produced 
from an electron beam of 1.5 GeV energy via bremsstrahlung using a copper radiator. The energy 
of the photons was determined by a momentum analysis of the scattered electron in a magnetic 
spectrometer (Glasgow photon tagger [11, 12, 13]). After collimation the beam impinged on the 
target (liquid deuterium or 3 He, respectively). The target was surrounded by a cylindrical plastic 
scintillator strip detector [14], which was used for charged particle identification, and the spherical 
electromagnetic calorimeter Crystal Ball [15]. This detector consists of 672 Nal crystals and 
covers 94% of 47r steradians. The hole in forward direction of Crystal Ball is closed by the TAPS 
detector [16, 17] which is made of 384 BaF 2 crystals (378 BaF 2 + 24 PbW0 4 crystals in the 3 Hc 
experiment). In front of every crystal a thin plastic scintillator element is installed as a charged 
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particle veto detector. As trigger condition a deposited energy sum of 300 MeV in the Crystal 
Ball and a total multiplicity of two or more hits in both calorimeters was requested. 

1.3 Data analysis 

Using the rj — > 77 decay of the 77-meson the final state products of the quasi-free reactions 7p— > r)~p 
and 7n— > ?m were selected by requesting 2 neutral clusters + 1 charged cluster or 3 neutral clusters, 
respectively, in the detector system. The ?y-meson was identified in the invariant mass of the two 
photons by applying a cut around the rj mass. Background coming from e.g. 777r-photoproduction 
was suppressed by cutting on the r/iV-coplanarity and the missing mass. 
Excitation functions were calculated as functions of 

1. the incoming photon beam energy E-y (converted for better comparison to the center-of-mass 
energy Wbeam = y/2E 1 m,N + m? N of the photon-nucleon system) 

2. the invariant mass of the final state ^-nucleon system W rec = M(r]N) by reconstructing the 
energy of the recoil nucleon by kinematics 

In case of 1) structures are smeared by the Fermi motion of the initial state nucleon inside the 
target nucleus (deuteron or 3 He). In contrast structures appearing in 2) are not affected by Fermi 
smearing and the resolution only depends on the energy and angular resolution of the detector 
system. 

The multiple hits in the photon tagger were treated as single events and weighted accord- 
ingly to the tagger-detector timing to achieve a statistical subtraction of the accidental coincident 
background. 

Finally, cross sections were extracted by normalizing the excitation functions with the target 
density, the photon flux and the detector acceptance obtained by a Geant4 based simulation of 
the detector setup. 




Figure 1.2: Very preliminary excitation functions for the quasi- free 7p— > f?p and 711—)- rm reactions 
obtained by the LD2 and the 3 He experiments. 



1.4 Preliminary results 

Preliminary results of the analysis are shown in figure 1.2. Since the absolute normalization of the 
data is not final all results were scaled to a n in the Sn(1535)-maximum of the deuterium analysis 
and are given in arbitrary units. 

The previous results are confirmed by both experiments: An excess in the neutron cross section 
is observed which is not seen on the proton. The width of the structure in a(W rec ), that is not 
affected by Fermi motion, is around 50 MeV and is mainly resolution dominated. 
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Furthermore the results of the LD2 and the 3 He measurement are in good agreement. The 
differences in a(Wbeam) seem to be only related to the larger Fermi motion inside the 3 He nucleus 
(smearing at threshold and at the shoulder of the Sn(1535)). The deviation in a(W rec ) is due 
to the fact that the kinematic reconstruction within the participant-spectator model is only an 
approximation in case of 3 Hc. Obtaining consistent results by using nuclei with different Fermi 
motion and n/p-ratio makes it very unlikely that the observed effect is caused by rescattering of 
mesons or FSI. 

The final results of the LD 2 measurement will include the analysis of the r\ — > 3ir° decay 
channel and will provide angular distributions with high statistics. 
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A THEORETICAL APPROACH 
TO rj'N SCATTERING 

A. RAMOS*and E. OSET* 



1.1 Introduction 

The rf meson has interesting properties associated to the underlying QCD dynamics of hadrons, in 
particular to the U(1)a axial vector anomaly [1]. Being close to a singlet of SU(3), its interaction 
with nucleons is supposed to be weak compared for instance with the case of its partner, the 77 
meson. Experimentally, there are only poor estimations of the rf N scattering length from analyses 
of the pp — > pprf cross section near threshold at COSY [2], establishing a value around |o^'jv| ~ 
0.1 fm. The rj'N interaction has also an infuence on photoproduction -fp — > rj'p reactions measured 
at ELSA and Jlab [3], which have been analyzed so far in the framework of resonant models [4]. 
Finally, some information on the rj'N interaction can be extracted from the transparency ratio in 
(7,77') reactions on nuclei [5]. 

We present a theoretical study of the rj'N interaction within a unitary approach in coupled 
channels [6]. Our study also includes the explicit effect of resonances lying very close to the rj'N 
threshold, which are generated dynamically from the vector meson-baryon (VB) interaction in 
s-wave [7] . Finally, we also consider a term in the Lagrangian coupling the baryons to the singlet 
component of the pseudoscalar meson, which is allowed by the symmetries of QCD [8]. We find 
that the elastic rfN cross section is very sensitive to the unknown singlet strength, while the 
inelastic cross sections are rather stable and constitute a genuine prediction of our model. 

1.2 Model 

The meson-baryon scattering dynamics is taken from the lowest order chiral Lagrangian reduced 
to the two meson fields needed in the process, 

= (Bi7"^[($0„$ - 0„$$)B - B($a„* - (1.1) 

where / is the pion decay constant, B represents the baryon octet matrix and $ contains the 
pseudoscalar meson octet as well as the singlet field rji/y/3 added to its diagonal elements. The 
physical 77,77' mesons are given by 77 = cos Op 773 — sin Op rj\ and rf = sin Op 77s + cos #p 771, where 
Op is the rji — rj$ mixing angle, for which we take the value Op = —14.34° reported in [9]. Our 
basis of pseudoscalar meson-baryon states (PB) is ttN, rjN, rj'N, KA and KY,. 

We have also implemented in our model the effect of a resonance appearing around 1970 MeV, 
close to the rj'N threshold at y/s = 1896 MeV, which is generated dynamically from the vector 
meson-baryon (VB) interaction in s-wave [7] and couples mostly to the K*A and K*H channels. 
As described in Ref. [6], this is achieved by still working within a basis of coupled PB states 
but including explicitly, as part of a PB — PB potential, the transition of the PB to the VB 
channels, the later ones interacting among themselves to produce the resonance. The PB — VB 
transition vertex contains the standard PPV couplings as well as the anomalous PVV terms since 
the standard ones are supressed by sin 6>p for transitions involving rj'N states. 

Finally, we also include a Lagrangian that couples the singlet component of the meson field to 
the baryons [8] and gives rise to new contributions to the transition potentials between rjN and 
rj'N states, V$ nN = Csm 2 P , V$ V , N = -Csm0 P cos0 P and Vfy >v , N = Ccos 2 6» P , where the 
unknown singlet strength C is written in terms of a parameter a obtained from C = a 2m v > (Eb + 
E B ,)/Afl/2M N . 
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With the obtained meson-baryon transition potential V we derive the scattering amplitude T 
by solving the Bethe-Salpeter equations in their on-shell factorization form, T = [1 — VG]~ X V, 
where G is the loop function for the intermediate meson-baryon states which we evaluate using 
dimensional regularization. 



1.3 Results 

Our results for the rj'p scattering length are shown in Table 1.1, while the elastic and inelastic rj'p 
cross sections, together with the cross section of the reaction n~p — > rj'n, are displayed in Fig. 1.1. 



a 


model 


a vp [fm] 


K'p| [H 





PB 


0.0017+10.0139 


0.014 


PB- 


-VB 


0.0210+10.0192 


0.029 


-0.126 
0.204 


PB- 
PB- 


-VB 
-VB 


0.073 + i0.019 
-0.072 + 10.020 


0.075 


-0.193 


PB- 


-VB 


0.098 + i0.020 


0.1 


0.256 


PB- 


-VB 


-0.098 + 10.020 


-0.333 


PB + VB 


0.149 + i0.020 


0.15 


0.352 


PB + VB 


-0.149 + 10.021 



Table 1.1: The rj'p scattering lengths for various models discussed in the text. 



The rj'p scattering length is very small when only the Weinberg- Tomozawa terms are included 
(a = 0, PB), since the particular structure of the Lagrangian of Eq. (1.1), $9^$ — makes 
the contribution of the singlet vanish. Including the VB states (a — 0,PB + VB) increases 
the real part by an order of magnitude an enhances the imaginary part considerably due to the 
decay channels K*A — > nKA and K*T, — > ttKY,. The effect of the VB channels is more visible 
by comparing the thin and solid lines in Fig. 1.1. In all cases, one sees a pronounced structure 
around p v i = 500 MeV/c (p^ = 1600 MeV/c) which corresponds to the energy 1970 MeV of 
the resonance. The strongest relative enhancements are seen for the elastic channel and for the 
tt~P — > rj'n reaction, but the results do not reproduce the available experimental information. The 
scattering length only amounts to \a v ' p \ — 0.029 fm, while the n~p — » rj'n cross section is six times 
smaller than the measured value of around 0.1 mb at a pion lab momentum of 1600 MeV/c [10]. 

Important changes are obtained when the singlet Lagrangian is considered with its strength a 
tuned to values of |ay p | around 0.1 fm. Each assumed value of |a^/ p | admits two types of solutions 
for a, one positive (repulsive) and one negative (attractive). For a given sign, the changes in a 
modify the real part of a^/p in a nearly proportional way. The imaginary part remains practically 
constant because the singlet Lagrangian induces the largest effect on the rj'N — > rj'N amplitude 
as it behaves like cos 2 #p. This also means that our model, in spite of using the rj'N scattering 
length as a free parameter, provides a robust prediction for the inelastic channels, as corroborated 
by the moderate changes seen in the inelastic cross sections displayed in Fig. 1.1. 

The optical theorem allows us to obtain the total rj'p cross section from the elastic scattering 
amplitude and, by subtraction, we can derive the inelastic one, a mc = a tot — er cl , which amounts to 
5 mb at a value of 200 MeV/c for the rf momentum. The preliminar analysis of the transparency 
ratio of rf photoproduction in nuclei [5] estimate the inelastic cross section to be around 10 mb 
for a rf momentum of 1 GeV/c. This is an upper bound for the one-body absorption inelastici- 
ties, since the transparency ratio can also be contaminated by multinucleon absorption processes, 
although they are expected to be small [11]. In any case, the rj' energies of the transparency 
ratio experiment are too high for our model, the validity of which we trust for momenta smaller 
than p^ — 600 MeV/c, implying rf kinetic energies smaller than 200 MeV. Data analyses with 
momentum cuts will be available in the future and, although the statistics will be lowered, we 
expect them to help in constraining the properties of the rfN interaction. 
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Figure 1.1: Elastic and inelastic rj'p cross sections, together with the cross section of the reaction 
7r~p — > rj'n, for various models discussed in the text. 
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1.1 Introduction 

The investigation of the exotic objects in the nuclear physics is a proven method for revealing many 
interesting properties of nuclear systems and for accessing to an unexplored areas of physics. The 
recent progress in the spectroscopy of pionic and kaonic atoms, as well as pionic and kaonic nuclei 
has permitted to obtain deeper insights into the meson-nucleus interaction and the in-medium 
behaviour of spontaneous chiral symmetry breaking [1]. 

Analogically to the other exotic nuclear systems, the investigation of the 77-mesic nuclei would 
provide many interesting informations about the 77-N interaction, N* in-medium properties [2] and 
would deepen our knowledge of the fundamental structure of the 77 meson [3]. The 77 meson is 
electrically neutral, therefore such a system can be formed only via the strong interaction which 
distinguishes it qualitatively from pionic atoms where the binding is the effect of the sum of the 
attractive electromagnetic force and the repulsive strong interaction. 

The search of the 77 - mesic nucleus was performed in many experiments in the past [4, 5, 6, 7, 
8, 9] and is being continued at COSY [12, 10, 11], JINR [6], J-PARC [13] and MAMI [9]. Many 
promising indications where reported, however, so far there is no direct experimental confirmation 
of the existence of mesic nucleus. In the region of the light nuclei systems such as ?7-He or 77- 
T, the observation of the strong enhancement in the total cross-section and the phase variation 
in the close-to-threshold region provided strong evidence to the hypothesis of the existence of a 
pole in the scattering matrix which can correspond to the bound state [14]. However, as it was 
stated in [15, 16], the theoretical predictions of width and binding energy of the 77-mesic nuclei is 
strongly dependent on the not well known subtreshold 77-nucleon interaction. Therefore, the direct 
measurements which could confirm the existence of the bound state, are mandatory. 

1.2 Experiment 

In June 2008 we performed a search for the 4 He — 77 bound state by measuring the excitation 
function of the dd — > 3 Hep7r~ reaction near the 77 meson production threshold using the WASA- 
at-COSY detector [17]. During the experimental run the momentum of the deuteron beam was 
varied continuously within each acceleration cycle from 2.185 GeV/c to 2.400 GeV/c, crossing the 
kinematic threshold for the 77 production in the dd — > 4 HeT7 reaction at 2.336 GeV/c. This range of 
beam momenta corresponds to the variation of 4 Hc — 77 excess energy from -51.4 MeV to 22 MeV. 
The experimental method is based on measuring the excitation function for dd — > 3 Hep7r~ and 
a search for a resonance-like structure below the 4 He — 77 threshold. The relative angle between 
the outgoing p — tt~ pair which originates from the decay of the N*(1535) resonance created via 
absorption of the 77 meson on a nuclcon in the 4 He nucleus, is 180° in the N* reference frame 
and is smeared by about 30° in the reaction center-of-mass frame (CM) due to the Fermi motion 
of nucleons inside the 4 He — 77 nucleus. The 3 He plays the role of a spectator and therefore 
its momentum in the CM frame is relatively low and determined only by the Fermi momentum 
distribution in the 4 He — 77. This signature distinguishes production of mesic helium from other 
reactions, in which the distribution of the 3 He momentum reaches much higher values. 

Figure 1.1 presents the not normalized excitation functions for a "signal-rich" region: p^ < 
0.3 GeV/c (top) and a "signal-poor" region: p^ 1 ^ > 0.3 GeV/c in 20 intervals in beam momentum 
(middle). Also, the difference between two regions is shown in the same Figure (bottom). 
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Figure 1.1: Excitation function for the dd — > 3 Hep7r~ reaction for the "signal-poor" region 3 He 
momentum above 0.3 GeV/c (top plot) and the " signal- rich" region, low 3 He momentum below 
0.3 GeV/c (middle plot). The excitation functions are not normalized. The beam momentum 
corresponding to the He — 77 kinematical threshold is marked as a red dashed line. The difference 
between the "signal-rich" and the "signal-poor" regions is shown at the bottom plot. 



The figure indicates no structure below the kinematical threshold where the signal is expected. 
In order to estimate the maximum upper limit for the cross-section for the production of the bound 
state via dd — > 4 He — r] bound — > 3 Hep7r~ reaction, the Breit-Wigner function along with the linear 
background was fitted to the normalized " signal- rich" region. The absolute normalization was 
obtained based on the dd — > 3 Hen reaction, whereas the luminosity dependence of the beam 
momentum was determined using the quasi-elastic dd — > pp(n sp n sp ) reaction. The integrated 
luminosity equals L = 117.9 ± 13.6 nb _1 . 



1.3 Outlook 

In November 2010 a new two-week measurement was performed with WASA-at-COSY. We col- 
lected the data with approximately 20 times higher statistics. In addition to the dd — > 3 Hep7T~ 
channel we registered also the dd — > 3 Hen7r° reaction. The data analysis is undergoing. After two 
weeks of measurement with an estimated luminosity of 4 • 10 30 cm -2 s _1 , we expect a statistical 
sensitivity of a few nb (a). A non-observation of this signal will significantly lower the upper limit 
for the existence of the bound state. 



141 



PrimeNet Workshop 2011, Jiilich 



W. Krzemieri 




-50 -40 -30 -20 -10 10 20 

Q [MeV] 



Figure 1.2: Normalized excitation function for the dd — > 3 Hep7r reaction for the "signal-rich" 

A-(-) 2 

region. The Breit-Wigner function f(Q) = -rp, — - \l — rrro along with the linear function was 

fitted. The parameters T and Ebe were fixed to the value of -20 MeV and 10 MeV respectively. 
The free parameter A is consistent with zero within the obtained accuracy. The dotted red line 
represents the linear background. 



1.4 Support 

This work has been supported by FFE funds of Forschungszetrum Juelich, grant No 41831803 
(COSY-107), by the European Commission under the 7th Framework Programme through the 'Re- 
search Infrastructures' action of the 'Capacities' Programme. Call: FP7-INFRASTRUCTURES- 
2008-1, Grant Agreement N. 227431 and by the Polish Ministry of Science and Higher Education 
under grants No. 2367/B/H03/2009/37 and 0320/B/H03/2011/40. 

Bibliography 

[I] S. Hirenzaki, Prog. Theor. Phys. Suppl. 168 (2007) 458-465. 

[2] D. Jido, H. Nagahiro, S. Hirenzaki, Phys. Rev. C66 (2002) 045202. 

[3] S. D. Bass, A. W. Thomas, Acta. Phys. Pol. B 41 (2010) 2239. 

[4] B. J. Lieb et al, Proc. Int. Nucl. Phys. Conf., Sao Paulo, Brazil (1989). 

[5] G. A. Sokol et at, arXivmucl- ex/ '9905006 (1999) 

[6] M. Kh. Anikina et al, arXiv.nucl-ex/0^12036 (2004). 

[7] A. Gillitzer, Acta Phys. Slovaca 56 (2006) 269. 

[8] A. Budzanowski et al, Phys. Rev. C79 (2009) 061001(R). 

[9] B. Krusche, F. Pheron, Y. Magrbhi Acta. Phys. Pol. B 41 (2010) 2249. 

[10] J. Smyrski et al, Phys. Lett. B 649, 258 (2007). 

[II] T. Mersmann et at, Phys. Rev. Lett. 98, 242301 (2007). 
[12] P. Moskal, J. Smyrski, Acta. Phys. Pol. B 41 (2010) 2281. 
[13] H. Fujioka, K. Itahashi, Acta. Phys. Pol. B 41 (2010) 2261. 
[14] C. Wilkin Acta. Phys. Pol. B 41 (2010) 2191. 

[15] Q. Haider, L.C. Liu, Phys. Rev. C66 (2002) 045208. 

[16] Q. Haider, L.C. Liu, Acta Phys. Polon. B 41 (2010). 

[17] W. Krzemieh et al, Acta Phys.Polon.Supp. 2 (2009) 141-148. 



142 



List of Participants 



Name 


Affiliation 


Contribution 


Adlarson, Patrick 


Uppsala University 


Page 21 


Amaryan, Moskov 


Old Dominion University 


Page 80 


Bashkanov, Mikhail 


Tubingen University 


Page 96 


Beck, Rcinhard 


HISKP, Bonn University 




Bergmann, Florian 


Westfalische Wilhelms Universitat Munster 


Page 11 


Berlowski, Marcin 


Soltan Institute Nuclear Studies 


Page 9 


Bijncns, Johan 


Lund University 


Page 68 


Caldeira Balkestahl, Li 


Uppsala University 


Page 18 


Coderre, Daniel 


Forschungszentrum Jiilich 


Page 77 


Czerwinski, Eryk 


Jagiellonian University Cracow 


Page 94 


Di Donato, Camilla 


I.N.F.N. Naples 


Page 35 


Ditsche, Christoph 


HISKP, Bonn University 




Dybczak, Adrian 


Jagiellonian University Cracow 


Page 111 


Fariborz, Amir 


State University of New York, Institute of Technology 


Page 74 


Frohlich, Ingo 


Goethe University Frankfurt 


Page 16 


Gauzzi, Paolo 


Sapienza Universita di Roma / INFN 




Giovannella, Simona 


I.N.F.N. Frascati 




Goldcnbaum, Frank 


Forschungszentrum Jiilich 




Goslawski, Paul 


Westfalische Wilhelms Universitat Munster 


Page 91 


Grzonka, Dieter 


Forschungszentrum Jiilich 




Gullstrom, Carl-Oskar 


Uppsala University 


Page 43 


Gumberidze, Malgorzata 


Institut dc Physique Nucleaire d Orsay 


Page 106 


Hanhart, Christoph 


Forschungszentrum Jiilich 




Heijkcnskjold, Lena 


Uppsala University 


Page 61 


Hcjny, Volker 


Forschungszentrum Jiilich 




Hodana, Malgorzata 


Jagiellonian University Cracow 


Page 45 


Hoferichter, Martin 


HISKP, Bonn University 


Page 122 


Hoistad, Bo 


Uppsala University 




Johansson, Tord 


Uppsala University 




Kampf, Karol 


Lund University 




Khan, Farha Anjum 


Forschungszentrum Jiilich 


Page 52 


Khoukaz, Alfons 


Westfalische Wilhelms Universitat Munster 




Klaja, Joanna 


Forschungszentrum Jiilich 




Kolesar, Marian 


Charles University Prague 


Page 30 


Krzemien, Wojciech 


Jagiellonian University Cracow 


Page 140 


Kubis, Bastian 


HISKP, Bonn University 


Page 13 


Kuc, Hubert 


Jagiellonian University Cracow 


Page 117 


Kupsc, Andrzej 


Uppsala University 




Lanz, Stefan 


Bern University 


Page 24 


Lersch, Daniel 


Forschungszentrum Jiilich 


Page 38 


Leupold, Stefan 


Uppsala University 




Machner, Hartmut 


Universitat Duisburg-Essen 




Masjuan, Pere 


Universidad de Granada 


Page 71 


Messchendorp, Johan 


KVI / RU Gronigen 


Page 83 


Metag, Volker 


Universitat Giessen 


Page 51 


Moskal, Pawel 


Jagiellonian Uniwersity Cracow 




Niecknig, Franz 


HISKP, Bonn University 


Page 63 


Nikolaev, Alexander 


HISKP, Bonn University 


Page 90 


Nuhn, Patrick 


HISKP, Bonn University 




Ozcrianska, Iryna 


Jagiellonian University Cracow 


Page 87 


Perez del Rio, Elena 


Tubingen University 


Page 100 


Prado Longhi, Ivan 


Universita degli Studi "Roma Tre" and I.N.F.N. 


Page 119 


Pricking, Annette 


Tubingen University 




Pszczel, Damian 


Andrzej Soltan Institute for Nuclear Studies Warsaw 





Continued on next page 



143 



List of Participants - continued from previous page 



Name 


Affiliation 


Contribution 


Ramos 0617107; Antrels 


Universitat He Barcelona 


Paee 137 


Ramalho, Gilberto 


Instituto Superior Tecnico Lisbon 


Page 128 


Ramstcin, Beatrice 


Institut dc Physique Nucleaire d Orsay 




Rcdmcr, Christoph Florian 


Uppsala University 


Page 68 


Pvitman, James 


Ruhr-Universitat Bochum and Forschungszcntrum Jiilich 




Salabura, Piotr 


Jagiellonian University Cracow 


Page 114 


Sawant, Siddesh 


Indian Institute of Technology Bombay 


Page 61 


Schadmand, Susan 


Forschungszcntrum Jiilich 




Schneider, Sebastian 


HISKP, Bonn University 


Page 63 


Sibirtsev, Alexander 


Forschungszentrum Jiilich 




Skorodko, Tatiana 


Tubingen University 


Page 103 


Stollenwerk Felix 


Forscrinn0'S7,ontrnm Julicb 


Prljrp 33 


Teilab Khaled 


CrOotho Tlnivorsitv 'Frankfurt 




Terschliisen, Carla 


Uppsala University 


Page 58, Page 40 


Thiol, Annika 


HISKP, Bonn University 


Page 125 


van Pee, Harald 


HISKP, Bonn University 




Werthmullcr, Dominik 


University of Basel 


Page 131 


Wilkin, Colin 


University College London 




Wirzba, Andreas 


Forschungszentrum Jiilich 


Page 5 


Witthauer, Lilian 


University of Basel 




Wurm, Patrick 


Forschungszentrum Jiilich 


Page 48 


Zdebik, Jaroslaw 


Jagiellonian University Cracow 


Page 54 


Zdrahal, Martin 


IPNP, MFF, Charles University Prague 


Page 27 



144 



